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Abstract 



This thesis is about the use of a novel, exact functional quantization method as 
applied to two commonly studied actions in theoretical physics. The functional 
method in question has its roots in the exact renormalisation group flow techniques 
pioneered by Wilson, but with the flow parameter not limited to the familiar 
momentum cutoff. Finding a configuration satisfying an expression for the exact 
effective action which does not vary with this parameter provides the basis for 
finding solutions to the physical actions we study. 

Firstly, the method is applied to an expression for the bare action of the pseudo- 
scalar axion used to explain the strong CP problem in QCD. When quantized, we 
find that the effective potential of the axion, when interactions are not considered, 
is necessarily flattened by spinodal instability effects. We regard this flattening as 
representing the very early stage in the development of the axion potential, when 
the Peccei-Quinn U(l) symmetry is spontaneously broken resulting in a double- well 
potential. Using commonly quoted values for the parameters of such a potential, 
we devise an expression for the energy density of the emerging axion potential and 
this is compared to dark energy. 

We then apply the functional method to the bosonic string with time varying 
graviton, dilaton and antisymmetric tensor (resulting in the string-axion) back- 
ground fields. We achieve a demonstration of conformal invariance in a non- 
perturbative manner in the beta functions, contrasting with conventional string 
cosmology where cancellation of a perturbative expansion is performed. We then 
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offer some hints as to possible cosmological implications of our configuration in 
terms of optical anisotropy. 

The research is closely related to the following three published papers. 

• Non-perturbative string backgrounds and axion induced optical 
activity, J. Alexandre, N. E. Mavromatos and D. Tanner; published in New 
Journal of Physics 10, 2008, hep-th/0708.1154 

• Antisymmetric-tensor and electromagnetic effects in an a'-non- 
perturbative four-dimensional string cosmology, J. Alexandre, N. E. 
Mavromatos and D. Tanner; published in Phys. Rev. D 78, 2008, hep- 
th/0804.2353 

• Quantization leading to a natural flattening of the axion potential, 

J. Alexandre and D. Tanner; published in Phys. Rev. D 82, 2010, hep- 
th/1003.6049 
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Chapter 1 



Introduction 



1.1 Motivation and Context 



This thesis is made up of two parts, one dealing with the behaviour of the QCD- 
axion under quantization, and the other dealing with the axion evident in the 
bosonic string through the antisymmetric tensor field. The uniting feature of 
both components of the thesis is the use of a functional method to derive non- 
perturbative evolution expressions for the effective actions: of the QCD-axion; and 
the bosonic string. Deriving expressions which may be compared to experimental 
results in particle physics typically involves commencing with a suitable bare 
action describing the theory; and then quantizing this (e.g. using path integral 
methods) and identifying and managing the divergences using a regularisation 
and/or renormalisation scheme. The result is usually expressed in a perturbative 
expansion with terms defined by increasing powers of a small parameter (such as 
h or a small coupling). In a useful effective theory, the higher order expansion 
terms in the series may be neglected in the low energy regime. However the use 
of perturbative methods may be limited in certain theories or energy regimes. 
For example, it is impossible to describe confinement in QCD with perturbative 
expansions, which are valid only at high energies. For this reason, it is useful to 
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investigate the the infrared regime of a theory, using a non-perturbative approach 
such as the Exact Renormalisation Group (ERG). The so-called Wegner- Houghton 
equation, derived in ERG theory, is built on Wilsonian renormalisation group 
methods where the behaviour of the parameters of the effective theory with energy 
level k is studied, where k is some energy scale less than the high energy cutoff of 



the theory (termed A in this thesis, Section (2.3)). The renormalisation group is 
the set of transformations k — > k + 5k in which the full quantum theory is described 
as k — > and the higher energy field degrees of freedom are integrated out leaving 
any cutoff dependence in the infrared regime. The Wegner-Houghton approach 
takes the limit ^ — > and results in an exact expression for the evolution of the 
effective action with k, at constant field configuration, which we derive in Appendix 



A. The method used in this thesis is outlined in Section (2.3.3). A key difference 
between our method and the Wegner-Houghton approach is that that the evolution 
parameter is not the energy scale k, but another parameter of the theory which 
when varied can be used to describe the evolution of the theory from classical (bare) 
to fully quantized, a feature which is desirable for example in string theory which is 
assumed to be scale invariant. We apply this new, non-perturbative method to two 
well known bare actions, both of which are outside the standard model and both 
of which may be important at high energy scales, potentially near the Planck scale. 



The QCD-axion was postulated by Peccei-Quinn, [33] and the theory subsequently 
developed into one which dealt with the CP problem. As a scalar field 
which acquires mass, it may also offer a contribution to the theorized existence 
of cold dark matter in the universe, [26] and indeed dark energy, (see for 
example, [IS]). A scalar field which solves the CP problem and also provides 
a meaningful contribution to the missing mass-energy balance in the observed 
universe would indeed be an elegant addition to the standard model should the 
axion be successfully detected. While much theoretical work has been done on 
the postulated axion, this elegance provides a motivation, in this thesis, for 
applying the novel, non-perturbative techniques used in shedding more light on 
its behaviour. 



1.2 Structure of the Thesis 
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String theory is the main candidate for unifying quantum field theory with 
general relativity. String cosmology is the application of string-based models 
to cosmological phenomena. String cosmology models with D = 4 spacetime 
dimensions are particularly useful in application to the observed physical world. 
A popular application of string cosmology models is in determining any preferred 
spacetime direction inherent in the universe, or anisotropy. Lastly, the energy 
scales of string theory and string cosmology are necessarily high and there is a 
question as to whether perturbative approaches are appropriate. This is a key 
driver for the use of our non-perturbative, exact functional method. Motivated by 
these themes, we apply our non-perturbative techniques to create a D = 4 model 
and investigate the role of the string-axion in anisotropic effects. 



1.2 Structure of the Thesis 

This thesis is divided into two parts, one dealing with a study of the QCD axion 
and the other with the bosonic string. (We note that the scalar field h arising 
from the anti-symmetric tensor in the bosonic string is often referred to as the 



string axion (see eq. (5.70). We also refer to this string-axion in this thesis as 
it is a key parameter in our results. We stress however that the QCD-axion and 
the string-axion in the two different sections of this thesis are not related in our 
research) . 

The two sections are united as the treatment of both is based on techniques 
originating in, but distinct from, the exact renormalisation group formalism 
pioneered by Wilson and others in the 1970s. In Chapter ^ a review of some 
aspects of theoretical physics shared by both sections of the thesis is undertaken. 
Effective quantum field theory concepts are introduced along with the effective 
potential. This naturally leads to a review of the Wilsonian formalism for 



renormalisation, the exact renormalisation group methods and in (2.3.3) we outline 
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the precise effective field theory manipulations which are the basis of the non- 
perturbative approach used in this thesis. These non-perturbative techniques were 
made extensive use of in work by my thesis supervisor Dr. Jean Alexandre and my 
second supervisor Prof. Nick E. Mavromatos in pH], [103], pH], [105], pH], [TUT] . 
The classical phenomenon of spontaneous symmetry breaking, which underpins 
many key particle physics theories, including axion theory, is introduced, as is the 
spinodal instability. 

In Chapter Q a review of topics specific to QCD-axion physics is presented. A 
summary of QCD from the viewpoint of the symmetries of the theory is presented. 
The concept of instantons in four dimensional non-Abelian gauge theories is 
covered, which naturally leads to a description of the charge-parity (CP) problem. 
The leading candidate for its resolution is the axion, postulated by Peccei and 
Quinn in 1977, |33j. Its development, the variety of potential axion models and 
axion phenomenology are presented. We should note here an excellent summary 
and history of axion physics from Kim, one of the world's foremost experts, [27] , 
which was relied on heavily for this section. 

Chapter Q presents our treatment of the axion. We initially build a justification 
for the quantisation of the axion and state the relevant actions and potentials we 
rely on. We then provide a detailed account of our computations leading to a 
non-perturbative evolution equation and show how the spinodal instability in the 
non-interacting theory leads to our formulation of the flattened effective potential 
of the axion. The work conducted in this chapter was conducted by myself and 
Dr. Jean Alexandre and published in 2010 in [76J while analysis of the /i = and 



A = cases and comparison with accepted dark energy values in section (4.4) was 
largely my own and is unpublished. 

Chapter ^ provides an overview of bosonic string theory with a focus on conformal 
invariance and the conditions for proving this. Two excellent and widely known 
sources from Polchinski, and from Zwiebach, [77], [78] were invaluable in this. 



1.2 Structure of the Thesis 
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A brief summary of critical string cosmology is presented with mention of non- 
critical and tachyon string cosmology noted. A brief section on optical anisotropy 
is presented as an area where our work on the bosonic string may be relevant. 



Chapter ^ presents our work on the bosonic string. We introduce the action 
we will be using and derive, using non-perturbative methods outlined in section 
(2.3.3) and exact equation for the effective action of string. We then show our 
solution to this exhibits conformal invariance non-perturbatively (our main result 
for this section) and offer some insights into its cosmology. The work in this 
chapter is based on work led by Dr. Jean Alexandre and published in 2007 and 
2008 in |106j . |107j . My solo contribution to this part of the thesis consisted of 
computations related to showing conformal invariance presented in sections (6.1.4) 
and fl6~3~2j). 



Chapter 2 

Topics in Quantum Field Theory 



The two parts of this thesis rely on some common themes and calculation methods 
within classical and quantum field theory. These are outlined here. 



2.1 The Effective Action and Potential 

In applying quantum field theory to a bare action, one of the most common 
regularization processes applies a high energy cutoff, A, above which the validity 
of the theory is not known (other methods include dimensional regularization). 
The resulting expression may contain divergences which need to be isolated and 
removed prior to arriving at a useful effective action at some energy scale fco where 
phenomenology can occur. In quantizing a scalar field theory such as </> 4 theory 



introduced in (2.33), we seek a function which when minimized provides an exact 
expression for the expectation value of the field, < >, with quantum effects 
accounted for. We also want an expression which agrees with the the classically- 
derived equilibrium of to lowest perturbative order. We define the generating 
functional, Z[J] which defines the full quantum theory up to energy scale A in 
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path integral form in Euclidean D dimensional spacetime: 

Z\J\ K = J D0exp- (s[<j>] + J d D xJ(x)<p(x^j , (2.1) 

Here, the Euclidean action is S[cf>] = f d D x(\(d(f)) 2 + U ((/>)) and J(x) is the source 
term which interacts with 0. Using Z[J] we can compute all correlation functions 
for processes by taking taking functional derivatives of Z[J] with respect to the 
source J(x) at the required space-time points. We further define a functional W[J] 
such that (now dropping the A subscript for brevity): 

Z[J] = exp(-W[J]). (2.2) 

where W and Z are both functionals of J the external source. Successive functional 
derivatives of W with respect to J generate the connected correlation functions 
(as opposed to Z which returns both connected and disconnected correlation 
functions). Thus W(J) is more useful in our motivation for deriving measurable 
quantities in a quantum field theory. Taking the derivative of W(J) (noting 
—W(J) = lnZ(J)), with respect to the source J and using jj = j^- x |j gives: 

sw 

5J{x 1 ) 



^ J P00(a;i)exp (-S[4>(x)\ + J d D x </>(x)J(x) 

= <f> d ( Xl ) =< (f>( Xl ) > . (2.3) 



This expression for C ; is the expectation value of the quantum field, < (f)(x) >. 
With source set to zero(that is a theory with only self interactions) we can take a 
second functional derivative of W[J], resulting in two terms. 



5 2 W 1 



sj(x 1 )sj(x 2 y J -° z 

-^J VijHftxJexp (-S[<i>(x)] + J d D x ( j ) (x)J(x) 



V(j)(j)(xi) 4>(x 2 ) exp -S[(f)(x)] + / d x(j)(x)J(x) 



This gives: 



P00(x 2 )exp ^-S[^(x)} + J d D x(P(x)J(x)^ . (2.4) 



5 2 W 

|j=o =< 4>{xi)(t>{x2) > -</>d(xi)<i>d(x2), (2.5) 



5J(xi)SJ(x2) 



2.1 The Effective Action and Potential 



19 



where it is noted that < 4>(xi)(j)(x 2 ) > is the two point correlation function defined 
as the amplitude of propagation of an excitation (or particle) between x\ and x 2 . 



fV4,^ Xl )^x 2 )exp- (jgj + fd D x.J(x)<p(x)) 

< ^ {Xl) >S JP»exp-(SM + Jd°»J(»W»)) ' <2 - 6> 

We consider the contributions to the expression < 0(xi)0(x 2 ) >. It is made up 
of two contributions. Firstly the sum of connected diagrams corresponding to 
propagation from x\ to x 2 . Added to this are contributions from disconnected 
diagrams corresponding to x\ multiplied by those to x%. These latter disconnected 
contributions to < (j)(xi)<p(x2) > exactly cancel with the term 4> c i(%i)4>ci{ x 2) ( or 



\j=o =< 4>{xi)4>{x 2 ) > connected i 

(2.7) 



alternatively < (f)(xi) >< (f)(x 2 ) >) in Eq (2.5) resulting in: 

_ 5 2 W 
SJ(xi)SJ(x 2 ) 

which is known as the connected correlator. As it is a propagation amplitude it is 
necessarily greater than or equal to zero and hence: 

5 2 W 

o<0. (2.8) 



8J(xi)SJ(x 2 ) 



One notes from (2.3) that the expression for the classical field 4> c i[J] is now a 
function of J, the source. Here we introduce the Legendre effective action T[(p d ] 
in four dimensions, which is defined by the Legendre transformation. 

= W[J\ - [ rf 4 x0 ci (x)J(x). (2.9) 



We take the first derivative of T[(f> c i] with respect to the classical field at a point 
xi, {<pci{xi)). 

sr[M _ sw[j] f "VL- J{xi) (2.io) 



5<p c i(x 1 ) 5(j)ci(xi) J " 5(j) c i{xx) 
and utilizing the chain rule, 

*TMJ _ f* tMx) 4M--j {xi) . (2.11) 



S4>ci(xi) J 5(f) d (xi) 5J(x) J S(p d (x 
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Using (2.3), that is: 
becomes: 



sw 

SJ(x) 



(> c ;(x), in the second term on the right-hand-side, this 



Taking a second derivative of (2.12) with respect to 
5 2 T 



'hi 



wi 

8J(xi 



-J(x 



(2.12) 



^(^2) gives the result: 



5(j) c i(xi)6(j) c i(x2) S(j) c i(x 2 ) 



XlX 2 



where 



5J(x x ) 



is a matrix in x\ and X2 representing the inverse of 



X\X2 



Using the relationship 



sw 

SJ{x 2 ) 

6 2 Tk 



(X2), we now have: 



( t 2 w[J] V 1 

\5J(x 1 )5J(x2)J XiX2 



(2.13) 



SJ{ Xl ) 



(2.14) 



S(j) c i(x 1 )6(j) c i(x2) 

We now define the effective potential, U e ff in terms of the Legendre effective action, 
as the part of the effective action not containing any parts of the expansion of the 
kinetic component of the bare action. 



r[0. 



cl\ 



d 4 x[U eff + Z(0 d )(«90 ci ) 2 + Oidfa) 



(2.15) 



Here Z((/> c z)(<9</>^) + 0(<90 c ;) 4 represents an expansion of the kinetic terms in the 
bare expression, which vanish with a field configuration constant with x. In such a 
configuration, when minimizing the effective potential with respect to the classical 
fields we can find the true vacuum state of a quantum field theory. In general U e jf 
cannot be computed exactly and approximations are deployed. With source set to 
zero and constant in space time, T[4> c i] 



n 



(FxU, 



eff 



W from (j^. ((2TJ) reduces to: 

v-u eff , 



(2.16) 



where V is the space time volume arising from the integral J d 4 x. One may see the 
use of the effective potential by considering Equs. (2.12) and (2.16) with source 
J = 0. 

0<Pcl{Xl) 

which provides the condition that the effective action T[0 C ;] has an extrema, with 



the solutions to eq (2.17) representing the stable quantum states of the theory. 
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2.1.1 Effective Scalar Theory 



In perturbative renormalisation theory, it is conventional to split the Lagrangian 
into two parts, one containing the renormalised (physical) fields and couplings (£\) 
and the other containing the counter terms (S£) which absorb the divergences 
generated by loop corrections. We thus have the expression 0: 



£ = £i(<P) + S£(cf>), 



(2.18) 



from which the Legendre effective action, valid at one loop only, can be computed 
(quoted without derivation) in Euclidean space, Section 11.4, pQ): 

8 2 £ x 



d x£\{ 



b d ) + ^ lndet 



d x5£{cj) c 



(2.19) 



To first order in 4 theory, [2] computes the leading order expression for T((f) d 
(where Si = f d 4 x£i) and counter terms are not written: 



y = Si(0d) + \tr\n[-d 2 + U"{(f> d )\ + 0{h 2 ). 



(2.20) 



In a constant field configuration the second term on the right hand side of (2.20) 
reduces to (in phase space): 



-trH-d 2 + U'\<p cl )} 



n 



d x 



d A k 

Wf 4 



In 



k 2 + U h 



k 2 



(2.21) 



where we have divided the factor k 2 + U"(4>ci) in the logarithm by a factor k 2 to 



avoid taking the logarithm of a dimensionful quantity. Combining this with (2.15), 
one obtains to first order and where B<fi 2 cl and C<p d are the counter terms quadratic 



and quartic in the fields, contained within 5£(0 c ;)-containing term in (2.19) 

k 2 + U"(<p d ) 



Ueff{.<Pcl) = U((j) d ) + 



d 4 k 



In 



k 2 



(2.22) 



The integral on the right-hand-side of (2.22) contains divergences quadratic in the 



cutoff which can be absorbed by the counter terms (again only first order in H 
terms are considered here). Here B and C are co-efficients containing divergent 
corrections to the powers of (f) d to absorb the cutoff dependence in the integral on 
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the right-hand-side. Evaluating (2.22), the following is obtained, with the integral 
taken up to k = A with h set to 1 for brevity (a convention largely followed in the 
remainder of the thesis). 

A 



Ueff((f>d) = U(4> d ) + 



:U"( 



l + Cfa. (2.23) 



32tt 2 ^ j 64tt 2 U"((j) c i) 

Later in this thesis we will study <p A theory in the context of spontaneous symmetry 
breaking and other phenomena, taking C to be: 

, 1 



c = \{d<i> d ) 2 + l -A 



4! 



(2.24) 



Here A(d(j) c i) 2 , B<p 2 d and C(p d are the counter terms and the bare potential can 
be considered as U(<j> d ) = -\^ 2 {<p d f + ^X(<j) d ) 4 + B(j) 2 cl + C<p, 
case of fi — where we have the condition 



d 2 U"(4>d) 



\4>ci=o 



d . We consider the 
0. This will be of 



interest later in the thesis as it represents the transition point between a purely 



convex potential to a double well shape. We can evaluate (2.23) in this context 
(to quadratic order in A): 



U, 



A 



e// |/i=0 



(~647r 2 



X + B 



(2.25) 



We want to absorb the dependence of the cut-off A in the counter terms B and C 
and impose renormalisation conditions. 



d 2 U, 



effWclj 



d 4 U, 



<:ff(<f>d) 



\<f>d=0 



p c l=m 







A(m). 



(2.26) 



The first condition in (2.26) implies that = (or the renormalised mass-squared 
term vanishes) at (f) d = in the 4 theory and that B = — (A 2 /647r 2 )A). In the 
second condition, since d Ue /J^ cl ^ is not defined at 



b d = due to the log term, we 
choose an arbitrary energy scale m. Here we are left with: 



u, 



e//l/i=0 



4! + (16vr 2 ^ 



which if plugged into the second of (|2.26|) gives: 

A(m) 



m 



X + K 1 X z ln- + K 2 , 
A 



(2.27) 



(2.28) 
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where K\ = and K 2 is a constant not computed here. These two conditions 
can eliminate the A dependence. (We do not derive the value of C required to 
remove the A dependence above). 

tftf/W = ^("Wi + (in § " f ) + 0(Mmn (2.29) 

where we now have the effective potential in terms of the classical field and the 
coupling A(m) for which we have a beta function relationship, which follows from 



(2.28) to order £>(A(m) 



m ^ = TtM™) 2 + 0{\{mf). (2.30) 
am lo7r z 



2.2 Spontaneous Symmetry Breaking 

Considering a particle at position q in one dimension defined by the classical 
Lagrangian, with velocity q and with k' and A two real constants with A > 0, 
we have the Lagrangian (here in Minkowski spacetime). 

= \e-U{q) (2.31) 

This exhibits a discrete reflective symmetry, invariant under the operation q — ^ —q. 
In both classical physics and quantum mechanics the equilibrium and ground states 
respectively are found by minimizing the potential. 

dU 

— = 2k' q + Xq 3 (2.32) 
oq 

If k' > there is one real minimum at <p = and the ground state respects 
the reflective symmetry and the system represents a damped harmonic oscillator. 
However when k' < there is a local maximum at q = 0, (corresponding to a 
potential of U max ) and two minima at q — ±y / ^-, which represents the double 
well or Mexican hat potential in one dimension. At energies less than U max , 
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no 

Figure 2.1: A one dimensional double well potential 

classically the particle must be in one of the two minima, thus breaking the 
reflective symmetry of the system. Since we have not added any terms by hand 
to do this, it is termed spontaneous symmetry breaking. (It should be noted that 
given a large number of particles q a initially with potential energy > U max this 
reflective symmetry will be restored in a sense as there will be equal probability 
that a particle will end up in either minima in this one dimensional model). 

We now promote the particle position variable q to N dynamical fields 4> n (x) with 
x the space time coordinate in D dimensions, n = 1, 2...N and k' replaced by —fx 2 



= \(d<p n ) 2 -U(<p n ). (2.33) 

This now exhibits a continuous O(N) symmetry with the fields transforming 
as an iV-dimensional vectors under the O(N) transformations. The potential is 
minimized at values of <p n such that 




(2.34) 



The length of the vector (f) n _ min is defined but its phase or direction is arbitrary. 
If one chooses the direction to be in the N direction, we have the set of shifted 
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fields such that: 

(pNmin = (0,0,0.. ..0 O ), (2.35) 

where here 4>o = The system can be perturbed around 0o such that 0at(x) = 
0o + cr(x), and define a new set of fields: 

<f>' N {x) = (<Mx), (0o + a(a;))). (2.36) 



We can now rewrite (2.33) in terms of cj)' N (x) and we obtain: 



(2.37) 

This can be interpreted as the perturbation field cr(x) with a mass factor 
a/2/i, with the fields 4> n (x) for n = 1,2...., (N — 1) being massless degrees of 
freedom. The O(N) symmetry has been spontaneously broken to an 0(N — 1) 
symmetry. Goldstone's theorem states that for every spontaneously broken 
continuous symmetry a massless bosonic degree of freedom results. In this case 
the O(N) symmetry (which has N(N — l)/2 generators of the symmetries between 
the N fields) is broken by setting the field in the N direction, resulting in 
(N — 1)(N — 2)/2 symmetries between the fields. The difference in the number 
of symmetries before and after (N — 1) is also the number of resulting Nambu- 
Goldstone bosons, pQ. 

The above description does not incorporate quantization. This will affect the terms 



present in (2.33) and (2.37) and thus will impact the symmetry breaking process. 
Considering the case of iV = 2 and breaking of 0(2) to 0(1), in the ground state 
we have (fix — and <p2 = <Po, with the (pi field representing the massless Goldstone 
boson. Considering quantum fluctuations about the ground state of <pi we can 
consider the mean square of these fluctuations, in momentum space, k. ([2], IV. 1). 

<(<M0)) 2 > = \ /p0(0!(O)) 2 exp(^(0)) 



Z „ 

x — z j 

f d D k expikx 
lim / - — — , , 2.38 
x^o J (2tt) d k 2 ' K J 
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where if 0i represented a massive field p- would be replaced by k2 + 2 ■ D is the 
spacetime dimensionality. There is an infrared divergence for D < 2 leading to 
the Coleman-Mermin- Wagner theorem which states that spontaneous symmetry 
breaking cannot occur for D < 2 dimensions, [I]. This illustrates the effects 
quantization can have on classical spontaneous symmetry breaking. We now 



consider a variation of (2.33) - a complex scalar field theory invariant under a 



U(l) transformation (denoted by — > expii 

C = d^d^cj) + /i 2 0> - A(0» 2 . (2.39) 

We now promote the global symmetry expi8 to expi6(x) (that is, to a locally 
variant gauge symmetry). The invariant Lagrangian is as follows, where is the 
gauge field and _D M = (d^ — ieA^cf) is the operator required to preserve symmetry. 

C = -\f^ v + (D^YDy + /i 2 0*0 - A(0*0) 2 . (2.40) 



Spontaneous symmetry breaking of (2.39) results in one massless boson. However 



a similar analysis, [2], shows that spontaneous symmetry breaking of (2.40) results 
in the gauge field A^ acquiring mass and one degree of freedom which it acquires 
from the massless scalar which disappears. This is the basis of the Higgs mechanism 
which is postulated to be the mass-acquiring mechanism for the particles in the 



standard model. Spontaneous symmetry breaking in a Lagrangian such as (2.33) at 
the classical level may be accompanied by further explicit breaking of the symmetry 
by quantum effects. Additional interaction terms involving the Goldstone boson 
may arise in the Lagrangian and it can acquire mass. This contrasts with the Higgs 
mechanism for gaining mass which is purely a classical effect. 



2.3 Wilsonian renormalisation Group Theory 



The descriptions below are based on pQ, [2], [3] and where noted. 
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2.3.1 The Renormalisation Group Equations 

The functional approach we use to derive exact evolution equations for the QCD 



axion and the bosonic string (described in Section (2.3.3) has its roots in the exact 
renormalisation group equations which in turn rely on the renormalisation group 
techniques pioneered by Wilson in the early 1970s. In this thesis we thus briefly 
introduce the latter and then the former prior to describing our technique. 



Initial misgivings about divergences in quantum field theories were dispelled by 
the application of regularization and renormalisation techniques in the second 
half of the 20th century to extract meaningful information from both infrared 
and ultraviolet divergences. In regularization, an infinitesimal, usefully definable, 
spatial distance e (or A ~ - in momentum space) is set and calculations performed 
in the limit e — > 0. In physical terms, regularization equates to not allowing 
the internal lines in Feynman diagrams to fluctuate at momentum k > A. If 
the resulting expression is definable it may contain terms proportional to A. 
Once regularization has been achieved, the renormalisation process refers to 
techniques to determine whether these A-proportional terms are significant (a 
non-renormalisable theory) or may be canceled (a renormalisable theory). A non- 
renormalisable theory such as gravity in four dimensions, for example, may not be 
tackled using perturbative quantum field theory techniques. For a renormalisable 
theory, a limited number of measurable field parameters can be calculated which 
are independent of A, which is an arbitrary energy scale above which even the 
bare theory may or may not be applicable. Until the late 1960s, it was not well 
understood why high energy self interactions in renormalisable theories had so little 
impact at lower, infrared energy scales where measurements could be taken, and 
thus why the techniques worked so well, renormalisation group theory emerged 
to address this and it describes how the physical couplings (e.g. \p in 4 theory, 
where P denotes physical) vary with momentum scale k where k < A. It is useful 
in investigating the high and low energy behaviour of quantum field theories. To 
look at this concept further, we consider the generating functional of a scalar field 
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theory in Euclidean space time with field parameter (source term here is zero 
for brevity): 

Z[J] A = J [D<j>} A eM-S[<f>}), (2-41) 

with [£>0]a = ri|fc|<A c ^ ) (^)- Wilsonian renormalisation takes shells of <p(k) at 
\k\ < A, integrates out these field degrees of freedom. The "group" transformations 
are energy scale transformations, k — > k + 5k, and the impact of these on the 
effective parameters of the theory, such as Ap, are studied. In a scale invariant 
theory the physics remains the same under such transformations. However, in all 
useful, fully quantized field theories of the standard model, scale invariance appears 
to be broken, pp. 

In implementing renormalisation group techniques, the starting point is a cutoff 



at scale k = A which represents the bare or classical theory, as in (2.41). To arrive 
at an expression for an effective action at some energy scale k < A the high energy 
field degrees of freedom 0(A) — > (p{k) are integrated out. To do this we divide the 
field degrees of freedom into two sharply defined sectors: 

<p(k) =$(k) for \k\ < A - 5 A (0 otherwise) (2.42) 
(j){k) = ip(k) for A - 5A < \k\ < A (0 otherwise), (2.43) 

where <p(k) = <3>(fc) + (f(k). The generating functional up to momentum scale A is: 
Z A = [ P(0)exp(-£ A [0])= / £>($) / 2%)exp(-£[$ + <d). 

JA Jo<k<A-8A JA-6A<k<A 

(2.44) 

This can be written as: 

Z A = [ D($)exp(-S A _ 5A ($)), (2.45) 

J0<k<A-8A 

where the Wilsonian effective action S A -sa{^) is defined as involving only the 
components of <f)(k) with \k\ < A — 5 A: 

exp(-5 A _ 5A ($)) = / 2%) exp(-5[$ + cp)). (2.46) 

JA-SA<\k\<A 
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In 4 scalar theory, the effective action S'a-5a($) includes corrections proportional 
to powers of the bare coupling A which contain information on the quantum 
interactions of the large k components in the sector <p(k). Wilsonian group theory 
leads to expressions for the variation of the effective parameters of a theory with 
energy scale k, with, as before k < A. The /3 functions are differential flow 
equations, where Aj are the effective couplings: 

kd k \i(k) = /3(Ai). (2.47) 

Fixed points of a theory are defined where kdkSk = 0, that is, the effective action 
is unchanged by a slight transformation, 6k in energy scale. The simplest case 
of a fixed point is the free-field scalar theory where the action is given by = 
J d D x^(d^(f)) 2 (also known as a trivial or Gaussian fixed point). Here the couplings 
vanish and in the vicinity, the form of the (5 functions can be studied. Couplings 
that grow in strength approaching this fixed point are termed relevant while those 
that die out are termed irrelevant. Quantum field theories with irrelevant couplings 
are non-renormalisable. 



2.3.2 The Exact Renormalisation Group Equation 



The equation (2.46), depending on the nature of S the effective action, is generally 
solved by perturbative methods, where one or two loop approximations may suffice. 
In this thesis, we are interested in very high energy theories where perturbation 
theory may not be well established and where exact expressions are useful. The 
exact renormalisation group (ERG) techniques, [8] build on Wilsonian methods. 
The resulting equations are non-perturbative in that they capture all quantum 
fluctuations in an effective expression and are based on applying a cutoff as in 



(2.42) and integrating out higher energy degrees of freedom. Expanding Sk(Q + <p) 
around $ gives, (where the phase space volume within p — > A is V = J £? b ) : 



(2.48 
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From (2.48) one can obtain an expression for the evolution of the Wilsonian 
effective action in a constant field $o configuration: 

0Sfc($o) - <Vi*($o)) 



d k S k (%) 



lim 



5k 



(2.49) 



Evaluation of the Gaussian integrals in (2.48) leads to the Wegner- Houghton 
equation, [6], where here \p\ = k and Qd is the solid angle in D dimensions. 



d k S k (<5> ) 



-ttnk D - 1 



V, 



D 



(27T) 



D 



In 



5 2 S k (<!> ) 
5$(p)5$(—p) 



(2.50) 



A full derivation of this is provided in Appendix |8.1| where we note that the term 

and the 



containing S $£^ vanishes as we consider a saddle point where g ^, ) 
field $ is constant with spacetime. A key point to note concerning the derivation 



of eq (2.50) is that we consider an infinitesimal momentum shell 5k « k and 
the equation becomes an exact expression as we take the limit y — > 0. We make 
two points about the introduction of the infrared cutoff, A. Firstly, if the cutoff 
is sharp, the Wilsonian exact renormalisation group methods may only be used 
to study the behaviour of the effective potential, as the derivative, kinetic terms 
will prove problematic at the sharp boundary. In this instance, a smooth cutoff 
will be required, [7]. Secondly, the introduction of an energy scale cutoff is not 
in general gauge invariant, as gauge invariance implies energy scale invariance. In 
our treatment, we note that we use a running variable other than energy scale or 
momentum, which may overcome gauge invariance issues associated with evolution 
with energy scale as in the original Wilsonian treatments, [6]. 



2.3.3 An Alternative Exact Approach 

Following this brief description of the exact renormalisation group, we describe an 
alternative non-perturbative technique for deriving an exact form for an evolution 
equation for the effective action. Motivated by the exact renormalisation group 
equations they were developed in [8] and utilised in |103j . |104j . |105j . |106] and 



[107] . These results, along with the results derived in Section 2.1 on the effective 



potential are key building blocks in the non-perturbative techniques used in this 
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thesis. We are interested in the evolution of quantized parameters like r(0 ci ). |1(J2] 
describes a functional method to arrive at an evolution equation of the effective 
action with a parameter of the bare theory and the equivalence of such an approach 
with the evolution using the exact renormalisation group. As an example, we 
consider a bare scalar theory, in Euclidean space <fi = <f>(x): 



S(<f>) 



dPx-dp 



+ -bm 



2 J.2 



A 
4! 



/.2\2 



) 



(2.51) 



to which we have added a unit-less parameter b to illustrate our non-perturbative 
evolution technique. When b = oo the theory is classical because the particle 
becomes infinitely massive. The decrease of b coincides with the increased 
importance of quantum fluctuations, the full extent of which are apparent when 
6—7-1. We note the analogy with the exact renormalisation group derivations 
where the cut off k is taken from A — > to capture quantum effects. Starting with 
(2.9), we take the partial derivative with respect to b and we consider b and 4> c i 
to be independent variables of T[0 C /] and consider a constant field configuration so 
0. 



that ^ 



d h T\ 



d b W[J] + J d 4 x 



dW[J] dJjx) 
dJ(x) db 



d 4 x(j) c 



dJ(x) 
db 



d b W[J}. (2.52) 



where we have used the definition of the classical field in (2.3) and performed a 



chain rule manipulation similar to that used in (2.11 ). We have from our definitions 



(2.1) and (2.2): 



d b W[J(x)] = -d b Z 



1 

z 



-S((f)) - / d n xJ(x)cj)(x) 



d b S(4>), (2.53) 



and can evaluate c^W[J] directly, using (2.3) with (2.53). 



d b W[J] 



m 

T 

2 

m 

T 



d 4 x < <p(x)(j)(x) > 

d A x I d i y5 A {x — y) < 4>{x)(j){y) > . 



(2.54) 
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Using the results (|2.3|), (|2.5|), (|2.14|) and (|2.52) the following is obtained. 



an 



m 

T 

9 

m 
T 



<i xd 4 yS 4 (x — y) 
rf 4 x(0dW) 2 + Tr 



8(j) d {x)5(j) d {y) 
S 2 T\6J 



+ (f>d(x)<f>d(y) 



-i 



(2.55) 



5(j) c i(x)5(j) c i(y) t 

where we define Tr[...] = J d 4 xd A y5 4 (x — y) [...]. We note that as with the exact 
renormalisation result obtained in (8.21) the equation (2.55) is an exact, non- 



perturbative expression for the effective action. While (2.50) uses the momentum 



cutoff k as an evolution parameter, the result (2.55) is not restricted to this form 
and uses any suitable parameter, the evolution of which captures all quantum 



fluctuations. Typically equations such as (2.55) are not directly solvable and an 
assumed form for the effective action is used in conjunction with this to investigate 
viable solutions. For example, the gradient approximation for scalar theory was 
described in the context of the exact renormalisation group in the Appendix in 
(8.10), which may be applied here, (in this <fi = <p c i). 

' ~ ' ■ ' " (2.56) 



Sefftt) = d u x\ -z^d^oy + o(d,<py + U eff ((f>) = r, 



where U e ff(<p) and Z(<p) contain quantum fluctuations and T is the Legendre 
effective action. Here, df,S((j)) e ff may then be compared with (2.55) and viable 
solutions found. This is the approach we take with effective theories for the QCD 
axion and bosonic string in this thesis. 



2.4 Convex Potential and Spinodal Instability 

The convex shape of the effective potential in a scalar field is a well known effect in 
classical physics and we begin this section with a brief discussion of the Maxwell 



construction, [9]. Figure (2.2) represents an isotherm predicted by the Van der 
Waals equation of state (the oscillating curve), where P is pressure and V is the 
volume occupied by the gas at that pressure. The middle section of the curve 
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where |^ > is not physical in that the pressure is decreasing with increasing 
volume. The two sections of the curve with|^ < bounding the green shading are 
known as meta stable states. In practice the two green sections are replaced by 




Figure 2.2: The curve is an isotherm in a real-gas model (such as that associated with a 
van der Waal's equation of state). The unstable concave region and the convex minimum 
are replaced by the straight line at constant pressure P. |9j. 

the straight line of constant pressure P = P eq where a phase transition from liquid 
to gas is occurring and both phases are present in the system. This two-phase 
system remains at constant pressure with increasing volume until all the liquid 



has vapourised (where the isotherm in Figure (2.2) intersects the P = P eq line on 
the right hand side). In the Maxwell construction which is used to resolve the 
unphysical nature of the equation of state isotherms, it is shown that the internal 
energy of the system (assumed to have a fixed numbers of particles) depends only 
on the volume for a given isotherm. It is shown that the change in energy between 
the two outer volume points where the P = P eq line intersects the isotherm may 
be calculated along either the P = P eq line of constant pressure or along the Van 
der Waals isotherm - both give identical results. The P = P eq line represents the 
physical reality and thus the concave point of the isotherm is flattened with the 
straight line. The method, which is inserted by hand rather than analytically, 
relies on the fact that the Gibbs free energy of the two phases must be equal when 
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they coexist and the two areas shaded in green are of equal area. It may be used 
in any thermodynamical system with P and V replaced with any pair of conjugate 
variables which can be used to express the internal energy of the system. We now 
consider this in the context of quantum field theory. 



\ 










i 1 

>2 03 



Figure 2.3: The curve indicates a possible form for the effective potential of a scalar 
field, with extrema at the three points 01,2,3- In our notation U e ff = V pQ. 




Figure 2.4: In practice the concave part of the effective potential is flattened. In our 
notation U e ff = V p]. 



PP notes that there is a field theory analogy for the Maxwell construction in 
thermodynamics in which the conjugate variables are the effective potential U e ff 
and classical field C ;. Figure (2.3) represents a form of £7 e //(0d) with an absolute 
minimum at 1; a local minimum at 03 and local maximum at 4>2, and a constant 



background field configuration. In eq. (2.17) it was noted that the extrema of 
U e ff {(pel) represent the stable quantum states of the theory, in this case 0i, 03 and 
02. In the case of our potential form in figure (2.3), it is clear that the stable 
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states are the minima: 4>i and 03; with 3 being a locally stable state which could 
decay to <pi via tunneling. 2 does not represent a stable state even though it is 
an extrema of U e ff(<f) c i). In a manner analogous to the Maxwell construction, if we 
consider a value of 4> c i between <pi and 3 , it can be described by a superposition 
of the two stable vacuum states (where U e ff is a minimum) such that, where a is 
between and 1. 

(pd = afa + (1 - a)4> 3 . (2.57) 
We may also describe the average value of U e ff at <p c i as: 

U eff (<f> d ) = aUefffa) + (1 - a)U eff (<j> 3 ), (2.58) 



whose line is represented by the bold line in Fig (2.4). Thus our requirement to see 
that the true vacuum state is represented by the minimum rather than maximum 
of U e ff leads to a flattening of the concave part of U e ff. In effect this implies the 
condition: 



5 2 U, 



eff{<Pcl) 



> 0. (2.59) 



It should be stressed that this scenario applies to a configuration with source J set 



to zero and in a constant field configuration as described by eq. (2.17). As with 



the theoretical form of the isotherm in Figure (2.2) yielding unphysical regions 



where phase transitions mean that the Maxwell construction is used to remove 



these by hand, unphysical regions in the effective potential's form as in figure (2.3) 
are removed by hand. The value of the local minima are unaffected by this. 



The shape of figure (2.3) is similar to the double well scalar potential which 



gives rise to spontaneous symmetry breaking as described in section |2.2| From 
the very qualitative account above it should be that when quantized, a concave 
scalar potential (with source set to zero and constant field configuration) should be 
flattened when the system is quantized. This is indeed a well known phenomenon 
and one that we utilize in our discussion of the flattening of the axion's effective 



potential prior to any interactions in section 4.4 In practice this flattening is 
due to spinodal instability effects which are briefly outlined, with this qualitative 
account taken from [ID], [II], [12], p3] and [33]. 
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We consider the potential: 



m = ~<? + ^ + £- (2-60) 



In (2.60), U(4>) has a local maximum at <fi = (when U" = — m 2 ) and two minima 
at = ±t> (f 2 = ^p-). The equation of motion for the scalar field fluctuations 
about the <fi = is quoted without derivation, where <pk are the momentum (k) 
modes of the field and U" is with respect to the fields and where the momentum 
squared in Euclidean coordinates can be expressed in terms of the time and spatial 
components, k 2 = co 2 + k 2 . 

{k 2 + U")4> k = 0. (2.61) 

While symmetry is still evident at <fi = and time t = 0, the field is massless for 
the quantum fluctuations <pk = -^exp(—iujt + ikx) (that is with m and terms 
inactive). As <fi deviates slightly from 0, the m term is "turned on", from when 
\k\ < m, the k term in (2.61 ) results in exponential growth of these modes implying 
a lack of restorative force driving the small quantum fluctuations. Qualitatively, 

— * 4 

[T2] notes that these \k\ < m wave modes grow in amplitude and reach energy ~ 
which is of similar size to the initial maximum of the bare potential U(0) = 
Thus a large part of U(0) is transferred to the kinetic energy of the field as it rolls 
to the minimum at U (</)). The flattening of the concave classical potential in the 
quantum model occurs to avoid exponential build up of negative mass squared- 



originating (hence tachyonic) undamped fluctuations (Figure 2.5). It is again 
stressed that at this point that this analysis is limited to non-interacting scalar field 



models. The graph in figure (2.5) was inspired by results computed numerically in 
[TU] . Here the authors combine the non-perturbative methods outlined in section 
(2.3.2) with the saddle-point (or steepest descent) approximation method for 
evaluation of divergent integrals. They found that when the spinodal instabilities 
occur near k 2 = m 2 , non-trivial saddle points appear in the integrals of the type 
expressed in (2.48). At tree level this gives k 2 — m 2 + |0 2 = 0, and a flat potential 
for any value of k down to the IR limit k — 0, within the two minima of Uk(4>)- 
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Figure 2.5: The flattening of a double well bare potential in one dimension by quantum 
effects, [Uj. In the context of exact infinitesimal Wilsonian renormalizatin group studies, 
the spinodal instablity is compensated, at each step, by the presence of a non-trivial 
saddle point in the path integral defining the blocking transformation, [TO]. The result 
of this cancellation is the flattening of the Wilsonian effective potential (the "average 
effective potential"), as the observational scale k goes to 0, as shown in the figure. It 
is noted that non-trivial saddle points in the path integrals defining the renormalisation 
group transformations imply a "tree-level renormalisation", since the bare non-convex 
part of the initial potential is removed by quantum fluctuations. 



Chapter 3 
Axion Theory 



3.1 Axion Related Topics 



Axion physics involves relationships between a number of concepts and phenomena 
in quantum field theory and particle physics which are reviewed here briefly both 
to introduce these relationships and terminology used. Note: in this and the next 
Chapters the term "axion" refers to the scalar field postulated by Peccei-Quinn to 
explain the CP problem, as opposed to the term "string-axion" used in Chapters 5 



and 6. Section (3.2.1.3) of this chapter briefly mentions string-originating axions, 
however. 



3.1.1 QCD and its Symmetries 



Quantum chromodynamics (QCD) is a highly symmetric quantum field theory and 
much of QCD research has centered around how these symmetries are broken or 
preserved and the resulting phenomenological implications. Initially the strong 
interaction was believed to be explained by a new quantum number, isospin, 
described by SU(2)j symmetry in which the proton and neutron existed in the 
fundamental representation and the thee pions, postulated to be Goldstone bosons 
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of the theory, in the adjoint representation of the SU(2)j group. Because the 
pions were shown to have a light mass, it was postulated that this symmetry was 
slightly broken explicitly. Modern interpretations refer to a more fundamental, 
quark-based model with associated symmetries of the Lagrangian and vacuum, 
the key ones of which are described. 

3.1.1.1 SU(3) Colour Symmetry 

This is the defining symmetry of QCD and results in a non-Abelian gauge theory: 

Cqcd = qSrd, - m)q t - g(q a ^ jqj )G; - ±G*»G%,. (3.1) 

Here qi{x) and qi(x) are quark and anti-quark fields of i flavours in the fundamental 
representation of 577(3). G° are the eight gluon fields which lie in the adjoint 
representation of SU (3) with the generators for this representation and 7 M are 
the Dirac matrices and m is the quark mass matrix. G a is the gauge covariant 
field strength tensor given by: 

= d^G?, — d p G a ^ — gf ahc Gb^G cv , (3.2) 

where g is the coupling constant of the theory and f abc are the structure constants. 
Throughout we use a trivial summation notation for the a,b,c and i,j indexes, 
with the '"up-down"' convention respected for the space time indexes. This gauge 
symmetry is an exact one which remains unbroken and along with electroweak 
SU(2) x £7(1) is one of the fundamental symmetries of unified theory. 

3.1.1.2 Baryon Number 

Baryon number is defined as B = \{n q — rig) where n q and rig are the number of 
quarks and anti quarks in a baryon. It is a conserved quantum number and hence 
exhibits a U(1)b symmetry, analogous to the U(l) electromagnetic symmetry. If 
baryon number is conserved precisely, it is an exact symmetry. We consider a 
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generic Lagrangian with a quark and anti quark doublet: 



u 




Q=[ d );QL=[ ;) ;Qr=\ , )■ (3.3) 




Here we have the relationships: 



Ql = [ Q; Qr = (±±P- ) Q, (3.4) 



2 

where 7 s is the chirality matrix: 

There is a conserved vector current j M = Qj^Q from the U(l) transformation Q — > 
(expi6)Q. If m u = = the left and right hand fields transform independently 
in the chiral symmetry that transforms Q — > exp(i(3 / y 5 )Q where /3 is a phase 
degree of freedom introduced by the symmetry. If we consider m u ~ 7^ and 
m u ,m,d Aqc£>, the quarks have a small mass meaning the chiral (also termed 
"axial vector") symmetry is explicitly broken as the quarks acquire mass. The 
chiral current in this model is characterised by: 

= 2tmQ 1 5 Q 1 (3.6) 

and the chiral current is conserved in the massless case. However, (see for example 
PQ, 19.2) when the above classical symmetries are quantized the Adler-Bell-Jackiw 
anomaly in four dimensions results in breaking of the chiral U(1)b symmetry even 
in the massless case of (3.1) such that: 

d^ = -^£~ 2 G^G^ (3.7) 

where iV is the number of quarks in the model, and is the dual of and 
equal to |e Q/3 ^G a/ 3. Detailed analysis [TJ] by 't Hooft and others in the early 
1970s revealed the chiral anomaly in a four space time dimensional QCD vacuum 
to be a one-loop effect, breaking what would have been a fundamental U(1)b 
chiral symmetry. The experimental absence of a pseudo-scalar of mass smaller 
than v^wii, plj, indicative of a spontaneously broken U(1)b chiral symmetry is 
known as the "U(l) problem". 
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3.1.1.3 Flavour Symmetry 

The idea of SU(2)j isospin symmetry gave way to flavour symmetry SU (Nf) where 



Nf is the number of quark flavours (denoted as % in (3.1)). At present it is thought 
there are six flavours of quarks. If all the flavours are massless there is also a chiral 
SU(Nf) symmetry. For our analysis of axion physics we consider only the lightest 
two quarks, u and d, such that m u ~ m d <C Aqcd, where Aq C d is the energy scale 
at which confinement occurs, ~ 200MeV. As with the U(1)b symmetry there is a 
SU (2) vector symmetry and if m u = ma = a chiral SU (2) symmetry. The vector 
symmetry is one of rotational invariance between m u = rrid and if the masses vary 
slightly it is an approximate symmetry. The vector and chiral currents are: 

f a = QYr a Q 

with t 0, the Pauli matrices. One may calculate the quantum anomaly contribution 
to the divergence of the SU(2) chiral current as, pQ: 

2 2 

= -^T-fi^G^ ■ tr[r a t c t d ] = -J^G^G, u tr[r a }tr[t c t d ] = 0, (3.9) 

as the trace over a single r a vanishes. It is noted that r a is an isospin matrix while 
t c is a colour matrix. Thus unlike in the case of the U(1)b chiral symmetry, in 
SU(2), the chiral anomaly does not break the symmetry. Due to the presence of 
a triplet of low mass mesons (the pions), it is postulated however that the chiral 
SU (2) symmetry is spontaneously broken. This is due to the fact that even when 
the mass of the quarks is set to zero, the QCD vacuum is thought to consist of a 
condensate of quark-anti quark pairs arising from the vacuum such that: 

< o|gg|o >=< o\q l q r + q r q l \o >^ o. (3.10) 

This non-zero expectation value of the quark-anti quark pairs in the condensate 
causes breaking of the SU (2) chiral symmetry, if we consider a model with two 
quarks, up and down. The quark-anti quark pairs act as scalar bound states and 
three massless quark-anti quark Goldstone bosons result - the pions. Since in 
reality the up and down quarks have a small mass the chiral 577(2) symmetry is 
also explicitly broken and these pions acquire mass. 
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3.1.1.4 Scale Invariance 

Scale invariance is a sub-symmetry of the conformal symmetry group. It occurs 
in a massless field theory, with field parameter ^{x^) and with no dimensionful 
coupling constants that is invariant under the following transformation, pQ. 

ip(x^) -> exp(-L>cr)<p(e~ ,7 £ M ), (3.11) 

where D is the canonical mass-dimension of the field (p and the space time 
coordinate transforms (exp — a)x^. (We write the constant in the 

form exp —a as a will be promoted to u{x^) in the context of later string- 
related discussions of conformal transformations). These are equivalent to the 



renormalisation group transformations discussed in Section 2.3 Classically this 



symmetry is evident in massless QCD, as in (3.1) with m = and where the 
coupling g is dimensionless. Classically, the physics of such theories remain 
invariant if the length scales are multiplied by a common factor. In quantum field 
versions of such a theory one can equate such invariance with length scales inverse 
to mass scales. Following renormalisation, the effective Lagrangian is dependent 
not only on the classical coupling parameters but also an arbitrary mass scale 
M introduced through the regularization of ultraviolet divergences. The scale 
invariance has potentially been broken by the chosen cut off. This is evident in 4 
theory by the flow in the coupling A with the arbitrary mass parameter M, pQ. 

x ^ x + Wf (N + s) 'i^' (3 ' 12) 

where O(N) is the classical symmetry of the theory. In general this leads to a 
formulation for the beta function of the coupling which describes the variation 
with energy scale /i. 

0(A) = fijp (3-13) 

If the beta function vanishes, the quantum field theory is scale invariant - no 
physically useful examples of these are known in the standard model, pQ. If /3(A) > 
0, as in QED, the coupling approaches zero at low energies to allow for perturbative 
solutions in this region. If /3(A) < 0, the theory is said to be asymptotically free, 
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a phenomenon discovered in the early 1970s as applying to non-Abelian gauge 
theories like massless QCD where the form of the beta function was found to be, 

M ~ _A) To^' (3 " 14) 

with /So a constant and g the strong coupling. Thus the coupling g is high at low 
energies making effective perturbation theory difficult in this region. Conversely, 
high energy behaviour of such theories is predictable via perturbation methods. An 
outcome predicted by asymptotic freedom is that the energy required to separate 
particles governed by such interactions could tend to infinity, potentially explaining 
why quarks are not seen alone in nature. 



3.1.1.5 CPT Symmetry 

CPT symmetry is now considered to be a fundamental symmetry of the standard 
model, including QCD. Parity is a discrete transformation which sends (x, t) to 
(— x, £), that is a reflection in space. A discrete time transformation on the other 
hand sends (x, t) to (x, — t). Charge conjugation is a discrete transformation 
which changes a particle into its antiparticle. The CPT Theorem, [T7], states 
that any quantum theory in flat space time is symmetric under the combined 
CPT transformations provided the theory respects Lorentz invariance, locality 
and conservation of probability (unitarity). Theories may break C, P, T or double 
combinations thereof individually. The couplings of the SU(2) gauge bosons in 
the QCD Lagrangian violates C and P symmetries individually 20.3) but the 
combined CP transformation is a symmetry of the QCD Lagrangian as it stands in 



(3.1 ), as is T, the time reversal transformation. There are, however, possible terms 



discussed in Section 3.1.2 which can break CP symmetry in the QCD Lagrangian, 
and the fact that this has not been observed in practice is known as the CP 
problem. 
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3.1.2 Inst anions 



In the context of a non-linear field theory, solitons are stable, well-behaved solutions 
to the classical theory, which are stable against decay (or topologically distinct 
from) to the trivial solution. Most solitons are exact and non-perturbative. The 
stability of these solutions arises, |3j, as a result of constraints imposed by boundary 
conditions of the coordinate space being considered. This boundary has a non- 
trivial homotopy group associated with it, which has a mapping with the coordinate 
space. This can result in an infinite number of topologically distinct solitons which 
can be degenerate. Solitons are often termed topological defects in this context 
and examples are kinks [19J, domain walls and 't Hooft-Polyakov monopoles, [20]. 
Mathematical discovery of the latter trigged searches for magnetic monopoles. 
They are non-linear solutions to gauge theories (such as U(l) electromagnetic 
gauge theory incorporating the Higgs field) which represent objects with finite 
energy, localized around a particular point in space, and can be shown to possess 
magnetic charge. Here, however, we focus on a form of soliton localized in space 
and time, hence instanton. The following is sourced from [3] and where noted. 



In the double well potential illustrated by the Lagrangian (2.31 ) where the constant 
k' obeys the condition k' < 0, one may apply the time independent Schrodinger 
equation simply as (with factors of Ti left out for brevity): 

= (2m(V(q) - J3))V>, (3.15) 

where ip is the particle's wave function, q its position and E its total energy and 
V(q) the potential. A solution is ip = exp(-ik'q) where k! = \/lvn{E — V). If 
k' is real one obtains a familiar plane wave solution, however if k' is imaginary, 
that is E < V(q), as is the case with the particle being in one of the two 
double wells, one obtains an exponential (i.e. exp(— k'q)) solution which represents 
quantum mechanical tunneling between the two wells with amplitude proportional 
to exp —(J^ dx \J 2m{y (q) — E). If we apply the path integral approach to this 
two dimensional example using a Wick rotation in which it — > r, one obtains as 
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the amplitude, with H the Hamiltonian: 

< — q\ exp(— Hr)\q > = J dq(r) exp(— Se) 



St; l T9) dr I ^ 1 \ .. . (3.16) 



r(-q) 



2 



where one notes V(q) changes sign under the Wick rotation and the double well 
is inverted. This transforms the process from — q to q as the particle traveling 
from the top of one maximum and rolling to the other. It is noted that this 



quantum field theory approach differs from that shown in (2.38) in that it does 
not involve perturbation theory. Rather this so-called instanton solution is based 
on applying the path integral approach in Euclidean form to a classical solution in 
a fixed time span and localised area. Instantons add alternative solutions to the 
vacuum structure of a theory, in this case suggesting that the particle may reside 
with equal probability at both — q and q, two areas in the theory disconnected 
classically and by conventional perturbative quantum field theory, but connected 
by quantum mechanical tunneling. [18] notes that the result in Euclidean space- 
time can be equated to Minkowski space time in the path integral approach to a 
good approximation. Instanton solutions to four dimensional non Abelian gauge 
(Yang Mills) theory were discovered [12] in the early 1970s shedding new light on 
the QCD vacuum. While a rigorous derivation is beyond the scope of this thesis, 
the major steps are outlined, drawing from [3], [15] and [18J. 



• Considering the Euclidean action in 4D of the kinetic term of the pure Yang- 
Mills theory, 0: 

Se = ~^J d 4 xTiG^ lu G IMU , (3.17) 

with G^ v = G a T a , and T a are the generators of the symmetry group, G® 
are the gauge fields, G^ is the gauge potential equal to G a ^T a and G^ v = 
\& vpa G pa . We take the symmetry group as the SU(Nf) flavour symmetry of 
a two quark model, where Nf = 2. This is a local, continuous and therefore a 
gauge symmetry. Unlike Abelian U(l) theory, this contains cubic and quartic 
terms representing self interactions of the gauge bosons, G". For Se to be 
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finite in the integral as \x\ — > oo, the potential must be pure gauge in such a 
configuration: 

Gf, = U~%U, (3.18) 

where U is an element of the symmetry group, (e.g. in SU(2), U = exp i9 a T a 
with T a the generators in the adjoint representation) and it is noted G^ v = 
with this gauge at the \x\ — > oo boundary. 

Euclidean space time in 4D has as its boundary the 3-sphere, or S 3 . 
Meanwhile the group SU(2) may be represented by: 

3 

[/ = Uo + iJ^UjVji ( 3 - 19 ) 

3=1 

where Oj are the Pauli matrices. U is unitary with Uq + U\ + C/f + C/f = 1, 
which is also the equation for the 3-sphere, S 3 . Thus one can say that 
the gauge potential at oo describes a map from group space to physical 
space S 3 — > S 3 , with the mapping defined by an integer q (known as the 
Pontryagin index). A solution with one value of q is termed stable if it cannot 
be continuously transformed into a different q solution, and the mapping is 
non-trivial. 



We can define a total divergence, d^K^. 

' =<V^. (3.20) 



Applying the classical equations of motion for (3.17), that is, D^G^ = and 
applying Gauss' theorem gives (where K± is the component of A M normal to 
the surface of the volume under consideration): 

J d 4 x Tr G^G^ = aJ tfx^Kj =±j d ^ K ^ ( 3 - 2 l) 

which implies that the integral J d 4 x Tr G^G^ depends only on the 
homotopy of the mapping S 3 — > S 3 . With our pure gauge condition on 
the \x\ 2 — > oo boundary (i.e. G M = £/ _1 <9 M £7), it can be shown that: 

\ [ d A xTiG, u G, lu = (3.22) 
9 J 9 
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and we define q here as equal to ^2 J d 4 xTr G^G^ as a measure of the 
degree of mapping of the 577(2) group space to physical space S3 — > S3, [3], 
which is the Pontryagin index, q noted above (sometimes referred to as the 



winding number), [18J. The solution of the equations of motion of (3.17) 
with the pure gauge condition at the \x\ 2 — > 00 boundary is an instanton. It 
represents the transition as Euclidean time evolves from negative to positive 
infinity, from one vacuum (represented by a homotopy class n — 1 to another 
in homotopy class n, with Pontryagin index, q in this case equal to n — 
(n — 1) = 1. The non-trivial mapping S3 — > S3 is represented by the set of 
integers which physically means there are an infinite number of identical but 
topologically separate vacuum configurations. The q — 1 case is known as 
the BPST instanton, [22]. The barrier penetration is given by exp(— Se), or 
exp(— ^|-) in this case. 

Thus the QCD vacuum is infinitely degenerate with non-zero transition 
amplitudes between the vaccua belonging to different homotopy classes. If 
I n > is a vacuum described by homotopy class n the real vacuum should be 
invariant under a transformation (termed a '"large gauge'" transformation) 
which maps the n vaccua onto one another. A gauge invariant vacuum state, 
parameterised by 9 is thus constructed as a superposition of the n homotopy 
class vaccua: 

\9>= y ' expin9\n >, (3.23) 



00 

£ 

n=— 00 



with 9 a phase with period 2n. 9 parametrizes the degree of tunneling 
between the n vacuum configurations which occurs in the true vacuum state. 
If 9 7^ 0, 2mr, instanton effects are present and the vacuum state is complex 
and is not invariant under CP transformations (discussed in more detail in 



Section 3.1.3). This instanton parameter 9 can be accounted for in the QCD 
Lagrangian by adding a 9 containing term: 

£e = -9-^TrG^ u G^. (3.24) 

9 parametrizes the degree of tunneling between the n vacuum configurations 
which occurs in the true vacuum state. 
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Considering the Lagrangian (3.1) with m = and % = 2, that is a massless 



two quark model, [15] computes the chiral anomaly associated with the U(1)a 



symmetry breaking as (3. 7) and notes this is a one-loop effect. Comparing 



this with the result in (3.22) gives: 



(3.25) 



where N is the number of quarks in the model and q is the Pontryagin index. 
This implies that in an instanton background there is non-conservation of the 
charge associated with axial current. This implies the possibility of decays 
which violate both baryon and lepton number (the U(l) axial charge relevant 
here), such as: 

p + n -+ e + + zv (3.26) 



The probability of such decays is small, however on the order of exp 



10 



-262 



0. 



While the computations leading to (3.22) involved the gauge fields only, they 



were specific to the gauge group SU(2) relating to a two quark model. In a 
more general SU(Nf) model, the corresponding SU{2) subgroup of this will 
produce the same result, [18] . If both quarks are massless there is also a 
chiral U(l) symmetry relating to conservation of baryon number. The U(l) 



problem mentioned in Section 3.1.1 may be resolved by considering that this 
U(l) chiral symmetry is dynamically broken by instanton effects resulting 
in the chiral anomaly. This is an inherent feature of the quantized theory 
and chiral U(l) is not a true symmetry of the theory and hence the lack of 
pseudo-Goldstone bosons whose mass vanishes in the limit m u = — > 0. 



3.1.3 The CP Problem 

As noted in [T7] , under the CPT Theorem, CPT symmetry may only be violated in 
the case where there is violation of either Lorentz symmetry, locality or unitarity, 
none of which we wish to consider here. We thus start by the assumption that 
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CPT symmetry is preserved in strong interactions. In the 1960s, CP violation was 
observed in weak interactions which implied T violation in order to preserve CPT 
symmetry. In strong interactions, instanton effects result in an additional phase 



degree of freedom to the QCD Lagrangian manifest in (3.24). This term potentially 



violates parity (as immediately seen by the four space time indices in the full form 
^e tJlUpa Gp fT G f11 ') while preserving charge symmetry, [28] and hence can violate CP 
symmetry. Why CP and separately T symmetry violation has not been observed in 
strong interactions, despite being theoretically possible under the standard model, 
is the CP problem. In a simplistic view, one may rotate away the CP-violating 9 



term via a chiral transformation in the massless m = QCD theory (3.1 ) and thus 
preserve the CP and CPT symmetry, [26J. 

N 

9^9-^a u (3.27) 

i=l 

where i is the flavour index. We consider a two quark model with m u = ^ 
where [25] notes that the quarks acquire mass from the Higgs mechanism and 
the quark mass matrix M y - need not necessarily be real nor diagonal. Chiral 



transformations similar to (3.27) can be performed to diagonalize the mass matrix 



to derive meaningful mass parameters, 

q iR -»■ exp(iai—q iR ) (3.28) 
q iL -»■ exp(-iai—q iL ). (3.29) 

However, the price paid for this is that CP breaking terms remain in the Lagrangian 
as the quark mass parameters transform as rrii — > exp(—iai)mi. The complex mass 
term can be factored into one of the quarks, |31j. 

M, = ['■;: ' I -{^ ' I (3.30) 




This results in the quark-mass phase oti being added to the 9 parameter, 

9 = 9 + NiOti = 9 + argdet My. (3.31) 
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Attempts to rotate away 9 by further transformations as in (3.27) will result in 



complex mass, and CP symmetry violating terms in the Lagrangian. 

The neutron's electric dipole moment (hence nEDM) is a measure of the separation 
of the centers of negative and positive charge within the neutron and should be a 
consequence of the CP-violating 9 term in the effective Lagrangian when evident 
as a complex quark mass term: [28J. 



^■CP-viol — i9m q 



' 75 _ 



(3.32) 



The nEDM can be computed by considering this 9 term as proportional to a 
one-loop correction in nn-meson coupling, |27J and a relationship of the following 
arrived at, [28J: 

d n ~ 2.7 - 5.2 x l(r 16 £, (3.33) 

with the range depending on precise couplings considered. [2S] in 2006 concludes 
that the phenomenological accuracy puts the nEDM, d n < 2.9 x 10~ 26 e cm. 
The measurement method compares the Larmor frequency of the neutron spin 
polarisation in applied electric and magnetic field when E and B are parallel and 
anti parallel. Thus the term 9 is limited by the lack of observational evidence 
of the nEDM to be of the order 9 < 10~ 10 . 9 is a phase originating as the sum 
of two unrelated terms {9qcd and the electroweak-QCD interaction related term 
argdetMjj). Having period 2tt it could feasibly take any value from ~ tt. Why 
it should be so close to is the CP problem. The solution pertinent to this thesis 
is a theorized additional U(l) symmetry and scalar field termed the axion. For 
completeness, several alternatives are outlined, [27], [28], [31]. 



Massless quarks: As pointed out above, in a massless quark model, one 
may rotate away with a chiral transformation the CP breaking 9 term to 



eliminate (3.24), implying that in this case 9 is not a physical parameter 



within the theory. [27] notes that it is sufficient that the mass of the up quark 



vanish (as evident in (3.30)), but also notes that Weinberg's up/down mass 



quark ration Z = ^ = | has historically ruled this out, and more recently 



[32J showed in 2003 through lattice calculations that Z = 0.410 ± 0.036. 
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• Spontaneous CP breaking : [25] postulates that the CP symmetry which 
is theorized to be broken by the 9 term is actually spontaneously broken. At 
the bare Lagrangian level one may set 9 = 0. However the same source notes 
that the 9 CP symmetry breaking term reappears at the one-loop level and 
complex Higgs vacuum expectation values are needed to set the quantized 
CP breaking terms to zero. 



3.2 The Axion 
3.2.1 Axion Models 

The leading candidate for solving the CP problem is the axion. The idea was 
first put forward in two papers by Peccei and Quinn in 1977, [33]. The theory 
proposes a new global U(l) "PQ" symmetry for the standard model (later termed 
U(1)pq with phase a'). [55J showed that the condition for CP conservation in 



eq. (3.34) below, i.e. 9 + a = 0, could be naturally achieved in the quantized 
Lagrangian, as the effective potential is minimized. The mass of the up quark is 
rotated back to the real plane (achieving CP symmetry) and there is no need to 
set m u = 0. U(1)pq is spontaneously broken and a Goldstone boson produced 
from one of the Higgs degrees of freedom. Here the CP-breaking phase is 9 and 
the phase associated with U(1)pq is a. 

M=h 0\ „\ (3M) 

\0 m d ) \ m d ) \0 m d ) 

Peccei termed this spinless scalar field the axion a(x) such that when U(1)pq 
is broken at energy scale f a (known as the scale factor, or decay constant) it is 
transformed as follows: 

a(x) — >■ a(x) + af a . (3.35) 

In effect, when added to it, the axion promotes the phase 9 (which is arbitrary) to a 
dynamical parameter, or equally, the axion as a dynamical phase can be redefined 
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to absorb 9. [35] notes that expressed in terms of the chiral anomaly (3.7), the 
invariant effective Lagrangian has the following 9 and a-containing terms (with 
axion interactions not included here). 



4?,a = A d ^ a + ^lf^ Tr GWGW + ij^TvG, u G, vl (3.36) 



where a kinetic term for the axion field has been added by hand. We note the 
negative sign in front of this kinetic term as a convention deployed in [35]. £ is the 



U(1)pq chiral anomaly co-efficient defined by (3.7): 

9 2 

®tijpQ = ~ ^ 167T 2 ^ r ^w^^vi (3.37) 

which appears when the U(1)pq symmetry is explicitly broken by QCD instanton 
effects, [31]. Mass acquisition by the axion provides a natural mechanism for the 
minimization of the now-dynamical 9 (as explained below). This initial model 
triggered a search for the axion and many variations of it, most of which involve 
the axion acquiring mass via coupling to other fields within the standard model. 
They key axion models are outlined below. 



3.2.1.1 Peccei-Quinn-Wilcek- Weinberg Theory 

This model arrived soon after Peccei-Quinn's initial papers following input from 
Weinberg and Wilcek, [33j . In addition to the standard model Higgs doublet, the 
PQWW model proposes an additional Higgs doublet, in which one field, 4> u couples 
to up quarks and the other, 0^ to down quarks with no cross coupling. The model 
requires that the quarks acquire their mass from the neutral components of the 
new Higgs fields, 0° and 0°. 

CpQWW-m = y u ULi4>uUR + ydd L <P d d R + h -c, (3.38) 
where m u = y u v u . The potential of the model is, [26] : 



U(4> u ,4> d ) = -fyl^u + h.c. 



(3.39) 
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With the U(1)pq symmetry, the Higgs, N quark fields (u and d) and 9 parameter 
transform as: 



u , <f) d ->■ exp(z2a u )0 u , exp(i2a d )0 d 
it, d — )■ exp(— za M 7 5 )w, exp(— 20^75 )d 
9 ->■ 0- iV(a u + a d ). 



(3.40) 
(3.41) 
(3.42) 



When electroweak symmetry is spontaneously broken the neutral Higgs compo- 
nents 4>° u d acquire vacuum expectation values (0° d ) and hence Nambu Goldstone 
fields. 

,.P,. 



v u exp^ 
exp(z 



(3.43) 
(3.44) 



In PQWW model, a linear combination of the two fields resulting from the 
neutral components of the Higgs fields results in the Z boson, while its orthogonal 
combination results in the axion field, a. 



sm(/3 v )P u + cos(/3 v )P d , 



(3.45) 



where (3 V is the angle between P u and Pa, the Goldstone fields. The axion couples 
to the quark fields resulting in complex quark mass which via the transformations 



in (3.40) can be transferred to 9 to give: 

9^9 



v v d v u > 



(3.46) 



a change which can be absorbed by a redefinition of a. [2S] notes that non- 
perturbative QCD effects explicitly break U(1)pq and result in the axion anomaly 



term (third term on the right-hand-side of (3.36)). [35] notes that it may be 



regarded as an effective potential for the axion 

U(a) ef f 



S £ 1e j ±L ^ fll/^ fjLU • 



(3.47) 



The addition of the axion to the Lagrangian (3.36) allows the promotion of 9 to a 



dynamic variable (9 + £ j-). This will have a minimum when the expectation value 
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of the axion field < a >= —^0. The fact that a is a dynamic variable provides a 
natural means for this potential and the phase (0+£j-) to be minimized, addressing 
the CP problem, where the < ... > is the vacuum expectation value operator. 



< ^sr- >= -7.4? < > i<— *» = °- (3 ' 48) 



Differentiating again with respect to the axion field, |35j provides an expression 
for the mass squared matrix for the axion: 

M 2 , d 2 U{a) eff f g 2 d ~ 

M a =< >= -J a ^Q- a < TrG, u G, v > \ <a>= _ re . (3.49) 

[25] and [37] show that due to inherent difficulties in computing low energy effective 
QCD quantities, axion mass computations are more practical if the axion degrees 



of freedom in (3.47) are transferred into effective interactions of the axion with 
QCD (the 7r and rj mesons), and the term quadratic in a equated to the mass. The 
PQWW axion mass has the following form (where m n and f n are the mass and 
scale factor of the 7r and v is the electroweak energy scale, ~ 250GeV), 



77V fir \A«,,W2,J 

ma-pqww = W V = 25KeV. 3.50 

v m u + m d 

The PQWW model and its variations are firmly linked to the electroweak scale 
as the field a is coupled with the Z boson. Experimental evidence soon ruled this 
set of models out, but the axion dynamics remain valid in general for subsequent 
models, which are outlined qualitatively below. 



3.2.1.2 Invisible axion models 

The PQWW model assumes that the U(1)pq symmetry breaks at the electroweak 
scale v resulting in a relatively heavy, coupled axion which was not found. If 
f a »v the axions are light (m a ~ j), weakly coupled and invisible. Making f a a 
free parameter allows the axion to be a candidate for cosmological phenomena: cold 
dark matter and to a lesser extent dark energy. In the Kim-Shifman-Vainshtein- 
Zakharov (KSVZ) model [3H], the axion is the phase of a new electroweak singlet 
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scalar field and couples only to a heavier quark, Qh with interactions of the type 
—hQhL<t>QhR — h*Q hR (fi*QhL, where <p is an electroweak scalar field singlet and h a 



coupling. As with all axion models, a chiral anomaly term originating from (3.37) 
arises. In KSVZ, rather than coupling directly to the (u and d) quarks (as in 
PQWW), it couples to a heavier quark and the axion couplings are then induced 
by the interactions of this heavier quark with other fields. 



The Dine-Fischler-Srednicki-Zhitnitsky model, [39], like the PQWW model re- 
quires a doublet of two non-standard model complex Higgs scalars. Like the 
KSVZ model it also has an electroweak scalar singlet which transforms under the 
U (1) pq symmetry and whose phase results in the dynamical axion field. This axion 
field then couples with the Higgs doublet and the complex degrees of freedom are 
transformed to the chiral anomaly term as in the PQWW model. The two invisible 
axion models share similarities. Firstly, they contain an electroweak (SU(2) x U(l)) 
scalar singlet which spontaneously breaks the U(1)pq symmetry at some arbitrary 
energy scale f a »v with the axion degree of freedom resulting from the phase 
a. QCD instanton effects explicitly break U(1)pq at some energy scale \i less than 

2 ~ 

f a resulting in a chiral anomaly term of the form a^^Tr G^ U G^ U which can be 
regarded as a potential for the axion field, [25] which minimises to eliminate the 
CP breaking 9 term. Crucially, the axion may acquire mass via direct coupling to 
heavy particles other than the light quarks at an energy scale less than fi. While 
not derived in this paper a result is quoted from [35] : 



_ m^U ym u m d 

m a -KvsF — — 7 ; — x 

f a m u + m d 



10 6 GeV 



eV, (3.51) 



fa 

where the lack of dependence on v, the electroweak energy scale is noted. The 
DFSZ axion mass has a similar form (not quoted here). f a is a free parameter in 



both and (3.51 ) may be expressed in the general form for an invisible axion model: 

m a = ^, (3.52) 

J a 

where \i is an energy scale related to QCD confinement Aq C d ~ 0.2GeV. This 
can set bounds for the value of f a via m a . These can be tested by considering 
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the interactions a QCD axion is likely to have and the resulting cosmological 
implications of these. 



3.2.1.3 String Axions 

String axion models arise from string compactifications generating PQ symmetries 
which can be spontaneously broken. These involve natural origins for the PQ 
symmetry unlike in non-string axion models. Model-independent string axions, 
[33"] . arise from the antisymmetric tensor field of the bosonic and heterotic string 



theories (we consider this axion later in the thesis, eq (5.70)). The properties of 
the string axion do not heavily depend on the details of the compactification. [4_3] 
computes the theoretical value of / as given by: 

/ = y/2—Mp ~ 10 w GeV, (3.53) 

where Mp is the reduced Planck mass, and ocu is proportional to the square of 
the unified gauge coupling of ten-dimensional superstring theory compactified to 
four dimensions. The value of / arises from the theory. The mass acquisition 
scale /i is a free parameter. [33] notes that in order for these models to address 
the CP problem, QCD instantons must be the dominant form of mass acquisition, 
thus /i ~ 0.2GeV. The authors note that higher energy scale instantons could 
also play a role. Model-dependent string axions arise from the zero modes of the 
antisymmetric tensor field, [43J. The / values are more variable in these models 
with a typical value of / ~ 10 l7 GeV noted, but with / ~ 10 15 GeV possible with 
fine tuning of the string action parameters. As with the model independent string 
axions, /i, the energy scale where the U(1)pq symmetry is explicitly broken so that 
the axion acquires mass, can be a free parameter of the model, depending on the 
energy scale of instantons responsible for explicit symmetry breaking. 
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3.2.1.4 Dark energy and axions 

This brief review of dark energy is sourced from [33], [35], [3Z] and where noted. 
Following the discovery of the acceleration of the expansion of the universe in 1998, 
dark energy in the form of a homogeneous energy density, contributing almost three 
quarters of the universe's mass-energy, permeating all space and exerting a negative 
pressure was postulated. A key revelation for theoretical physics of the newly 
observed phenomenon was that it appeared as if particle physics developments 
of the early universe were effecting current-era cosmology. [H] and [32] note a 
required density for dark energy in the current era of Ude ~ (10~ 3 eV) A to fit with 
observations of the known mass of the universe and acceleration of the most distant 
objects. A positive cosmological constant, interpreted as a universal vacuum 
energy, of = 0.7 was immediately proposed as the simplest explanation of the 
observed accelerating expansion. The source of this vacuum energy in the light of 
the much larger fundamental energy levels associated with quantum theory remains 
unclear. A vacuum energy originating from quantum theory, if one considers 
very early universe energy levels close to the Planck scale would be 120 orders 
of magnitude higher than the this required level. One may, as in renormalisation, 
introduce counter terms to cancel the high vacuum energies but this requires ad- 
hoc fine tuning. These and other difficulties have led to a second class of theories 
being proposed grouped under the term " dynamical scalar field models" , the most 
well known of which are quintessence models, although others include tachyon 
fields and dilatonic dark energy. [3H] notes that current observations are unable 
to rule in favour of either a cosmological constant or dynamical scalar fields (or 
another form of theory) as the cause. Quintessence models currently are the most 
favoured in theoretical dark energy research. They are represented by an scalar 
field coupled to gravity with a potential U(<f>) which may explain the dynamical 
aspects of dark energy and perhaps other dynamical aspects of the A-CDM model 
of cosmology, such as inflation, [4"7] . 




(3.54) 
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Initial quintessence models used a potential such as: 




(3.55) 



with a a positive parameter, and M an energy scale. The use of energy scales 
observed in particle physics such as M = lGeV can result in the required 
energy density of the field without the need for further fine tuning, and provide 
dynamical variation. We do not attempt to describe the details of the full range 
of quintessence models here, but focus on axion-related dark energy theories. 



Axion based quintessence models have emerged in the last decade. The 
quintessence axion model [H], (1999, 2000) uses four new pseudo scalar Goldstone 
bosons created by additional U(l) symmetries. Two of these relate to axions, 
the other two make contact with hidden sector quarks to provide mass to the 
axions. The two axions, f g ~ A p i an k and f a ~ 10 12 GeV, describe quintessence 
and the conventional CDM-QCD axion respectively. Mass acquisition occurs at 
/i ~ Aqcd for the f a axion and at fi ~ 10~ 12 GeV for the f q quintessence axion. 
The latter results in an ultra light mass of m ~ 10~ 33 GeV which is equated with 
quintessence. The mechanism of the explicit symmetry breaking is via the two 
additional bosons which have hidden sector interaction at the intermediate SUSY 
scale (10 l3 GeV) and electroweak scale (10 2 GeV) respectively. The quintaxion [50] 
(2002, 2009) builds from the quintessence model and seeks the qualities of a very 
large / value and a slow roll of the potential to current times. It also relies on 
a number of pseudo scalar Goldstone bosons, three in this case. Two of these 
represent invisible axions and one is a model-independent string axion, aui = o, q 
and the other a composite axion a comp which is the QCD axion with / ~ 10 12 GeV 
and \i ~ Aqcd- o>mi is the quintaxion with / ~ Api anc k and m < 10~ 32 GeV and 
slow roll potential U ~ A(^), where A ~ flrn 2 ~ 10 _1 GeV '. 

Several variations of the so called false vaccua theory postulate that the axion 
field does not correspond to its true value, and this false vacuum can act as dark 



3.2 The Axion 



59 



energy provided its lifetime is longer than the age of the universe. [53] suggests an 
"unstable axion quintessence" model in which the minimum of the axion potential 
is negative. 

3.2.1.5 Heavy axion models 

The heavy axion model in [51] is motivated by the lack of observational evidence for 
the axion mass in the ranges predicted by the invisible axion and suggests axion 
physics could fall under a superstring force, dubbed "QC'D" which operates in 
parallel and with similar properties to QCD. This model puts a lower bound on m a 
of ~ 5 x 10 2 GeV, with the explicit symmetry breaking scale Aqc'd = yU ~ 3GeV 
and / ~ 2 x 10 2 GeV. [55] (1997) builds on this idea. A toy GUT model with 
S(U)5x S(U)5 gauge symmetry is considered with the second SU(5) a mirror of the 
first, but breaking at lower energies and resulting in a QCD scale A mirror > Aq C d- 
The axion acquires mass from mirror interactions and is equivalent to a PQWW 
model in this mirror sector. The mass of the axion in this mirror PQWW-like 
model is given by 



where v and v m are the Higgs VEV in the standard model and mirror sectors and 
m aPQww is the axion mass as calculated by the PQWW model. A m a < 10 3 GeV, 
with A mirror < 10 5 GeV is proposed. The value of / in this instance we have 
/ ~ 10 7 GeV. The model in [56], (1993) uses as its basis a "Walking Technicolor" 
model, which results in sextet quark-axion state, the tjq, which the author notes 
has the properties of a conventional Pecci-Quinn axion, but with higher color 
instantons providing additional mass contributions. It suggests a m a ~ 60GeV 
with \i ~ 10 2 GeV suggesting a value of / also of this order, although a value for / 
is not specifically referred to. In [58J, (1992), as with other heavy axion models [56], 
[57] . breaking of the PQ symmetry occurs at just above the electroweak scale such 
that / ~ 2 x 10 2 GeV. Also as with other heavy axion models this uses non-Higgs 
EW symmetry breaking, in this case via a heavy top quark and includes four quark 



3 




(3.56) 
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flavours in total. A value of \i ~ / ~ 2 x 10 2 GeV is used and a m a ~ l^GeV is 
computed. 

3.2.2 Axion Physics 
3.2.2.1 Axion interactions 

We focus here on invisible axions which involve an electroweak singlet and thus 
experience the electromagnetic and weak as well as the strong nuclear forces and 
exhibit, depending on model, relevant interactions in the 9 Lagrangian. In a 
comprehensive review of axion physics [27] a generalisation of the ^-containing 
terms in the Lagrangian is presented, (here the overline notation previously used 



is dropped so that the 9 of section (3.2.1) is now taken as simply 9 and where it 



is understood it is now dynamic and incorporates the pseudo-scalar axion field a 
such that we may set 9 = j- in the following descriptions). 

£e = \fld"9d,9 - ^G^CT + {q L iD qL + q R iD qR ) 

ci(<9 M 6>)g7 M 75<? - (q L mq R exp(ic 2 9) + h.c.) (3.57) 
9 ~ 9 ~ 

Qj j^g^_2 ^-^G fj,t/G fii/ ~\~ CQjj-^^—^'DlFhi/Fhi/ -\- H'leptons ,8 

where the term c\ is the coupling of the interaction term derivative in 9, c 2 is the 
phase in the quark mass matrix, C3 is the coupling in the CP symmetry restoring 
term, and the coupling of an electromagnetic anomaly term analogous to 



(3.37) with F^y the electromagnetic field strength tensor. Ci ep tons,e contains axion 



interactions with leptons. ci, c 2 and C3 are couplings below the scale f a and the 



quark mass matrix m is real. (3.57) is constructed as a general expression and 
by assigning a non-zero values to combinations of Ci, c 2 and C3, well known axion 
models result (e.g. the PQWW axion is given by c\ = 0, c 2 7^ and c 3 = 0, and 
the KSVZ axion [38] by c\ = 0, c 2 = and C3 7^ 0). [27] notes there are also 
axion couplings to the electroweak bosons of the form 9WW and 9ZZ which are 
not shown here. It is noted that through 9 = j- all axion couplings are relatively 

J a 
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weak, suppressed by a large f a . [27] notes that the CP symmetry restoring term 
can be represented by a three quark instanton diagram of the kind first suggested 



by 't Hooft [15J. Couplings are represented graphically, (Figure 3.1). 




Figure 3.1: Axion couplings represented by Feynman diagrams, with G the gluon, 7 the 
photon, q the quark and I leptons, [27 1. 



Hadron couplings: These are the basis of low-energy laboratory based 
axion detection efforts and axion physics in supernovae. They are governed 



by Ci and C2 containing terms in (3.57). 

Photon couplings: pQ (19.2) describes the one loop QED chiral anomaly 
dfj,j^ 5 = — jf^TrF^F^ where e is the electric charge and F M „ the 
electromagnetic field tensor. In a manner analogous to the derivation of the 
axion-gluon coupling term, a axion-photon coupling term arises, as noted in 



03.57) : 

9 ~ 

■£#77 = C07777P; 9 Tri^MjyFui/. (3.58) 

OZ7T 

This may further be expressed in terms of the axion field a, and the electric 
and magnetic fields E and B, [26J: 
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with: 

c a77 = — Ti {3.60) 

with c 7 containing the model-dependence and a the electromagnetic fine 
structure constant. For the DFSZ invisible axion model, c 7 = 0.36 and for 
KSVZ, c 7 = —0.97, [26]. This coupling results in the a — )• 77 decay, which 
estimates has a lifetime of: 

r(a— + 77) = 2^-3, (3-61) 

c a77 a em 'n2 a 

and notes that for c a77 ~ 1, an axion with m a = 24eV has a lifetime of the 
order of the age of the universe Tu ~ 4.35 x 10 17 s. 



3.2.2.2 Effective Potential of the Axion Field 



In this research we utilize without formal derivation the commonly quoted form 
(e.g. eq 29, [26J) for the effective potential of the axion field where m is a mass 
scale: 



U(a 



eff 



m (1 — cos 



fa 



(3.62) 



In (3.57),|27J provides an expression for the 9 dependence in the axion's Lagrangian, 



and notes that the cosine effective potential of the form (3.62) is determined by two 



fundamental properties of the axion - periodicity with period 2irf a and minima at 
a = (0, 2nnf a ..), with n an integer. However, we note that (3.57) does not represent 



an effective Lagrangian for all the 9 terms, for example, the 9 photon and lepton 
coupling and kinetic terms have been added by hand. [27] notes that while the 
cosine form for the effective potential may be a simplification the vast majority of 



cosmological axion models use the initial cosine term as in (3.62). 



The history of computations of the effective axion potential is reviewed briefly. 
Prior to the initial suggestion of a new dynamical axion field in 1977, [33], 't 
Hooft, pj)] and Weinberg, [23] described the non-perturbative nature of the QCD 
vacuum. Based on this, Peccei-Quinn began with the following Lagrangian, where 
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g' and h are coupling constants and ip and <fi are fermion (quark) and Higgs scalar 
fields respectively. 



C = -\G%G a ^ + nPD^ + iP 



+ # 



2 

^ | 2 -fx 2 \<f)\ 2 -h\<f)\ A (3.63) 



Peccei-Quinn find the effective potential of this axion model can be stated to a good 
approximation to leading order in g' and h and g'X, where 0's vacuum expectation 
value is defined by (0) = Aexpz/3, with (3 and A being real constants. 

V e (<t>) = U{<j>) -K\gf<f>\ cos 9, (3.64) 

where K is a real and positive constant, and 9 is the dynamical axion phase. [33] 
notes the limitations on this result in that g' and h and g'X are to leading order 
and the region of validity being that these three constants be small. However the 
authors note that for constant and using the dilute gas approximation (where 
a dilute gas of instantons is approximated by considering a superposition of one- 



instanton solutions at great distances from one another) (3.64) is valid to all orders 
of g'X. In recent research [53] the authors note that exact analytical solutions to 
forms for the the QCD axion effective potential are limited by the presence in 
the action of strongly coupled terms, while numerical methods are limited by 
the imaginary nature of the action. They strive to demonstrate using the so- 
called "interacting-instanton-liquid- model" a cutoff independent formulation for 
the axion potential and mass. It is noted that difficulties with these calculations 
arise at low energy levels. [59], (2002) presents an analysis of the form of the 
axion effective potential in light of new understandings of the QCD vacuum using 
supersymmetric gauge and brane theories. The authors note that the potential 



as in (3.62) may have higher cosine powers and in general is a smooth periodic 
function of 9 with period 2ir. They conclude that with axion models (including 
invisible axion models) using light quarks such that m q <C ^QCD-cutoff (where 
A-QCD-cutoff is the cut off used to integrate out the QCD degrees of freedom), the 
form (3.62) is valid. We make comments on our use of this form in our research in 
section 
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3.2.3 Axion Phenomenology 
3.2.3.1 Axion Production 



Section 3.2.1 describes how the axion field arises as a result of spontaneous breaking 
of the U(1)pq symmetry at some energy scale f a and then acquires mass close to 
A.qcd- The production of massive axions from this method is known as vacuum 
realignment, [26J, with the axions characterised as non-relativistic and largely non- 
interacting. If f a » Kqcd-i as in invisible axion models, the axion may be 
considered in the context of the ACDM model of the universe as a candidate for 
cold dark matter. Considering the cooling of the universe following the big bang at 
To the key temperature milestones are inflation reheating at Tp, the temperature at 
which the U(1)pq is broken at Tpq and the QCD scale temperature at Tqcd when 
the axion acquires significant mass. [2E] notes two additional plausible methods 
of axion production from the vacuum which may contribute to an axion density 
with the condition Tp > Tpq. These are axion strings and domain wall decay. 
We comment only on the former here for brevity. Axion strings [H], [60J, are 
theorised to arise from the U(l) pq breaking as topological defects in a similar 
fashion to cosmic strings arising as a result of a global U(l) symmetry breaking. 
Axions are produced in the string oscillations until the strings decay. This decay 
continues until the axions acquire mass at the QCD temperature scale [62] , [61]. 
If Tp > Tpq, the axion field is not homogenized by inflation and axion strings are 
produced from Tpq. 

We here consider only vacuum realignment as a production mechanism. At the 
Tqcd the axion acquires a temperature dependent mass m a and an effective 
potential given by [26] , [62] , [63] : 

U eff = m 2 a (T)f 2 a (l-cos6). (3.65) 

If f a is very large and with suitable choice of a(x) the oscillations can be damped 
and the potential varies only with m 2 a {T). [26] considers the solution to 9{x) = 
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a(x)/ f a equation of motion in the FRW metric. 

6 + 3H(t)9 - -^V 2 ^ + m 2 a T(t) sm(9) = 0, (3.66) 

where T(t) is the time dependent temperature, S(t) is the scale factor, H(t) the 
Hubble parameter equal to | and dot is derivative with respect to time t. For 



T R > T P Q, [26J derives a value of f a ~ \Q l2 GeV which from (3.51) results in 
a lower bound for m a ~ QfieV. Other models and scenarios for early universe 
axion production [26], [62], [75], [65] present a range of bounds. In axion models 
based on U(1)pq symmetry breaking occurring after inflation (or in inflation-less 
cosmological models) estimates place f a in the range 10 s ~ 10 12 GeV, with the 
upper limit of f a constrained by existing levels of cold dark matter observed. If 
Tr < Tpq estimates of the value of m a range from meV to neV levels (resulting 
in f a in the range 10 9 ~ 10 15 CreV), thus raising the energy scale of f a to pre- 
inflationary levels (A/tvfl ~ 10 14 GeV, |67j). Given an existing density of axions 
remaining from vacuum realignment, jSH] notes that as well as the axion photon 



conversion, the following processes (based on couplings in (3.57)) can result in 
axion emission in astrophysical objects (stars): 

• hadron-hadron bremsstrahlung-type interactions: H + H—tH + H + a; 

• photon-electron interactions: 7 + e _ — > e~ + a; 

• electron- nucleus bremsstrahlung: e~ + N — e~ + N + a; 

• photon fusion: 7 + 7 — > a. 

Observations of stellar energy loss result in a lower bound on the invisible axion 
of f a ~ W 9 GeV. 



3.2.3.2 Axion Detection 



In the original PQWW scheme for axion production, the mass would have been of 
order lOOkeV, [25], and thus within the bounds of laboratory testing. Experimental 



3.2 The Axion 



66 



studies of the PQWW axions have consisted of examining radioactive decay 
and results to date [26] have not detected an axion-like particle at this mass. 
Subsequent axion detection has focused on the search for cold dark matter- 
motivated invisible axions. Detection efforts have centered around the axion- 
producing version of the Primakoff effect. It is predicted that axions with 
energies of a few keV may be produced in the interaction a < — > 7 + 7 in the 
presence of electric and magnetic fields within the solar plasma. Helioscope 
experiments consist of dipole magnets orientated towards the sun to catch the solar 
axions which would be converted to photons in the form of x-rays. The largest 
helioscope experiment is the CERN Axion Solar Telescope (CAST) project [68] 
operational since 2003. Initial results provided an upper limit for the coupling g a y 
of 8.8 x 10 _11 GeV^ _1 applying to an axion mass of m a < 0.02eV. Phase II results 
increased the sensitivity of results to m a < OAeV, with sensitivity of m a < 1.2eV 
predicted by 2011, [58] . [69] . There are several searches underway specifically aimed 
at axions being responsible in part for cold dark matter. The Cryogenic Dark 
Matter Search (CDMS), [71] looks for possible solar-axion conversions to photos 
or galactic-axion conversions to electrons within germanium crystal detectors. The 
Axion Dark Matter Experiment (ADMX), [70], uses a microwave cavity detector 
to search for CDM axions in the Milky Way galactic halo. The mass detection 
zone for ADMX is 1.9 < m a < 3.4/xeV. In [72], the authors refer to the theory 
that axion-photon interactions can reduce the attenuation of very high energy 
gamma rays (> lOOGeV) traveling over cosmological distances, but such effects 
have not been observed within the Fermi gamma ray experiment. [66] summarises 
the recent laboratory, astrophysical and cosmological limits placed on the invisible 
axion's mass and / values and regards the most conclusive evidence of axion mass 
to come from stellar data. Restrictive limits arise from observed neutrino signals 
from the SN 1987A Supernovae. If the axions are above the m a ~ lOeV level, it 
is argued that they would have been observed in the Cherenkov detectors used to 
count the neutrino output. In a 2010 survey of data from the Fermi Gamma Ray 
Space Telescope, looking at the high energy spectra from two specific sources, the 
authors conclude that there is no evidence for an axion-like particle attenuation 
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effect. To sum up, while axions have not yet been conclusively observed, Figure 
3.2| cites the range in which they cannot be ruled out: O.lmeV < m a < lOmeV 
corresponding to a range for / of 10 12 > / > 10 9 GeV. 
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Figure 3.2: The bounds of the scale factor f a and mass with axion production and 
detection methods, [27]. 



Chapter 4 

Full Quantization of the Axion 



4.1 Introduction 



The axion is proposed as an additional degree of freedom necessary to explain the 
lack of CP symmetry breaking observed when the action based on the bare QCD 



Lagrangian (3.1 ) is quantized resulting in a one loop anomaly term (3.24). The full 



axion Lagrangian is illustrated in (3.57). We note that the kinetic and non-QCD 
interaction terms for the axion in this have been added by hand. Full quantization 
of the axion should then realistically involve quantization of the axion degree of 
freedom. There have been some recent attempts to do this (see for example |74] ) . 
in axion models and we propose to further this in this research to shed new light 
on fundamental physics arising from the axion in the early universe. Finally, it is 
acknowledged that a large portion of the derivations in this section were published 
in [7U], (2010) by Alexandre and Tanner. 
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4.2 Quantization of the Axion 

We begin with an action in which the non-axion degrees of freedom have been 
integrated out. We thus assume that for some general action containing axion and 
non-axion fields, the following operation has been performed. 



Z[J}= V[(f>}eMSM- s em, 



(4.1) 



where <fi represent the non-axion fields (gauge bosons, leptons, hadrons) and S[9] 



is the 9 (axion) degree of freedom of the action. Section 3.2.2.2 describes attempts 



made at such computations and general results obtained. We do not attempt 
such detailed computations in this research and denote the resulting 9 dependent 
potential as U (9) with the conditions only that it is periodic in 9 with period 2ir 
and U(9) = at 9 = 0, 2n, ....2nir, [27]. We consider this our bare action in a 
general form. 



So 



d x 



oo 



a^e + y a n (i-(cos 



2 " ' ^ 

n=l 




(4.2) 



This is based on discussions outlined in Section 13.2.2.21 in which the most 
generalised form of the axion potential is represented by a cosine power series, 
Full quantization will determine the coefficients a n and powers n of the 



cosine. (The axion potential expression used commonly (3.65) is returned if only 



a\ is non-zero). We place the / dependence in (4.2) in the the kinetic term and 
make several remarks about our quantization approach. (Note: in this chapter, 9 is 
the phase of the dynamic axion field defined in section (3.2.2.1), i.e. a(x) = f a 9(x), 
and also we hence drop the a subscript such that /a = /). 



We define an ultraviolet cutoff A to impose on (4.2), which will, along with 



/ will be a parameter of the theory. We note at this point that A is distinct 
from the cutoff used to integrate out the non-axion degrees of freedom in the 
operation conducted in (4.1 ). This we term ^QCD-cutoff and we assume that 



^QCD-cutoff < A. Further we point to recent research outlined in section 



(3.2.2.2) suggesting the form of the periodic axion potential is independent 
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of AqcD-cutoff an d we assume here that it is not a parameter of our theory. 
We further assume that the cutoff A defines the upper limit of the theory 
and as such / < A. 

In terms of the the normalized scalar axion degree of freedom a = f9 in 



(4.2) we wish to deploy non-perturbative methods derived in section 2.3.3 



Rather than an evolution with an energy scale k (as in the Wegner-Houghton 



equation, (8.21 )) we utilise another parameter of mass dimension, in this case 
/, to describe the exact evolution equation. Such an approach is used in |103] 
and |104j and in the string cosmology section of this thesis, in the context of 
the bosonic string with the evolution parameter being the string mass scale. 
While / is generally not considered a variable parameter within conventional 
axion theory, it is not fixed in a physical sense, and we use the evolution of / 
as a mathematical technique to arrive at our evolution equation. With similar 
logic to |102j , as / — > oo the theory describing the axion represents the bare 
theory as the kinetic term is large compared to the potential term. As / 
decreases, quantum fluctuations begin dressing the system. Thus evolution 



with / in an exact equation similar to (2.55) should serve our purpose of 



capturing all quantum effects in an effective expression. 



We utilise results on effective field theory outlined in Section 2.1 for our 



quantization. In this, the starting point is the generating functional, Z[J] 



which defines the full quantum theory in path integral form (2.1). Then 
via the definition of the connected graph generating functional Z[J] = 
exp(— W[J]) we arrive at an expression for the effective action ( 2.9[ ) which 



is expressed in terms of the classical field <p c i. We want to deploy the same 
notation and results for the scalar axion field a(x) = fO(x). We state the 
partition function (or generating functional), Z[J] in Euclidean 4-D space 
time. 

Z\J\ = J D6exp- \ S[6] + J d D xJ( y x)9( y x)^j = exp-W[J]. (4.3) 
We note that 9, while a phase periodic in 2tt, is also a function of space time 
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and make the following assumption about computing (4.3). 



J D6 = IlJ o de ^ 



(4.4) 



4.3 Detailed Calculation 



Taking the derivative with respect to / (represented by an overdot) of (2.9) we 
obtain (with <p c i(x) replaced by 9 c i, the classical field for 9 for which we now drop 
the cl subscript and denote simply as 6): 



an 

d'x^J(x) = W[J], 



(4.5) 



following the methods used to arrive at (2.52). We have used 8W[J}dJ(x) = 6 as 
the definition of the classical field 9 and consider a constant field configuration so 



that If = 0. As was done in (|2.54|), combining (|4.2|), (|4.3|) and (|4.5|) we obtain: 



of 



T[6] = W[J] = -Z = f I <Tx< dffi&'O > . 

Zj 



(4.6) 



Which as derived in (2.55|) can be computed as: 

no] = f 



d 4 xd^6j 



Tr 



d d fS 2 T[9] 



d x d y \56J6 y 



(4.7) 



where Tr[....] = J d 4 xd 4: y[....]5 4: (x — y) and accounts for all quantum corrections 
up to the distances which are the space time equivalents of A, our regularization 
cutoff energy scale for the quantized axion theory. We now assume a generic form 
for T as follows, with / dependence also in U e ff(9): 

P 
2 



v[e] 



d x 



-d^9 + U eff (9) 



We consider an axion field constant in Euclidean space time, 6 = 9q. Using (2.16) 



and taking the partial derivative with respect to / of T[6] we have (with V the 
Euclidean space time volume): 



f [0] = VU ( 



eff- 



(4.9) 
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which we combine with (4.7) to get: 

VU eff = fTr 



d d f5 2 T[9] 



8 X dy \59J9y 



-i 



(4.10) 



We evaluate 



S 2 T[8] I 
66 x 50 y l o 



from (4.8) using a derivation similar to that performed in 



arriving at (8.23): 



6 2 T[9] 



59J9 y 



fd^]5\x-y), 



(4.11) 



'M = P"ff{d ) 

where the " refers to a second derivative with respect to the axion field 9 X . We 



evaluate the expression under the trace in (4.7) and convert to momentum space. 
We use: = -p 2 , F(p) = J d^xe' 1 ^ F(x) and / d 4 xe" ip - x 5 4 (0) = V/(2n) 4 . 



fTi 



d d f5 2 T[9] 



d x d y \ 59 x 59. 



\e<> 



, X UVy 



f / d 4 xd*y5\x - y) 



d d 



d x d y (U" 



f 2 d^) 



f / d 4 xd*y5\x - y) 



d 4 p d 4 q —p^q^ip + q) 
(2tt) 4 (2vr)M^o) + /V! 



-i(px+qy) 



Vf 



d A p 



(2tt) 4 (U» ff (9 ) + /V) 



(4.12) 



where we have equated 5 4 (0) in momentum space to the spacetime volume V. 
Equating this (4.12) with T as shown in (4.9) gives: 

fp 2 



U, 



ff 



dp 



(4.13) 



(270*17^(00) +/V 

We evaluate the integral by considering / = J dQ ' D ( ^)D dp = j^fo J p D ~ 1 dp 1 
where the solid angle in four dimensions ^4 = 2n 2 . We use the substitution x = p 2 
and J p 3 dpf(p 2 ) — of P 2 d(p 2 )f(p 2 ) = 2 / xf(x)dx to compute the integral (4.13), 
where A is our high energy cutoff, above which our theory is not defined. 



U, 



a. 



2(2tt) 4 
1 



A 2 



dx 



fx 2 



16tt 2 





A 4 
2/ 



xf 2 + U t 
(0 )A 2 



"ff( e °y 



U" 



P 



+ fS — ln 



/ 2 A 2 



1 + 



• (4-14) 
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4.4 Flattening of the Axion Potential 



We refer to our description of the flattening of any concave potential in Section 



2.4 In this it was argued that given the result expressed in (2.59), the shape of 



the effective potential must be convex, with concave features quickly suppressed 
by tachyon-mode fluctuations as the system undergoes quantization. We note that 



the result arrived at in (4.14) is generic and is in line with the derivations used to 



arrive at (2.59). We also point to our demand that the constraints on the effective 



potential in (4.8) are that it be periodic in 9 and U e ff{9) = at 9 = 0, 2ir, ....2mr, 
with n an integer. The periodicity and requirement that it be convex thus requires 
U e ff in our theory to be flat, or only a ^ in (4.8). We thus require U'Jjj(9) 
and (4.14) reduces to: 







U, 



A 4 



eff 



327T 2 /' 



(4.15) 



We make several comments on this result. We can interpret U e ff, the effective 



potential of the axion in (4.15), as a vacuum energy density associated with the 
axion at some point in its evolution. It has general solution: 



/ 



dU eff (f) 



A 1 



327T 2 



and if we take the boundary condition that U e ff 

A 4 



U e f f {f) =U A + 



32tt 2 



In 



k dJ_ 
/' 

- U\ when / = A: 

/ 
A 



(4.16) 



(4.17) 



We note that the condition (2.59) we have imposed on the solution to (4.14) applies 
to a non interacting scalar theory, such as spontaneous symmetry breaking in the 



double well potential as outlined in section 2.4 We note that interactions are 



certainly present in axion theory as detailed in (3.57) and Section 3.2.2.1, (It is 
also acknowledged that the axion coupling to the metric g^ u is not considered in 
this thesis). In a dynamical axion theory, these terms arise following spontaneous 
symmetry breaking at energy scale / and subsequent mass acquisition by the axion 
following contact with the QCD energy scale. We use the following qualitative logic 
in the interpretation and relationships between U\, the spinodal instability effect 



and the perturbative expansion of the quantized bare action (4.2). 
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The flattening of the potential induced by the spinodal instability is a 
quantum effect which by nature includes all quantum corrections in the 
resulting constant value of the effective potential, i.e. it is a tree level, exact, 
non-perturbative value. 

In our evaluation of U\ we are considering the potential at energy scale near 
the cutoff where / = A — > oo. In this regime, in the bare action represented 



by (4.2), (and also by our scalar field approximation of this, represented by 



the Lagrangian density introduced in the next paragraph, (4.19)) the bare 
action may be considered to be quadratic in the field 9 with the interactions 
negligible, and hence free. As such the second derivative of the bare potential 
is zero and the one loop correction is an exact expression of the full quantum 
theory, as with the spinodal instability effect. 

We thus interpret the flattening of U e ff as an effect which provides a 
contribution to the evolution of the effective potential in the early phase of 
axion development, prior to mass acquisition and significant coupling effects. 



We consider the axion scalar field expressed by 9 — a/f, and are interested in the 
evolution of this field near the cutoff when / — y A. Thus with small 9, we may 



express the Euclidean potential U{9) in ( |4.2[ ): 

oo 

J^a n (l - 



COS I 



n=l 



^e 2 - V + o[9 6 } 

2 4! 



ki f a. 2 



2<f 



4!7' 



(4.18) 



where k\ and /c2 are constants which are combinations of the factors a n in (4.2) 
(note the units of a n and hence k\ and &2 are quartic in mass units). To terms 
quartic in %, (4.18) resembles the double well potential in (2.33). We may consider 
the axion initially as arising from a U(1)pq symmetry which is spontaneously 
broken, a system whose classical Lagrangian for a scalar field = <p(x) = f9 is: 

A 



/2\2 



(4.19) 
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where we now equate ^ = jj 2 and |f = A. The one loop effective potential of such 



a Lagrangian can be evaluated by reference to the Euclidean version of (2.21) 



C/ e// (0) = f/(0) + ^Tr \n[-d 2 + C/"(0)], (4.20) 

where is the classical scalar field and U(<f>) is the classical potential and we thus 
have when A = /: 

U A = U{<j>) + ^Trln[-<9 2 + U"(<f>)]. (4.21) 
For completeness, we consider two cases, where \i = and A 7^ 0; and vice versa. 



At n = 0, A 7^ (i.e. a vanishing renormalised mass squared condition) the 
system can be considered on the verge of spontaneous symmetry breaking as 
/i becomes > 0. As such, it may be considered that the system is in the initial 



stages of the development of the axion field, discussed further in section (4.5). 



In the derivation leading to (2.29),[2J provides a one loop effective potential 



of (4.19). We state rather than derive the result here: 



V'ffM = ^ + j^y^ (ln^ - f ) + <D[\U, (4.22) 

where m is an arbitrary energy scale being considered and A m the energy- 
scale dependent effective coupling. Here is the classical field denoted by 
cj) d =< > elsewhere in this thesis. The term quadratic in A m is the first 



order correction. The flattening of the potential exhibited in (4.17) is a tree 



level effect and we identify the term U\ with the first order correction in 



(4.22). In our initial assumptions stated in the text following section (4.2) 
we consider a constant axion field configuration 9 = 9 . We further set a high 
energy limit to our theory of A and found that spinodal instability effects 
result in a flat effective potential. With these conditions in mind, we consider 
a constant and also that 0, m « A, thus satisfying our approximation in 



(4.18). We further now make an assumption that m ~ based on the fact 
that both are arbitrary for the purposes of our reasoning and both are small 
compared to A. We note we are considering the case where A = / giving 
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from eq. (4.17) U e ff = U\. We do not take the energy scale m (and therefore 



<ft) as equal to A as we wish to explore the behaviour at energy scales lower 
than the cut-off of our theory, A. We wish to keep m in our theory as a 
variable representing the energy scale below A at which we are investigating. 
We further note that in this logic we should include the ^-containing term 



in eq. (4.17) but we assume that the variation of / from A is not significant 
leaving the assumption U e ff = 



We now have the following. 



Ueff 

and with / = A: 



A 4 . ff\, 4 /l, 25 2 



U, 



eff,f=A 



m 



rA, 



25 



A 



(4.23) 



(4.24) 



4!'"" 96tt 2 ' 

In an analysis of the parameter A m , we note that the use of the classical 



potential associated with (4.19) with // = in quantization of the axion is 



similar to the approach taken in [71] where it is shown that invisible axions 
form a Bose-Einstein condensate and in and [73] where the axion field 
evolution is considered. In [71] (here the conventional, invisible QCD axion 



is considered as in section (3.2.1.2)) the authors use a 4 scalar model for 



the axion and compute the effective scalar coupling constant as follows: 

,2 



A 



m 2 a m 3 d + m 3 u 



0.35- 



(4.25) 



f 2 (m d + m u ) 3 p 1 
where m a is the axion mass and and m u are the up and down quark 
masses. The authors state that this formula is obtained by using current 
algebra methods to derive an expression for the axion effective potential and 



equating the fourth-order coefficient to A. While A m in (4.23) represents the 



scalar coupling at the energy scale m and A in (4.25) represents the coupling 



at the QCD energy scale, as a cursory approximation we take A = A m , which 
(with / = A) gives: 

4 

U effJ=A ~ 4 x 1(T 2 —ml (4.26) 
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which, as our energy scale m approaches the cut-off A = f reduces to 
~ 10~ 2 m 4 . Here we have assumed A 2 << A. The mass of the invisible axion 
has been experimentally reduced to a bound of 10 _4 eV < m a < 10 _1 eV, 



[26] and Figure 3.2 [TT] notes a phenomenologically required energy density 
for dark energy in of order Ube ~ (10 _3 e^) 4 , which is representative of 



commonly quoted values. While the result in (4.26) is limited in usefulness 



by the inconclusive nature of the ratio yt < 1, it could be that it does not 
vary by orders of magnitude with accepted values for Ude- 



The beta function for the coupling constant in <f) 4 scalar theory can be 
expressed ( 2.30[ ), (([I], 12.2), with m the energy scale: 

5X 3A 2 



/3(A) 



m- 



5m 167T 2 



0(X ; 



(4.27) 



The result (4.24) depends on X m and through this the energy scale m. Using 



(4.27), fine tuning of (4.26) is likely necessary to arrive at a more accurate 



residual axion energy density. 
We now consider the case of A 



0, \i 7^ 0: Here we consider the following 



bare potential, from (4.19): 



U{tf>) 



-J' 



2 J.2 



(4.28) 



where /x is a mass parameter. In we similarly consider the flattening of the 



potential as a tree level effect, the A = case for (]4.20|) may be solved directly. 

(4.29) 

(4.30) 



U A 



f/(0) + -Tr{ln [(-d 2 + U W)]} 



1 A 4 
2^ + 64^ 



We note \x is a mass scale associated with the bare axion potential and we 
make the assumption that this is at least several orders of magnitude less 
than the overall cutoff of our theory such that /i << A and we heglect the 



|/i 2 2 term in eq. (4.29). We thus have: 



Us 



—A- 

327T 2 \fi 



(4.31) 
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and: 



U, 



KM 1 - 

327T 2 \fi 



(4.32) 



Here we have taken 6 the classical field such that 6 < A and assume <t> A << 



A 4 in the evaluation of U\ above. If, as in (4.26) we refer to [H] where \i is 
equated to the axion mass, we obtain, with A a free parameter of the theory: 

A 4 



U. 



eff,f=A 



32tt 2 



In 



A 



(4.33) 



We provide an interpretation of these results in section (4.5). 



4.4.0.3 Axion Interactions 



The description of the flattening of any concave potential in a scalar effective 
potential when quantized and its resulting flattening of the scalar potential was 



discussed in section 2.4 and refers to an effective potential where the source J{4>) is 
zero, that is, a theory with only self interactions. As the axion makes contact with 
the QCD and electroweak energy scales a range of interactions emerge as outlined 
in section 3.2.2.1 and in (3.57). This means that J(4>) is not necessarily zero and 
the convexity condition (2.59) (that is — ^ < 0) may not hold and it may 
be avoided. For example we consider a non-QCD interaction, that of the axion 
with leptons: gfOip^ip. Integrating out the lepton degrees of freedom results 
in an expression for U e ff which contains the one loop correction to the lepton- 



axion interaction term. In this, for simplicity, we assume the form of (4.8) is 
if l — (cos#)), where a\ = m 4 , the infrared cutoff is k', the coupling is g and we 



m 



neglect the mass of the leptons and make use of (2.22). 



U eff = m 4 (l - (cos9)) - 



d^k 



v k 2 WfQ 

~k 2 



In 



2 f2n2-\ 



(4.34) 
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We evaluate the second derivative of U e ff at 9 = it where the m 4 (l — (cos 9)) form 
is concave. 

S 2 U eff (0) , = _ 4 ^ 2 / 2 fc /2 fc /2 + 3g 2 /V 

C^ 2 m 167T 2 fc' 2 + 3# 2 / 2 7r 2 

^ /2 (3, 2 / V)lnfl + -^y (4.35) 



167T 2 V g 2 f 2 7r 



For a high value of the cutoff, when k' >> gf the result (4.35) is less than zero 
and the spinodal instabilities arise, flattening the potential and it can be said 
that the lepton interactions add to the concave potential. If k' < gfir, a realistic 
assumption given the accepted value for / ~ 10 12 GeV and if k' ~ 10 3 GeV, the 
electroweak scale, the spinodal instability effects can be avoided. 



4.5 Summary and Discussion 



In section (4.4) we described how the spinodal instability in the axion's early stage 



quantized potential flattened it. We thus derived a non-perturbative expression 



for the evolution of the effective potential with / as in (4.15). We then assigned 
a boundary condition to this differential equation such that the potential when 
/ = A is some value U\ leading to: 

u °>'W>=&> iD {£) +u - {4M > 

Physically, at energy level A, which represents the upper limit cutoff of our theory, 
we consider that the evolution of the axion field in its earliest describable form. In 
this sense, for very small 9, the commonly used cosine form of the axion potential 



(3.62) may be equated to a double well scalar potential as in (4.18). We consider 
that such a double well potential represents the origin of the spontaneously broken 
U(1)pq symmetry responsible for the axion field. An initial admission on this 
analysis is it does not incorporate any time parameter describing evolution of the 
field, but refers qualitatively to sequences of configurations. We consider that our 



result (4.36) describes the system on the verge of spontaneous symmetry breaking. 
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In the double well scalar potential (4.19) this occurs when \x changes from negative 



(representing a symmetrical "U" shaped potential) to positive, which is a form 
where spontaneous symmetry breaking can occur. Thus we explore the \x = 
point as when the axion field is emerging. At this point the axion has yet to evolve 
fully and the effects of the spinodal instability serve to flatten the potential, prior 
to any interactions. We thus regard the potential as being represented by our 



result (4.36), the form of which is impacted by the spinodal instability which is a 



tree level quantum effect. We consider that the term U\ in (4.36) is a one loop 



correction to this. We identify it with the one loop correction to the quantized 



form of (4.19) at fi = 0, which we compute using an established result quoted in 
[2]. The remaining parameters of our form for U e ff(f), when we take A = /, is 



then A, the scalar coupling constant of the quantized form of (4.19) and an energy 



scale m used to compute the one-loop correction to the scalar field Lagrangian as 



in eq. (4.22). We rely on work [74] in which the authors compute (using current 
algebra methods) the effective scalar coupling constant in terms of / and the axion 
mass m a . For A = /, the result which can then be expressed in terms of m a and 

4 

the ratio of yr where m is arbitrary and less than A and / is the scale factor equal 
to A (thus the ratio is small). 



U, 



4 x ID" 2 ^mi, 



(4.37) 



where is an undetermined ratio less than one (but may be close to one as we take 
m — > A). We interpret this as an energy density associated with the axion at the 
earliest phase in its development, prior to significant interactions, and influenced 
by spinodal instability effects. (At later stages the axion acquires mass and resolves 
the CP problem as well as providing a candidate for cold dark matter). With the 

4 

caveat that the ratio jr is not determined by our theory, it could be that it is not 
orders of magnitude away from required values for dark energy (Ude ~ (10 _3 eV) 4 , 
jllj) . We make several further comments and refer to recent related research. 



The steps leading to (4.37) start with quantizing the dynamical field of the 



axion. We next consider the quantum effect of the spinodal instability on the 
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resultant field in its earliest phase of development. Finally we use and existing 
result for the quantization of the double well potential approximating the 
early phase axion field as a boundary condition to add numerical estimates 
to our result. Thus our theory is very much a full quantum theory of the 
axion. 



Section (3.2.1.4) outlined quintessence theories involving a dynamical scalar 
field changing in space time, as evidenced by the accelerating expansion of 
the universe in the current era. Key work in this area was conducted by 
Kim and Nilles, [48], [49]. This focuses on linking quintessence with an ultra 
low mass axion whose potential has "slow rolled" down to a level associated 
with the required dark energy value. The ultra low mass is obtained by 
considering a high /-valued string axion or an axion which acquires mass 
through contact with some hidden sector quark of ultra low mass. In contrast, 



our result (4.37) does not rely on any non-standard model physics other than 



the proposed QCD axion. 

The well cited works by the same authors in [63J and [64] describe a family 
of particles termed "pseudo-Nambu-Goldstone-bosons" (PNGB), of which 
the axion is an example. These particles exhibit spontaneously broken £7(1) 
symmetry at a scale / and further explicit symmetry breaking at a lower 
scale /i, and acquiring a mass ~ f^ 2 /f- [63] treats the neutrino as a PNGB, 
and attempts to link its dynamical field to an effective cosmological constant 
for several expansion times in the universe. We note that this approach links 
the mass of the neutrino-PNGB at certain eras to achieve required energy 
densities. [11] builds on the well-cited work in [63] but, as in our analysis, 
considers spinodal instability effects on the cosine form of the neutrino's 
effective potential resulting in a flat energy density of M 4 , where M is the 
mass of a light neutrino, corresponding to a dark energy like effect. In 
contrast our work, with the QCD invisible axion as the PNGB, results in 
m 4 being proportional to an energy density which we compare with dark 
energy. 
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We have characterized the early phase axion field as being flat due to spinodal 
instability effects. Quintessence models require a dynamical scalar field, [U] . 
We have discussed only qualitatively the evolution of the axion field. Further 



analysis of our result is necessary to determine how (4.37) could be shown 
to evolve into the current era. In terms of the axion mass, m a which is a 
free parameter of our theory, [53], for example, provides a discussion on how 
the axion mass may evolve with temperature scale T in the early universe 
evolution (the result m 2 a = a a A 4 / f%(T/A) n ) is quoted where a a = and 
n = 6.68). Additionally the parameter A can be considered as a running 
coupling whose evolution with energy scale is governed by a beta function as 
in gj7b. 



Chapter 5 

Background to String Cosmology 



Since its emergence in the 1970s and revival in the 1980s, string theories have 
occupied a central position in theoretical physics in the quest to unite gravity with 
the quantum field theories associated with the standard model. Our motivation 
within the string context is to implement the non-perturbative approach outlined 



in section (2.3.3) to the simplest string-based quantum field theory - the bosonic 
string including the string-axion - and examine any cosmological predictions the 
solutions suggest. As such, this thesis makes no attempt to review the breadth 
of topics within string theory such as advances in superstrings, D-brane and M- 
theories. Rather we focus on the bosonic string as a conformal field theory and 



its applications in cosmology. In section (3.1.1), scale invariance was discussed 
in the context of QCD where it was noted that the breaking of this symmetry 
results in the particular form of asymptotic freedom demonstrated by the strong 
interactions. Scale invariance is a specific case of wider conformal invariance, 
and this symmetry provides important constraints in forming many viable string 
cosmology configurations. The following account is derived from [77] and [78J , and 
other sources where noted. 
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5.1 Bosonic String Theory 



5.1.1 Introduction to the Bosonic String 

String theory provides what many consider to be the leading candidate for 
incorporating the highly successful quantum field theories of the standard model 
with Einstein's equally successful theory of gravity. [77] notes areas where the 
picture of our universe based on these two theories is incomplete. Firstly, the 
standard model implies a particular pattern of fields and particles with certain 
initial conditions required. These conditions are not naturally derived from the 
standard model's equations but rather need to be set by hand or from observation. 
Secondly, a quantum field theory of gravity is non-renormalisable in four spacetime 
dimensions. Thirdly, classical gravity breaks down at its singularities. With 
regards to the second point, the non-renormalisable nature of quantum gravity 
can be qualitatively seen by dimensional analysis. We consider quantum gravity 
as the exchange of gravitons between two massive particles with associated higher 
order quantum corrections. The one-graviton exchange is equivalent to the classical 



(a) 



(b) 



ZJULSLMJULZ 



(C) 



CJUliUULSLSLil 



Figure 5.1: Simple illustration of quantum gravity: (a) two free particles (zero graviton), 
(b) one-graviton correction and (c) two-graviton correction, [77J 



theory of gravity and its amplitude is proportional to Gn, Newton's constant. The 
ratio of the amplitudes of (b) to (a) must be dimensionless and thus proportional 
to GjyE 2 ^ 1 ^ 5 where E is the energy scale, as this is the only dimensionless 
combination of the parameters available. With H = c = 1 one has by definition 
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Gn oc Mp where Mp is the Planck mass. Thus the one- loop to tree level amplitude 
ratio is proportional to (]^) 2 , meaning it is weak at the low energies of particle 
physics effective field theories but perturbative quantum gravity field theory breaks 
down at high energies > Mp. The two-graviton to tree level ratio similarly diverges, 
proportional to (;J^) 4 - Such divergences have not been shown to be controllable 
via renormalisation as was successfully done with QED and other field theories 
in the 1970s and a non-zero fixed point for the quantum field theory of gravity is 
not known. Alternatively it can be considered that new physics emerges beyond 
the Planck energy scale, and string theory is a leading candidate. To avoid the 
high energy divergences, one dimensional strings (as opposed to point sources in 
conventional particle physics) are considered to be the basic building blocks of 
nature, moving in a two dimensional worldsheet. This one dimensional string lives 
in D dimensional spacetime, called target space. Very qualitatively, high energy 




particle 

Figure 5.2: Particle, String and Brane in Two Space Dimensions 

divergences may be "smeared" out over the string to reduce them. The length of 
a string corresponds to the Plank length lp which is considered fundamental. In 
classical, non-string physics a zero-dimensional point source particle moves through 
spacetime and its motion can be described by D — 1 functions of time t = Xq, (or 
X l (X°) in Minkowski notation). In a covariant form, the motion can be described 
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by D functions of X^(t) where r is a newly introduced worldline parameter. A 
key point is that the choice of r, as an additional, non-physical degree of freedom, 
should not affect the physics. Thus: 

X"V(r))=X"(r), (5.1) 

with the simplest Poincare-invariant action given by (where a dot indicates a 
derivative with respect to r): 



S P = -m j dT\J-X^X^. (5.2) 
This has the two local symmetries of worldline reparameterisation-invariance (i.e. 



X'^{t'{t)) = X M (r)) and Poincare invariance. The action (5.2) may be expressed 
in an equivalent form with the introduction of an additional parameter r\ (t) where 
rj = v/— 7 Tr (r) where 7 TT is the worldline metric. Now the action can be expressed: 



1 



S'p = tt / dT(TT\X»Xp - r]m 2 ), (5.3) 



2 

where m can be considered as the particle's classical mass. While S' P contains an 
additional parameter rj, it is quadratic in the derivatives of X^ and its quantization 



in terms of a path integral can be more easily achieved. The action in (5.3 ), however 
turns out to be non-renormalisable in four dimensions as discussed qualitatively 
in terms of dimensional analysis above. In [79], Polchinski notes that an elegant 
means of removing these divergences in that instead of a point we consider a 
one dimensional object as fundamental, a string, with an additional worldsheet 



parameter a as well as the r in (5.3), as illustrated in Figure (5.3). With a string 
as basic, a two dimensional worldsheet is traced out and one may define the an 
induced metric such that: 

h ab = daX^X^ (5.4) 

where a and b run over the two world sheet coordinates a and r. The parameters 
r, a form a two dimensional vector we will denote as £ a = (r, a). It is termed 
"induced" as it is superimposed on spacetime target space via the metric g^ u and 
may be written: 

h ab = d a X^d b X u g^(X a ). (5.5) 



5.1 Bosonic String Theory 



87 




Thus the induced metric is defined by the spacetime manifold in which it is 
embedded and whose metric may be dependent on the spacetime coordinate X a . 



It is noted that in eq (5.15) the index a is not summed over and thus X a acts as a 



scalar in this equation. The simplest Lorentz invariant action (analogous to (5.2) 
for a worldline based theory), the Nambu-Goto action is given by: 

SNc(h) = J drda^-det(h ab ), (5.6) 

where a' is a parameter of the string theory with units of inverse energy squared 
(or length squared in a H = c = 1 convention) which is inversely proportional to 
the tension of the string, (a' is termed the Regge slope). One may express the 



action in a manner analogous to (5.3) in the Polyakov action: 

Spo = J drdai-^h^daX^X,, (5.7) 

where •y ab is a newly introduced independent worldsheet metric with 7 = det7 a h 
and 7 afe = — 7^. It is introduced in order to make the resulting Polyakov action 
more suited to the path integral formalism, while it introduces an additional degree 
of freedom compared to the Nambu-Goto action in that 7 a b is proportional to h a f, 
with the proportionality constant to be fixed by hand. The Nambu-Goto action 
may be recovered from the Polyakov action by the relationship: 

hab(-h)-3 = 7o6 (-7)-!, (5.8) 
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and the definition of h a b in (5.8). As evident in the definition for h a b in eq (5.15) 



jab also has an indirect dependence on the spacetime coordinate X a . The form of 



(5.7), like (5.3), does not contain a quadratic under a square root and is easier to 
quantize. There are three important symmetries which are manifest in the form 
of Sp in terms of the invariance of the underlying spacetime physics. 

• Poincare invariance of the spacetime coordinates (called target space) in 
which the string is embedded (where /j, = 1....D, in D dimensional space 
time). The worldsheet metric jab remains the same in any Lorentz frames. 
This is a global symmetry of the action Sp a and affects the spacetime 
parameter and worldsheet metric as follows. 

= A^X V (C) + A:". (5.9) 

V«6(0= 7-6(0. (5-10) 

• Diffeomorphism (or diff) invariance of the two worldsheet dimensions 
implies that physics is independent of the choice of worldsheet parameters 
r, a and can be expressed by the following transformations which leave the 
action invariant: 

X'"(0 = *"(0, (5-11) 

while the worldsheet metric transforms as a covariant tensor under a shift in 
coordinates £ a (r, a) — > ^ /a (r', a'): 

7a6(0 Ya*(0 = f^T^O- (5-12) 

Diff invariance implies that the physics on the target space is independent 
of the choice of worldsheet parameters. 

• Weyl invariance of the two worldsheet dimensions is a local rescaling of the 
worldsheet metric 7 a &(0 which leaves the action and physics on the target 
space invariant: 

7a6(0 lUO = exp(2w(0)7«6(0. ( 5 - 13 ) 
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and 

X">(Z) = X"(t). (5.14) 

This symmetry is not present in the Nambu-Goto action, and arises in Sp Q 
as the worldsheet metric jab is determined up to a scale factor represented 



by the proportionality relationship in (5.8) which can be rephrased 



h ab = f(Olab, (5.15) 

with /(£) a non- vanishing worldsheet function. The fact that this propor- 
tionality is governed by f 2 means that the signatures of the metrics h a b and 
7 a t, agree. Choice of /(£) is analogous to choice of gauge in quantum field 
theories. In the Nambu-Goto action, /(£) is a constant and the symmetry 
is not evident. (We note here that Weyl invariance is a generalised form of 



the scale invariance of QCD expressed in (3.11), which is the global version 



of (5.13) where u(r,a) = constant). 



Thus through the diff and Weyl invariance in the worldsheet metric 7^, 



the action (5.7) contains three redundant degrees of freedom. On the other 
hand, the quadratic form of the X M "fields" in this action as compared to 
the Nambu-Goto action makes path-integral quantization possible. These 
redundancies need to be fixed however, in the course of the quantization 
process. 

The metric 7o&(6) is a symmetric two by two matrix meaning it has three degrees of 
freedom. One may, via suitable choice of the diff transformation, fix two of these, 



as noted by the contractions in the transformation (5.12). The remaining degree 
of freedom, associated with the Weyl invariance, may be set by choice of oj(£). So 
one may define the worldsheet metric: 

7a6 (£) = exp(2a;(0)7a6, (5.16) 

where here both 7 afe and are fixed. A common choice of 'jab is 5 a b, the flat 
metric, known as the fiducial gauge. The string parameterised by r and a in Figure 
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(5.3) is an open string. If the string has a defined length, /, one has for an open 
string: 

-oo<r<oo, 0<a<l, (5.17) 

and the worldsheet is an infinite strip of width I. One may also consider a closed 
string which forms an infinite cylinder such that it has the constraints: 

X"(r, 0) = X>(t, 1); lab (r, 0) = lab (r, I), (5.18) 

that is, the parameter a is periodic with period I. So far in our discussion, the 
worldsheet metric •y ab has a Euclidean signature (+,+). However, one may have 
a worldsheet metric with a (— , +) Minkowski signature and the r parameter may 
be considered analogous to the time coordinate, X° in spacetime. In this case, the 
indices a, b may be raised or lowered by 7^ in conventional Minkowski fashion. As 
a final general point [77] notes that all string theories contain closed strings but 
not all contain open strings as the former can be created from the latter. 

As a note on convention, in this thesis, worldsheet indices will always be denoted 
by a, b with spacetime indices left in the conventional Greek nomenclature (/i, v...). 



5.1.2 Classical String Actions and Constraints 

In the Polyakov action for the string, the spacetime coordinate and the 
worldsheet metric •y ab are treated as independent variables of £. In effect, A M 



can be regarded as a scalar field, (with /i an inactive index in (5.7)), coupled to 
the metric 7^. One may obtain two sets of constraints on the classical string via 
minimising X^ and 7 a b with respect to the action. 

SS Po = f^S lab = ^6X^ = 0. (5.19) 
d-f ab dX» 

We consider first the variation of Sp Q with X^, |78j . 

5S Po = J drda5X»d a ((-7)^^) , (5.20) 
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with total derivatives set to zero and left out here. This results in the constraint: 
da ((-i)h a %X^ = (- 7 )2V% = 0, (5.21) 
where V 2 = d a d a . An additional surface term will result if boundary conditions 



such as (5.17) are placed on the string. This surface term vanishes if: 

<W(r,0) = d a X^(r,l) = 0, 



(5.22) 



i.e. the derivatives of with a vanish at the boundary (known as the Neumann 
conditions, [77]). In the closed string the surface term vanishes naturally. Variation 
of 7o6 with respect to the action defines the stress energy tensor for the worldsheet, 

m 



T ab (r, a) 



4tt dS Po 
7=7 d% 



1 

~a' 



ab 



d a x^d b x^ - - 7 ab (a c x^9 c x A , 

2 



Diff invariance means that: 



V a (d a X fl d b X ll ) = V a T ab = 0. 



(5.23) 



(5.24) 



Weyl invariance as governed by the transformation in (5.13) means that an 



infinitesimal variation of the parameter 7 a & is (to first order in the expansion of 
exp 2cu): 

Sj ab = 2w7 ab , (5.25) 



which when substituted into the equation of motion for ^ ab in (5.19) gives: 

dS Po 



2u% 



ab 



fry, 



o. 



ab 



which, since u is arbitrary, leads to: 



dS Po 

lab- = 0, 



#7, 



ab 



(5.26) 



(5.27) 



finally leading to the trace of the stress energy tensor vanishing, a direct 
consequence and condition of Weyl invariance in the classical string. 

dS 



lab 



#7, 



ab 



0^T a a = 0. 



(5.28) 



5.1 Bosonic String Theory 



92 



Thus we have, in solving the equations of motion of the worldsheet metric, 
'jab deployed two key symmetries (diffeomorphism and Weyl invariance) of the 
Polyakov action to derive results. Further, while our discussion relates to bosons 
only, fermions may be considered by the addition of super-symmetry, and this 
forms the basis of the current generation of superstring theories (which will not be 



discussed in this thesis). The Weyl and diff invariance of (5.7) are fundamental to 
the quantization of the classical string and in setting boundary conditions for use 
in string cosmology. 



5.1.3 Conformal Invariance 

Continuous conformal transformations on a manifold rescale the metric but 
preserve locally angles (and hence shapes). Considering D dimensional geometry 
with metric and line element ds 2 = g^ u dx^dx u , under a change of coordinates 
x — > x' one has: 

9^ -> V(x') = dxffldx> „ 9ap(x). (5.29) 

The conformal group is the group of transformations which leave the metric 
invariant up to a position dependent scale factor such that g^ u — > g'^ u (x r ) = 
u(x)g^ u , [80J. These are the transformations which preserve the angle between 
two vectors on the manifold. The infinitesimal generators of the conformal group 
can be found by considering the transformation x M — > x^ + e M , where e M is small. 



To satisfy (5.29) one must have (quoted here without derivation), 



2 

<V„ + dytp = —{d ■ ejr)^. (5.30) 



The following types of transformations leave the metric invariant as in (5.29 ) (where 
a M and 6 M are constant vectors): 



translations independent of x: e 



rotations: e M = a£x u ; 



scale transformations: e M = \x^ 
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• special conformal transformations: e M = b^x 2 — 2x tl b ■ x. 



The first two are a representation of the Poincare group (u> = 1), the third are 
scale transformations (u = A~ 2 ) while the special conformal group has u(x) = 



(1 + 2b ■ x + b x ) . When dimensionality D = 2, (5.30) results in : 



<9iei = d 2 e 2 \ d x e 2 = -d 2 £i, (5.31) 

which are the Cauchy-Riemann equations providing a sufficient condition for a 
differentiable function to be holomorphic. Thus we may write the above coordi- 
nates in the complex plane: z, z = x l ± ix 2 and the two dimensional conformal 
transformations are a representation of analytic coordinate transformations in z, z: 

z _►/(*), (5.32) 

I ft f I 2 

and ds = dzdz — > ^ rfzrfz with w = | ^ | . The local symmetry group is infinite 
dimensional in the special case of D = 2 which provides for the complexity of a 
two dimensional conformal field theory such as the bosonic string theory under 
discussion here. (We note that D here refers to the dimensionality of the geometry 
in general and not necessarily physical spacetime dimensionality). [HQ] notes that 
a combination of diff and Weyl invariance in a two dimensional field theory with 



a Euclidean metric as in Section (5.1.2) is equivalent to conformal invariance on 



a complex worldsheet such as the analysis leading to (5.31). We may restate the 



Polyakov action by expressing the worldsheet coordinates (r, a) in the complex 



plane. First we consider the Polyakov action (5.7) with a flat Euclidean worldsheet 



metric 'jab = 5 a b = (1,1), [77j . and hence the change of sign compared with eq. 



(|5_7|): 

S = -i- / d'ridxX^dtX^ + d 2 X»d 2 X fl ). (5.33) 



Arra' 

The choice of -j a b = 5 a b fixes the local gauge of the worldsheet metric to be 

flat, a commonly used configuration for quantization. The complex worldsheet 
coordinates are: 

z = a 1 + ia 2 , z = a 1 - ia 2 , (5.34) 
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and 



d z 



1 



(<9i - id 2 ), d- z = -(dx + id 2 ) 



with terminology d = d z and d = d z used. Now (5.33) becomes: 

1 



5 



27ra' 



d'zidX^dXf, 



The equation of motion obtained by minimising the field becomes: 

ddX»(z,z) = 0. 



(5.35) 



(5.36) 



(5.37) 



In (5.28) Weyl invariance was shown to result in a traceless stress-energy tensor. 



In the the coordinates (5.34) this becomes: 

T - = 

± zz — 



0. 



while the equivalent of (5.21) becomes: 

Bt zz = 



dT- z 



0. 



(5.38) 



(5.39) 



We denote T zz = T(z) and = T(z), and for a free massless scalar theory the 
analogy of ( 5.20[ ) leads to the following forms for the stress-energy tensor. 



T 

A zz 



T 



T(z) = --dX^dX"^ 
a 

f(z) = --BX^dX^r,^ 



a 1 







(5.40) 



which are equivalent to the requirements of diff invariance that the stress energy 



is conserved (5.24) and from Weyl invariance that it is traceless (5.28). Invariance 
of the target space physics with respect to conformal transformations of the 
worldsheet metric is a highly desirable feature of a string theory attempting 
to describe the physical universe. We have stated that a test for conformal 
invariance in both the real and complex planes of the worldsheet space relate 
to the conservation of and trace-vanishing properties of the stress energy tensor. 
Discussion thus far has remained at the classical level. We now consider 
implications of quantization on the above results, and the constraints imposed 
on this process by the conformal conditions. 
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5.1.4 Quantization of the Bosonic String 

Rather than offer a comprehensive account of the derivations, methods and tools 
involved in quantization of the bosonic string, common results are simply quoted. 
Following the above account of how conformal invariance presents key constraints 
in the classical theory, our aim is to build a brief background leading to an 
explanation of Weyl anomalies in the quantized theory. In a quantum field theory 



of the classical bosonic string, the fields are promoted to operators. In (5.36) both 
and 7 become operators and are functions of £. The theory may be quantized 
using the path integral formalism, where the partition function of an action S[X, 7] 
is: 



Z= J[dXdj] exp(-S[X,j}), (5.41) 



where this integral runs over all worldsheet metrics 7 and all values of X M (£) 
(where the index /1 is not integrated over). The integration over the world sheet 
g?7 is conducted over the three degrees of freedom in the metric 7 ab , and it is 
here that anomalies arise which set a limit on the spacetime dimensionality of 
the effective theory if conformal invariance is to be preserved. In the quantized 
theory, T a a does not necessarily vanish as it does in the classical theory, but is 
proportional to the central charge of the theory, c. The value of c is found through 
the summation of the central charges arising from various sectors of the theory. 
For the bosonic string, the central charge has contributions from the spacetime (of 
dimension D) coordinate acting as a scalar field (cx = D) and from the "ghost" 
sector (c g = —26) so that (results stated rather than derived here): 

c = c x + c g = D - 26. (5.42) 

The setting of c = provides for the well known D = 26 constraint of the free 



bosonic string. In quantizing the Polyakov action (5.7), we use a Euclidean metric 
7 afc of signature (+, +) and with dadr replaced with d 2 £. 



(5.43) 
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We generalise (5.43) by adding a term which is allowed by the Poincare and diff 
invariance. This results in a term (Sr) which may be added to Sx'- 



S 



R 



1 

47T 



(5.44) 



Here is the Ricci scalar constructed from the metric 7 ab and A is a coupling 



(which we will later relate to the dilaton). The term (5.44) is invariant under 
Poincare, diff and Weyl transformations on a worldsheet without boundaries, [75] , 
resulting in our most general classical action, (which we denote simply as S here): 

S = S x + S R = J d 2 ^ {r'daX^X, + Xa'R^) . (5.45) 

The expression ^- J d 2 ^^ R^ is known as the Euler number, x, and depends only 
on the topology of the worldsheet. Under a path integral operation as in (5.41 ) this 
term produces a factor of exp(— A%) which impacts only the relative weightings 
the various topologies of the worldsheet hold in the path integral summation. 
Interactions on the worldsheet, including string self interactions are determined 
by the topology as defined by 7 a6 , which in turn has a unique Euler number. 



This Euler term in the path integral quantization of the action (5.45) controls the 
coupling constants in string theory, with the couplings being related to exp(A). 



The path integral (5.41) contains redundancies resulting from configurations of 
(X, 7) and (X', 7') related by the diff and Weyl transformations describing identical 
spacetime physics and it needs to be gauge fixed to avoid unphysical divergences. 
This can be done by the Faddeev-Popov method, [ST], used commonly in gauge 
theory, where slices of each identical gauge are taken only once. A common gauge 
choice is to fix the metric such that 7^ = 5 a b — (1, !)• 



In quantum field theories, Feynman rules for calculating cross sections denote 
a particle at external point x with momentum p by exp(— zp.x). Likewise, string 
interactions are usefully described with the aid of vertex operators. In conformal 
field theory, ([77], p 63) there is an isomorphism between the set of allowed states 
in the quantised theory and the set of local operators describing the quantum field 
theory. While details of this concept will not be developed here, as an example, 
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the lowest permissible energy state of the closed string spectrum can be described 
by a vertex operator in the string action: 



V = 2g c J <p£Syexp(ik.X), 



(5.46) 



where X = X^(£) and k is the momentum in worldsheet space, and g c is the 
string coupling. In this particular example, the state is characterised by having 
negative mass squared and is known as the tachyon. In this language, the lowest 
states of the closed string act as local operators and can be regarded as interacting 
background fields all of which depend on the X M (£) parameter (where the \x index 
remains inactive thus far). The free bosonic string theory represented by the most 



general version of the Polyakov action (5.45) thus becomes an interacting theory 



with the vertex operators representing the lowest quantum states playing the role 
of interacting fields. In the next section we will state the full action we will use 
for string cosmology, based on the expansion of the Polyakov action to include low 
level energy states of the string represented by local vertex operators. 



5.1.5 The Bosonic String in Curved Spacetime 

We now replace daX^d^X^ with g^ u (X)d a X fl dbX u and thus are considering the 
string in a curved target spacetime. (In the notation g^ v (X) we leave out the \i 
index on X as it is still inactive). One may consider the curved spacetime metric 
g fMU (X) = r]^ u + Xfii/(X), where Xnv{X) is an incremental variation in the metric 
at point X. The term Xnv{X) can be equated to a vertex operator representing a 
coherent vibrational mode of strings known as the graviton, whose dynamics are 
denoted by s^ v and where g c is a coupling: 

X^iX) = -4ng c exp(ik.X)s^ u . (5.47) 

Expanding this concept to include the backgrounds of other massless open or closed 
string states we state the following bare action, denoted again here simply as S. 
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+a'R (2) <S>(X) + 4WT(X) j , (5.48) 
Within this the following vertex operators represent string states acting as coherent 



fields (the tachyon vertex operator T{X) was defined in (5.46)). 



g^iX) = T}^ - 4:7ig c exp(ik.X)s^ u , (5.49) 
B^ U (X) = 4:7ig c exp(ik.X)a^ u , (5.50) 
$(X) = Aixg c exv(ik.X)<p. (5.51) 



We should make several comments at this stage. 



It was noted in (5.45) that an additional term in the free string action is 
allowed by the Poincare and diff invariance of the Polyakov action. This term 
is described by the Euler number j- J (P^y/j and a coupling constant 
A. A basic premise of string theory is that there are no predetermined 
parameters which need to be set by hand. The coupling A in curved spacetime 
is thus promoted to a dynamical parameter, known as the dilaton, a scalar 
field notated as <&(X), which appears naturally in the spectrum of the 
massless closed string and plays a key role in the couplings associated with 
the background fields of the quantized string. 



In (5.48) there are two fields coupling to the derivative term d a X fl dbX u . The 
term / ~f ab g flu (X) is present in the classical action when the spacetime metric is 
generalised to curved spacetime and this term represents the graviton string 
state. 

The term containing i e ab B^ v {X) is known as the Kalb-Ramond field and is an 
antisymmetric tensor field (that is, antisymmetric in two spacetime indices 
fi and v). The fact that this term is complex arises from the Euclidean 
continuation. It is analogous (in a tensor form) to the vector gauge field 
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of QED which couples to the fermion fields in the Yukawa fashion, that is: 



notes that the quantised state at fixed momentum of the closed string 
can be represented by an arbitrary matrix Mu of size D — 2 (where D is the 
spacetime dimension) multiplied by basis vectors representing the massless 
string states. Such a matrix Mu may further be generalised as the sum of 
symmetric traceless part, a antisymmetric part and a scalar trace - giving rise 
to the three background fields noted above - the graviton, the antisymmetric 
tensor and the dilaton. 



The action in (5.48) is the starting point for our treatment of string cosmology, 



with the caveat that we do not consider the tachyon field in our study. 



5.1.6 The Weyl Anomaly 



Equation (5.28) states that the a condition of conformal invariance in the classical 
bosonic string is a traceless stress energy tensor, T a a = 0. In our discussion of QCD 
it was noted that when renormalised, scale invariance (the global version of Weyl 



invariance, as shown by (3.11 )) is broken resulting in an effective theory dependent 



on energy scale. In standard model QFTs this is not a problem as one is concerned 
with the infrared scale, and one may compute beta functions such as /3(A) = 
(where \i is energy scale and A is the coupling) to monitor the variation with energy 
scale of the theory's parameters. However in the quantization of the bosonic string, 
maintaining Weyl invariance is usually conducted so as to derive invariant target 
space results. (A class of string theories (for the bosonic string) where the D = 26 
condition, is by-passed are known as non-critical string theories). In the quantised 
bosonic string, T a a = vanishes only under certain conditions and these form the 
basis of constraints for a critical, effective string theory actions. Classically, the 
stress energy tensor is derived from the variation of the action with the metric, 



(5.23). In the quantum field theory approach, the stress energy tensor is related 
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to the infinitesimal variation of the partition function with the metric, [77| : 

An 5 



rpab 



(5.52) 



ab 



where Z is computed as in (5.41). In the quantised theory, T a a must be computed 
and may not vanish. [79] shows that in a quantised theory, one has for a free string: 



(5.53) 



where c is a constant proportional to the central conserved charges of the fields in Z 
and is the Ricci scalar associated with the metric y a b- In the free bosonic string 



theory represented by (5.45) the central charge from the ghost fields arising from 



the Faddeev-Popov methods to fix the gauge is —26. The charge arising from the 
X^ L fields is simply D, the spacetime dimensionality. This leads to the commonly 
quoted result that free critical bosonic strings require 26 spacetime dimensions to 



display Weyl invariance. The string action (5.48) contains additional fields which 



may break Weyl invariance. We quote without proof the result, 

S 



8g* 



<S> 



(5.54) 



where S is the action in (5.45) and (3 l are the beta functions for each coupling i 



describing the flow of the couplings with the Weyl rescaling parameter u(a). The 
notation < ( ) > is defined by: 



< 



•) >=- 



1 

Z 



[dXdy]{ )exp(-5[X, 7 ]), 



(5.55) 



where Z is defined in eq. (5.41 ). Here g % is the set of couplings associated with the 



background fields g^ u , $, in the action (5.48). As in the case of the free string 



Weyl invariance occurs when the (3 l expressions vanish. The beta functions are 
perturbative expansions in the a', the Regge slope and we quote here from [77] : 



1 



1 



:/"' = -2^7 ab < 9 a XIX 9 b X" > -^(3^ ab < d a X»d b X» > --^R^. (5.56) 
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To zeroth order in a', the beta functions are as follows 
P* = - l -V K H KIXV + d K <j>H KIXV + 0{a!) 

P* = ^^-\v 2 <P + dyd,<p-^H, up H^ + 0(a'), (5.57) 

where R pu is the target space Ricci tensor. H^ K \ is the field strength of the 
antisymmetric tensor such that H^x = d^B K \ + d K B\ p + d\B^ K . From the form 
of it can now be seen that with background fields the space time dimension D 
must no longer necessarily be 26 to preserve conformal invariance as terms of a' 
or higher order could potentially be made to cancel the g ~ 26 term. In addition, 
the diff symmetry on the world sheet means that the beta functions are invariant 
under the following transformation, which is a field rescaling: 



pi^fii = pi + S(3\ (5.58) 



where we defined 5(3 % in some detail in (6.49). These modified f3 l functions are 
known as the Weyl anomaly coefficients and are often used in place of the l 
functions. 



In setting the /3 or (3 functions to zero (conditions of conformal invariance) one is 
effectively writing down equations of motion which set constraints on the field 
variables, and from this one can determine an effective action in space time 
expressed in terms of the dilaton field, H^x and the spacetime Ricci tensor (all of 
which are functions of X^), [9"E] . 

Seff - ~ J d D X^e~ 2 * (r + 4(0 M *0"<&) - ^Hx»„H^ - + <?(«')) , 

(5.59) 

where k is the gravitational constant in D space time dimensions, and here R is 
the spacetime Ricci tensor (note the target space Ricci tensor has been denoted by 
i?( 2 )). We note at this point the minus sign in front of the dilaton kinetic term in 
this frame (know as the u-model frame). The term "frame" in this context refers 
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to the spacetime metric configuration the action exists in. The positive sign is 
restored for the dilaton kinetic term in the physically relevant Einstein frame, by 
making an allowed transformation of the metric. A common practice is that the 
spacetime metric g^ v is transformed in a spacetime Weyl transformation such that: 

9^ -> 9%, = e'^g^u- (5.60) 



Under this transformation, the curvature scalar R in eq. (5.59) will have the 
standard co-efficient (under general relativity) - i.e. ^r. In the Einstein frame the 
effective action is, 196 1: 



-eA*!l£^) +OM l. (5.61) 

In the Einstein frame the transformed parameters superscripted with an E. In this 
frame, the spacetime curvature scalar is denoted R E and is canonically normalised 
with the coefficient 7^2 so as to be relevant to physical observations. Further, one 
may see that with D = 4 the dilaton kinetic term changes sign and is now positive. 

With this brief and largely qualitative introduction to string theory, we are now 
ready to describe the basics of string cosmology. 



5.2 String Cosmology 



5.2.1 Introduction 



The following is based on [90], [93], [95], [96] and sources where noted. The study of 
the early universe is highly concerned with the study of the gravitational field and 
its evolution. If one is to take string theory as representing the fundamental fabric 
of our universe and uniting gravity with quantum field theory, it would appear 



logical to take a bare action such as in (5.48) as a starting point in studying the 
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evolution of fields and particles in the universe. Likewise, cosmological observations 
are also one of the few feasible testing grounds for string theory, and currently show 
more promise than earth-bound collider experiments in testing its predictions. 
Conventional cosmology is based around the A-CDM model incorporating the big 
bang, an inflationary scenario and cold dark matter/dark energy both of unknown 
origin at present. Most string cosmology models seek to incorporate the dynamics 
of the A-CDM model and shed light on some of the unexplained phenomena such 
as dark energy. Moreover, the last twenty years has seen advances in string 
theory and allow it to more easily make contact with real-world observations. 
One such advance is moduli stabilisation, [93]. Conformal invariance dictates that 
critical string models may need to live in more than four spacetime dimensions. 
String models deploy gravitationally coupled scalar fields, known as moduli, to 
compactify the additional dimensions to a level where they cannot be detected by 
current observations. However, the moduli fields should also present cosmological 
observables and moduli stabilisation is the term given to theories which counter the 
effects of moduli to non-detectable levels. A further challenge of string cosmology 
is in isolating its predictions compared to those feasibly arising from conventional 
QFTs. Given the degenerate nature of string solutions one tactic, [93], is to limit 
string cosmology theories to those which produce outcomes which could not arise 
from conventional quantum field theories. One such area (which we will not discuss 
in this thesis) are string-based topological defects, known as cosmic strings. Prior 
to a description of our research, we will outline a variety of string cosmology 
approaches. 



5.2.2 String Cosmology Approaches 

There are many string cosmology models which address a range of scenarios such 
as inflation, the graceful exit and pre-Big Bang (see for example jHU] )• We do 
not attempt a comprehensive review in this thesis but rather introduce an early 
string cosmology model which can describe a static and also expanding universe 
in four dimensions. In the previous section, we described the path taken from a 
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bare worldsheet action such as the Polyakov action in (5.36) representing a free 



string to an action representing the bosonic string in curved spacetime though to 
the effective action in D spacetime dimensions useful for physical observations, 



(5.61). The derivation and evaluation of such an effective action motivated by 



early universe phenomena forms the basis of string cosmology. As the observed 
universe is isotropic and homogeneous, one is generally interested in configurations 
where the spatial co-ordinates are fixed with the time component varying. The 



conformal invariance conditions (i.e. the terms in (5.57) vanishing) determine the 
equations of motion for the background fields in target space. A common starting 



point is the the bare action (5.48) with the background fields displaying variation 

only in the time (A ) direction with fi, v — 0, D — 1. We exclude for the time 

being the tachyon field, (but make some comments on this in the next section). 

1 



S n 



<:t 



M 



cry^[(r b 9,u(X°) + ie ab B^ lu ( A ) ) d a X^d b X v + a'R^{X )}. 



(5.62) 

Here we again note that is the worldsheet Ricci scalar, and it is noted that the 
dilaton field $(A°) is one order of a 1 higher than the g^ u (X°) and -B Mi ,(A°) fields. 
We present without derivation and as an illustrative example one configuration 
used in string cosmology, the so-called linear dilaton background which satisfies the 
conformal invariance conditions and also demonstrates unitarity of the spectrum, 



(5.63) 



g^(X°) = r^, 5^(A°) = 0, $(A°) = -2QA , 

where Q is a constant. Q is often referred to as the central charge deficit, [97] 
and it describes the evolution of the dilaton in A . The stress energy tensor in 



the notation of (5.40) where z is the complex worldsheet coordinate, in the linear 



dilaton configuration of (5.63) is given by: 

1 



T 



^A^A<V(A U ) + Q^A U . 



(5.64) 



The Weyl anomaly (central charge) associated with this tensor is: 



c = D- 12Q 2 = 26. 



(5.65) 
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Thus the spacetime dimensionality to satisfy vanishing Weyl anomalies is now also 
dependent on Q (unlike the free string case where D = 26 was a condition to cancel 
the Weyl anomaly) which represents the dilaton field strength. String cosmology 
models where the spacetime dimensionality need not necessarily be D = 26 (for 
the bosonic string) due to the presence of a non-zero Q value are often termed 



non-critical string theories, [98]. In eq (5.65), it is noted, [98], that the spacetime 



dimensionality D may be expressed as D = d + cj where d is the number of 
large physically apparent dimensions and cj is the internal charge related to the 
compactified dimensions; further, cj and 12Q 2 need not be integers. In the linear 
dilaton background the line element is, [96] : 



ds 2 E = e^rj^dX^dX", (5.66) 

which may be expressed in the form of the Robertson- Walker metric, where t is 
cosmic time: 

ds\ = -(dt) 2 + tfdX'dX^ij, (5.67) 
where the scale factor, a(t) = t. 

D - 2 _____ ( 2Q „„ 



Since a(t) = d ^ — this describes a non accelerating, flat universe in which the 
dilaton field in terms of t is given by: 

«&(*) = (2 -D) In (5.69) 

This model may be generalised to one incorporating a curved spacetime background 
and a non zero value of B^X ). One can set the spacetime dimensionality in 
terms of observable dimensions to four and assume that any remaining dimensions 



as required by the constraint in (5.65) are compactified. The antisymmetric field 
strength may be expressed in terms of a scalar field h(x) (sometimes termed the 
pseudo-scalar axion field in this context, differing from the QCD-axion) 



H mv = exp(2$)e plluX d x h. (5.70) 
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The equations of motion associated with the Weyl anomaly in equations in (5.57) 
are as follows: 

1~ 1 XO 1 



R 





8c 



2 p 2 
V 2 h + 2d u kd»$ 



12Q 2 



-3e~ 



-R 



V 2 $ + -(<9$) 2 



,2$ 



(5.71) 



The effective action (termed here simply as S e ff) obtained when the (5.71) are 
satisfied, in these conditions is: 



S, 



eff 



d Xy g E 



R L 



,2$ 



(dhf 



(5.72) 



2 V p ' 2 v ' 3 
where this is in four non-compactified spacetime dimensions and where 8c is defined 



in the third line of eq (5.71). This admits as a solution the line element in polar 
coordinates in the Einstein frame: 

-dt 2 + a{tf 



ds\ 



dr 2 



+ r 2 {d9 2 + sin' 



(5.73) 



1 — nr 2 

where a(t) is the scale factor and k is a parameter describing spatial curvature. 



Eq. (5.73) and the solutions to (5.71) can be expressed in terms of the Hubble 



parameter H(t) = ^1, 
quoted here from [96]): 



which describes the evolution of the universe (equation 



H + QHH 



4/c 
a(t)2 



H 



H + 3H 2 + 



2k 
a(t)2 



-AH + 



4k 



{5c) 2 h 



2U2 




a(t) 2 36a(t) 6 (H + 3H 2 + 



2k 
a(*) 2 



(5.74) 



There are two types of solutions to this equation, both of which have k as non- 
negative and therefore describe closed de Sitter universes. One is the static Einstein 
universe and the other is a linearly expanding universe, with a(t) = t, as in 
the linear dilaton case. The above description of the linear dilaton non-critical 
string cosmology model is illustrative of an early model of an expanding and 
homogeneous universe in four dimensions. We note our model differs in that we 
derive a logarithmic relationship between the dilaton and X°, which nonetheless 
can describe an isotropic and linearly expanding universe by choice of our free 



parameters as described in section (6.3.3). 
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5.2.3 Other String Cosmology Approaches 
5.2.3.1 Liouville String Theory 

String theory is a two dimensional field theory where a parameter of the theory is 



a length scale / as in (5.18). Theories which are independent of the length scale 
are termed critical string theories and are Weyl or scale invariant, [77]. (Note: the 
term "non-critical" string theory is generally used in the context where D does 
not have to be limited to 26, such as Liouville models described here where the 
conformal invariance condition is relaxed. However we note that our model is not 
limited to D = 26 but we demonstrate conformal invariance). When the theory is 
dependent on length scale, the Polyakov classical Euclidean action becomes, with 



comparison to (5.43), and where the length scale I = [l: 

SnC= aL^I ((V*^ftX M ) + a*) • (5.75) 

Diff invariance remains in the theory but now the metric 7ab(£) may now only have 
two of its three degrees of freedom fixed, unlike in our gauge fixing exercise as in 



(5.16) where the fiducial metric was fully defined. Now we have: 

7ab (£) = exp(20(£))7 a6 (£), (5.76) 

where 0(£) is the scale factor of the Weyl transformation of the worldsheet metric. 
In quantizing the critical string, the gauge was fixed and the Faddeev-Popov 
procedure applied to generate the action to be quantized in the path integral 



approach as in (5.41). In applying this to (5.75) one has [77J, section 9.9: 



S= aL ; J d2 ^ l)l + ^exp(20) + ^-{^d a( j>d h( j> + R<f>) 

(5.77) 

While details are beyond the scope of this thesis, the scale factor now has a kinetic 
term and is promoted to a dynamical field in the worldsheet parameter, </>(£), 
known as the Liouville field (not to be confused with the dilaton field also denoted 
by in this thesis). In non-critical strings, this is kept as an extra degree of 
freedom to be integrated over in path integral quantization. This idea was initially 
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proposed in the late 1980s, [HH], jlOOj - [77] notes that the target space physics of 



(5.77) remains the same as compared to critical string theory where the metric 
is fully fixed which can be intuitively seen by the fact that the gauge fixing of 



(5.16) is arbitrary. Non-critical string models for free strings are not limited by 
the D = 26 constraint and the quantized model allows for alternative means of 
calculating string interaction amplitudes. In |101j it is proposed that the Liouville 
field is a dynamical parameter which acts as a renormalisation group scale in the 
quantized string theory flowing from an infrared fixed point to an ultraviolet one 
which represents a critical string vacuum state. The field is further identified with 
physical time. Liouville string cosmology using the non-perturbative quantum field 
theory methods employed in this thesis has been researched and demonstrated in 



5.2.3.2 Tachyon Cosmology 



Here we comment briefly on the tachyon field T(X) referred to in eq. (5.48) and 
the surrounding text. While the graviton, antisymmetric tensor and dilaton fields 
represent massless modes in the lowest state of the most general quantized bosonic 
string spectrum, the tachyon is also in the lowest state but is not massless. While 
we do not derive the result here, [77] notes that its mass squared is proportional 



to 2 — D, and thus in D > 2, is negative. In section (3.1.2) on instantons, it 



was described how negative mass squared states in a two dimensional scalar field 



theory (which is analogous to the tachyon in (5.48)) result in divergent fluctuations 
around the minimum of the field. Many, if not most, string cosmology models rely 
on string actions where the tachyon fields have been removed by adding constraints 
on the closed string, [78] to avoid such instabilities. Tachyon- containing string 
models with their instabilities controlled (by considering that T(X) = is not 
the true vacuum of the theory, [84J) have been deployed recently in early universe 
evolution modeling. We utilise a closed bosonic string action and assume that it 
has been derived so as the tachyon does not arise. In the cosmological discussion 



of our results in section (6.2.1 ) we note that our configuration of the bosonic string 
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without the tachyon cannot give rise to an inflationary scenario as our solutions 
lead to a power law dependence of the scale factors with cosmic time t rather than 
the exponential relationship a(t) ~ exp(Ht) required for inflation. This arises 
from the dilaton factor in the relationship between the string and Einstein metric 



expressed in eq. (5.60) with which one cannot obtain an inflationary scenario 
unless the value of the dilaton field is zero throughout. We note here that the 
non-perturbative methods employed in this thesis and using a tachyon- containing 
model is shown in [82] . In this, the authors demonstrate that a bosonic string 
model incorporating graviton, dilaton, (0), and tachyon, (T), background fields 
can be shown to be conformally invariant in four spacetime dimensions. Further 
this is shown, with the cancellation condition 2(f) + T = 0, to lead to an Einstein 
frame action in the canonical form with no dilaton nor tachyon pre-factors and 
hence leading to a scale factor of the form a(t) = aoexp(Ht) where H is the 
Hubble parameter. The authors note a merit of this model as being feasible in 
four spacetime dimensions and thus not requiring complications associated with 
compactif cation of additional dimensions during the inflationary phase. While this 
thesis does not consider a supersymmetric extension of our model, it is noted that 
tachyon modes would necessarily be eliminated in such a model and we do not 
consider such modes in this thesis. 



5.2.4 String Cosmology and Anisotropy 

In this thesis we are interested in any fundamental spatial anisotropy, or preferred 
direction, which may have originated in the evolution of the early universe as hinted 
by our string model. Analysis of any such anisotropics are a key motivation behind 
the study of the cosmic microwave background from the Wilkinson Microwave 
Anisotropy Probe, WMAP, and other data sources. The main phenomenological 
motivation of our string solution involves coupling it to an electromagnetic field 
and considering the rotational effect on the electric field. The phenomenon known 
as cosmic birefringence, |112j . |115] refers to additional rotation in the plane 
of polarisation by a electromagnetic wave traveling through the universe, after 
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Faraday rotation (caused by galactic magnetic fields) has been accounted for. An 
optical anisotropic or birefringent substance refers to one in which the refractive 
index varies with spatial direction, |119j . In this context "substance" refers to 
the intergalactic distances over which electromagnetic radiation travels prior to 
observation on earth. Such an affect will cause a phase difference or rotation rj of 
the electric field of the radiation along the optic axis. Such an effect is typically 
studied by viewing extra-galactic radio waves emitted by pulsars or other distant 
objects. This may be evidence of cosmological anisotropy on a vast scale, |113] 
and the rotation 9 of a polarised wave of wavelength A is defined by: 



where the a term is the Faraday rotational effect and rj is the wavelength- 
independent cosmic birefringence effect of interest to this study. Polarised 
synchrotron radiation was predicted and observed from distant extra-galactic 
sources from the 1950s and datasets emerged, |113j with information on the 
redshift, z and the angular orientation of the source-galaxy optic axis (which 
sets the source polarisation angle, tp). The fitting parameter a depends on the 
density of electrons and magnetic field lines in the line of sight of the radiation. 
|114] notes that studies to determine cosmic birefringence originating in large 
scale anisotropy account for Faraday rotation on a galaxy by galaxy basis in 
the line of sight of observation, and the residual angle r] is deemed to be due 
to cosmic birefringence. In |113] (1997) a systematic rj angle effect is reported, 
correlated with angular positions and source distances, with dependence on redshift 
z ruling out local effects. More recent work, |117] (2009) notes that several studies 
of the BOOMERanG 2003 flights data (which are devoted to measurement of 
polorisation of the CMB radiation) may show cosmological birefringence effects 
of T) ~ —4.3° ± 4.1°. Meanwhile |118] (2010) reports on further tests to look at 
cosmic birefringence from eight radio galaxies at redshift z > 2, and reports no 
rj rotation within a limit of rj ~ —0.8° ± 2.2°. While experimental research into 
cosmic birefringence remains ongoing, we consider theories which predict it with a 
focus on string related theories. 



0(A) = a\ 2 + rj, 
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[118] notes that the potential for the universe on a large scale to rotate (other 
than through Faraday rotation) the plane of polarisation of light could have a 
number of possible origins: the presence of a pseudo scalar condensate, neutrino 
number asymmetry and the violation of certain fundamental symmetries such as 
CPT. |118] notes that the starting point for modeling cosmic birefringence due to 
pseudo scalar effects is the interaction between the photon and a scalar field (not 
the dilaton here). 

C = -\f, u F^ - \d^d^ - U{4>) - ^F^, (5.79) 

with F fJiV F llu defined in a similar manner to the axion interaction Lagrangian in 



eq (3.57). (One notes the negative sign in front of the pseudo-scalar kinetic term 
^d^4>d^(j) is a matter of convention deployed in this paper, and also by Peccei in [35] . 
We deploy the opposite convention in our treatment of the pseudo-scalar axion in 
( 4.2[ )). |118] notes that in the lowest order in fluctuations, the photon is coupled 



to the time derivative of the scalar field <ft and thus different configurations of i](r), 
where r is conformal time, will lead to differing results for the cosmic birefringence 
angle. The angle r\ is related to variation of 0(r): 

^) = f [0(to)-0(t)]. (5.80) 

|112] and |115j consider the cosmic birefringence effects of the the antisymmetric 
tensor field B^ v in low energy effective string theories. In |115j . the field strength 
of the antisymmetric tensor is coupled to an electromagnetic gauge field to 
produce an "auxiliary" field T^ vX : 

T fluX = VGd^B uX] + ^VGA^F uXh (5.81) 

where F^ u is the electromagnetic tensor and G is Newton's constant. In this it is 
assumed that H^ uX = d[ fl B uX ] varies only with cosmic time, r, and H^ uX = e^y^d^h 
where h is the string axion pseudo scalar field. With a flat spacetime cosmology 
configuration with h = h'r + h (where h! and h are values of h at times r 
and r = respectively) the authors compute a cosmic birefringence angle i] is 
proportional to h'r, that is the cosmic time that has elapsed between the values 
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h! and h . In another scenario associated with large values of r, |115] finds the 
rotation for the radiation and matter dominated time periods goes as tan -1 - and 
tan -1 -3 respectively. In a follow up paper, [1 12j using similar methods but with 
a heterotic string theory, similar results are noted for the flat spacetime case. Our 
treatment of the bosonic string axion similarly couples to an electromagnetic field 
and we likewise consider optical anisotropy effects. 



Chapter 6 

Non-perturbative String 
Cosmology 



6.1 Our Approach to the Bosonic String 



6.1.1 Introduction 



In our treatment of the quantization of the axion, we used an alternative approach 



to the Wilsonian exact renormalisation group (ERG) as outlined in section (2.3.3) 
and implemented in section (4.2). Our aim in this section of the thesis is to apply 
these same techniques to a conformally invariant bosonic string model including 
the string-axion, motivated by cosmology. As with our treatment of the the QCD 
axion, the rationale for applying ERG-based methods resulting in non-perturbative 
equations is that we are examining theories which approach the high energy limit 
in terms of using perturbative techniques. We commence with a worldsheet action 



based on (5.48), without tachyon fields. As with our treatment of the axion, 



where the scale factor / was used to construct our final exact equation (4.14), in 
the treatment of the bosonic string, we utilise A which is a parameter of string 
theory, (where A = in the full quantum theory, and a' is the Regge slope). 



6.1 Our Approach to the Bosonic String 



114 



We seek a A-independent solution to an evolution equation of the effective bosonic 
string action in four spacetime dimensions. We demonstrate that our resulting 
configuration is conformally invariant by showing the beta functions can vanish to 
all orders in a'. This contrasts with most conventional critical string cosmology 
models which are generally limited to the beta functions vanishing at two loops or 
less. We achieve this via a two step process. Firstly, we show the beta functions 
of our solution are self-homogeneous to all orders in the time (X°) dimension. 
Secondly, utilising established methods which state the beta functions may be 
transformed, with the physics unaffected, by arbitrary linear transformations of 
the spacetime coordinates of the string background fields, we show that they 
can be made to vanish. It is stressed that both of these steps are required 
for our demonstration of conformal invariance to all orders. We then look at a 
particular area of cosmological implications. This involves the coupling of the field 
strength H plxv to an electromagnetic field, solving the resulting equations of motion, 
and examining any resulting optical activity demonstrated by the configuration, 
following work done in |115j . We initially fix one index of H pfiu (an anisotropic 
configuration of the antisymmetric tensor) and then we consider the general case 
with all indices free (isotropic). It is stressed at this point that the key result 
of our research in this area is to demonstrate a technique for showing conformal 
invariance to all orders in a' for a bosonic string cosmology model. The predictive 
powers of the model may be incomplete as the resulting metric does not support 
an inflationary scenario. As a note on content, it is acknowledged that a large 
portion of the derivations in this section were published in [106] and |107] both of 
which built on earlier work by Alexandre and Mavromatos in |103] and [104] . 



6.1.2 The Effective String Action 



We commence with the following bare action of the bosonic string in background 

l 



a pv and the dilaton field 



fields, based on (5.48), with A 
varying in the time X° dimension only (i.e. an isotropic and locally flat universe) 



We consider in Euclidean format and hence the imaginary component of (5.48) 
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reverts to real: 

Sx = ~Lj d ^ { A + d « x * dbXV + R {2) ®( x °)}- (6-1) 

This is the bare theory. Here are the microscopic fields over which the path 



integral in (6.2) is taken to compute the partition function. We consider flat 
spacetime rj^ and hence we assume no dependence by the worldsheet metric r j ab 



on X M . We will use techniques outlined in section (2.1) to construct a quantised 



effective theory and define the Legendre effective action V. We note that these 
techniques may only be employed if the evolution parameter is multiplied by terms 
quadratic in the field X M . We allow for X° (i.e. time) dependence in all background 
fields in an assumed form for the effective action later. First we define the partition 



function, Z, and connected diagram generator, W, for our theory, using (2.1) and 
associated results. 



Z = J V[X"\ exp - {S x + S s ) = exp(-W), 



(6.2) 



where we have added a source term Ss which interacts with the fields X M (where, 



comparing the notation deployed in section (2.1), J u = y/jR^rj^V^). 

Ss= ~LS d2 ^ R{2) ^ v " xu - 



Following from (2.3) the classical fields X£ are defined by: 

SW 1 



6W 



(6.3) 



(6.4) 



Hence as in (2.5) taking the second functional derivative gives (where the notation 



< > is as used in (2.6), and £ and ( are worldsheet parameter variables of 

integration): 



5 2 W 



XZ(t)X*(Q-<X"(Z)X'>{Q> 



(6.5) 



As in (2.9) we then introduce the Legendre transform of W: 

T = W- J (ft-y/jR^VX^. 



(6.6) 
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Using the result (2.12) we take the first derivative of T with respect to the classical 



field, and then the second as in (2.13|). 

1 



5T 



1 



5 2 T 



5 2 W 



(6.7) 



We now deploy the methods described in section (2.3.3) to derive an exact equation 



for the evolution of the effective action, T of our theory with the parameter A. We 
make several comments on this method first. 

• As with the use of the parameter / in the treatment of the axion, we use 
the parameter A as a tool in setting up evolution equations for our effective 
action, following which we can examine solutions with no variation of the 
action with our chosen parameter, A. 



Considering the bare action (6.1) the progression A — > oo corresponds to 



a' — > and thus approaches the classical limit of our theory where the 
kinetic term dominates and dilaton mediated interactions are negligible. 

As A — > ^7 we approach the full quantum theory as the $(X°) interactions 
become important. 



Where a dot signifies a derivative with respect to A and noting the result W = T 



as derived in (2.52) we have an exact equation for the evolution of T. 



r = ^y^V7(7 a V + £ a V) <d a X"d b X v > 



4ir 



-Tr 



(7 a V + ^ a V) 



d_d_ 
WW 



In this we define the trace as: 



Tr[ ] = J d^d\Vlfk^~0[- 



(6. 



(6.9) 
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Using the results ( |6.4[ ) and (6.7) we have an exact evolution equation. 



+— Tr 

47T 



d d 



5 2 T 



(6.10) 



As with our ansatz for the assumed form of the Legendre action of the axion in (4.8 ) 
we now also assume a form for the Legendre effective action of the bosonic string 
so we may equate two separate, A-dependent expressions for T. In the ansatz for 
the bosonic string's effective action we assume a form for the gravitational tensor 
which allows for curvature but varies only in the X° time dimension, with D the 
total spacetime dimensionality, which we will later assume to be four as it becomes 
a free parameter of the theory. In this we have an isotropic universe in terms of 
curvature. 

9lw = diag (k{X°), niX ), r D ^(X )) . (6.11) 

We denote the antisymmetric tensor B^ v again varying only with the X° time 



dimension. Using (6.11) the assumed form (as opposed to the exact expression 



given in eq. (6.8)) for the effective action is as follows. 



1 

Air 



D-l 

d 2 £V7 { r b <x°)d a x°d b x° + 7 afe n{x Q )d a x*d h x l 

+ e ab B„ u {X°)d a X^d b X u + R^<f)(X ) 



(6.12) 



where K,Ti,<p are A-dependent functions of X°. We thus have established two 



expressions for the Legendre effective action T for our theory in (6.12), an assumed 



form, and in (6.10) an exact non-perturbative differential equation. We utilise 



these, along with the constraints of conformal invariance expressed through (5.57) 
to search for meaningful solutions. 



6.1.3 Detailed Calculations 



We wish to evaluate our exact evolution equation (6.10) by substituting on the 
assumed form for T in eq. (6.12). We take a flat worldsheet metric (the fiducial 
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gauge) 7 = 5 and configurations of the classical fields constant in spacetime, 
-^cKO = x>1 - ^ or brevity we now drop the cl subscript denoting a classical field so 
that X^ = X^ = x^. The second functional derivatives of T in (6.10) are ( with 
no summation on i) and A is the worldsheet Laplacian operator. 

S 2 T 



5X°(()5X°(0 
8 2 T 

5 2 T 



|XM=0 



Tl A5 2 (e-C) 



5 2 (£-C), 



2^ 

z^j 



^(C)^(O = u (6 ' 13) 

where a prime denotes a derivative with respect to x$. We are thus considering 
a constant field configuration x M where the off diagonal spacetime terms play no 
role. As such, the antisymmetric tensor B^ v does not play a role in the evolution 
equation for T\, as its contribution is proportional to: 



e ab d a d b 5^(^-0 = 0. 



(6.14) 



With regards to conversion to momentum space on the worldsheet, we note that 
our closed string model is described by a spherical worldsheet with non zero 
curvature which we neglect for the purposes of taking Fourier transforms required 



for momentum space formulation using F(p) = J d 2 ^e~ ip ^ F(£) (in eq (6.15) below 
£, £ become p, q). For the purposes of evaluating eqs. (6.13) we consider a small 
but non zero R( 2 \ In phase space on the worldsheet, (6.13) read as: 

5 2 T 



SX°{p)5X°(q) 
6 2 T 



Air 



2Kp 2 + R^ ( p")S 2 (p + qy i 



n v ,2 



S 2 (p + q). 



(6.15) 



5X i (p)5X l (q) 2tt 

The area of a sphere with curvature scalar is 8tc/R( 2 \ so with the constant 
configuration X^ 



x^ we have 



T = 20 A (x o ), 



(6.16) 



where the A in signifies a A dependence. The term containing the trace in (6.10) 
is computed as follows. The results are in phase space using the results in eqs. 
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(6.15) and where [ ] in eq. (6.17) represents 

1 



a_a_ ( 8 2 r 

HUVdpdq \8X»(p)5X^(q) 



Air 



-Tr 



d 2 p 

(2tt) 2 \2k P 2 + RM(f) 



1^1 



2 <r-f Ti 

i=l 



8vr 



A 2 



+ 



2k 2 



In 1 + 



2A 2 /t 



To calculate this trace, we rely on the result: 

5 2 {p = 0) = world-sheet area 
We have used the definition of the trace as Tr[ ] = 



•57T 



(6.17) 



(6.18) 



Jd 2 pd 2 q^f q 5 2 {p -q)[}. 
In this, A is the worldsheet high energy cutoff used to regularise the theory. We 
note here that A does not appear in any of our final results, nor do we attempt to 
evaluate it. The evolution equation for T, expressed in effect as an evolution of 



is finally obtained by putting together results (6.8), (6.16) and (6.17): 

D-l 



A 2 



2RW 



2A 2 k 

WW' 



(6.19) 



We would like to find solutions independent of A (i.e. where = 0) and 
configurations which preserve conformal invariance. A fixed point solution is 
possible for k = F<fi" where F is a constant: 



A 2 



2RW 



' 1 + K §^ 



if the following condition is also satisfied: 

RW 



K 



i=i 



-1 + 



T; 



2A 2 F 



In 1 + 



4F 



2A 2 F 

~rW 



In 1 + 



2A 2 F 



constant. 



(6.20) 



(6.21) 



For large A, this constant is negative, representing a Minkowski signature in target 
spacetime, a desired configuration. We have then the following conditions to 
satisfy, in order to have an a -independent solution: 

k(X°) oc </>"(X°) 
1 1 



D-l 



We now test for and apply the further constraints of conformal invariance. 



(6.22) 
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6.1.4 Conformal Invariance 



6.1.4.1 Introduction 



We will base our demonstration of conformal invariance in the vanishing of the 



beta functions in (5.57) for our theory. A key point of our derivation resulting 



in (6.19) is that it originates from a non-perturbative expression. We thus want 



to demonstrate conformal invariance to all orders in a'. This contrasts with the 



common approach in string cosmology where the expressions in (5.57) are usually 



considered to first or at most second order in a'. We will do this by demonstrating 
first and second order a! conformal invariance and by induction showing this, under 
certain assumptions, can be applied to all orders. 



To first order in a , the beta functions (5.57) for the bosonic world-sheet cr-model 



theory in graviton, antisymmetric tensor and dilaton backgrounds are, [HE] and 
[109] : 



fJ,U 

B(l) 



Ruu + 2VuV„c 



1 



pa 



1) 



V fi H plu , + d"(f>H mv 

V 2 + d"4>d P 4> - —H, vp H^ p 



2 

D - 26 1_, . . I 



6a' 



24 



(6.23) 



We note that the terms in (6.23) are expressed in terms of R^ u and H PfMU rather 
than the tensors g> M „ and B^. We now consider an ansatz for H ppu . 



H ppv = w ■ (X°) m £ P pua(X°)d a h, 



where 



h = h x X x + h 2 X 2 + h 3 X 3 



(6.24) 



(6.25) 



j=i 



(Note: the dot in uq ■ (X°) m indicates multiplication of Uq by (X°) m whereas 
£pp,va{.X Q ) indicates an (X°) dependence in e ppua ). We consider the spacetime 
metric function of the X° coordinate: 



g pu = diag( K (X°), n(X°), r,(X°), r k (X )). 



(6.26) 
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In (6.24), e P n Ua is the Levi-Civita tensor in four dimensions which has some X° 



dependence. In the linear sigma string cosmology model, [96J, the linear function 



h plays the role of the string-axion (also noted in (5.70)), (with u(X°) = exp(2$)) 



We initially limit our choice of u ■ (X°) m e pill/a d a h such that: 



h n X 3 



co -(X 



0\m 



(6.27) 



where where ho is a constant and thus the antisymmetric tensor field strength is 
anisotropic as the X 3 spatial direction is preferred. For this by definition (eom) 2 = 
det(g Mi ,(X )) and thus we have for the field strength: 



H 0i j = ^KTiTjT k 



Tk{X») 



(6.28) 



We have taken an anisotropic configuration of the antisymmetric tensor field 



strength. The more general isotropic case with the full form of (6.24) will be 



considered later in section (6.3). 



6.1.4.2 Structure of the First Order Beta Functions 



The first order beta functions we will use are noted in (6.23). We assume a X° 
dependence for 0, k and 7$ given our constraints on these parameters in (6.22), 
where O an d k are constants. Prime denotes derivative with respect to X°. 

0'(X°) 
t,(X°) 



(XO)n' 

/to 

( X 0)n+1 

-kq 
(X°) n * ' 



(6.29) 



We now show that n = 1 to ensure the functions in (6.23) are homogeneous in 
powers of X° given the properties of the spacetime (i.e. the form of the geometric 
parameters and operators R^, R, V M etc). This is obviously essential to ensure 
any meaning to the summations of the terms in the beta functions. For the 
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spacetime metric g^ u (X°) = diag(/t(X°), ^(X ), r £) _ 1 (X )), the non-vanishing 
components of the Christoffel symbols are (without summation over the space 
index i) 



1 Of 



r ,, = - 



Ta 



II 

2n 1 00 2ft' " 11 2k 

so that the non vanishing covariant derivatives of the dilaton are 



(6.30) 



V O V O = 

2ft 

2ft 

The non- vanishing components of the Ricci and Riemann tensors are 



D l 



R... 3 



2n 



T-T ■ 

i J 

4ftT,- 



4ftT, ; 



2n 



R 



00 



D-l 

E 

2=1 

I 

4ft 



4ftr,- 



ft 



In 



2n 



+ 



2n 



(6.31) 



(6.32) 



Together with the condition ft oc 0", with power laws for the metric components, 
Rqo is homogeneous, in terms of powers of X°, to VqVq0 only if 



0(X°) 
ft(X°) 



o ln(X°) 

^o 
(X ) 2 ' 



(6.33) 



and thus we have n = 1 in (6.29). Given (6.28) the non- vanishing components 



of H P11V are (such that p ^ fi ^ u, and listing only one permutation here for 
brevity) : 



Hqu — 



= e 0123 d 3 h = -K h (X )^- 



-m-na-2)/2 



(6.34) 
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From our condition h = h X 3 we have singled out the X 3 direction, thus by 
symmetry n\ = (although we leave these denoted as n\ and n 2 for the time 
being). We may now evaluate the first order in a' beta functions in terms of nx, 
n 2 , n 3 , the time dimension X° and constants h , kq and </>o with a view to finding 



beta functions homogeneous in X°. The equations in (6.23) may be evaluated, 
(without summation on the index i). 

3 



A 



5(1) 
00 



1 



4(X 



0^2 



2h 2 {X 



0\n 3 



E 







5(1) 



n,; 



5(1) 
33 



4(X )^ 

4(X°)™3 
h, 



2hi(x°r+^ + Yl 



71, 



i=i 



1,2 



i=i 



/j^0\(n3-ni-"2)/2 



n 3 



_ii + io_ 

3a' 4k 



i=i 



/> 2 
''-o 

4«n 



(X u ) 



0\n 3 



(6.35) 



In order for these beta functions to be homogeneous, it is necessary to have n 3 = 0. 
The second equation of (6.22) can then be satisfied with ri\ = n 2 = 2. The non- 
vanishing components of H pflL/ are: 



K h 



012 



(X ) 3 

Meanwhile the non-vanishing components of the antisymmetric tensor B pu are 



(6.36) 



B 



12 



-B 



2K h 



21 = (x ) (n3 " ni - n2)/2 . (6.37) 

n 3 - nx - n 2 

Thus we have a configuration for the beta functions in which they are homogeneous 
to first order in a': 



B^(X°) 



diag 



(X ) 2 ' (X ) 2 ' (X ) 2 ' 



(^1^2 - S^Sui) 2 (xo)2 
0oln(X°). 



(6.38) 
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We note here that the components of the matrix g^ u (X°) are dimensionless, with 
no a constant and the time-variable X° rescaled so as to be dimensionless. We now 
check homogeneity in the beta functions to the next order in a'. 



6.1.4.3 Structure of the Second Order Beta Functions 



The expressions in (6.23) may be extended to the next order in a'. The precise 
expressions for the second order beta functions depend on the renormalisation 
scheme employed, [109J. However, the X° dependence will remain unaffected by 
constants and/or signs in front of the terms. We consider a representative list 
of possible terms (expressed by contractions of the Ricci or Riemann tensors, the 
dilaton field and the antisymmetric tensor field strength). We rely on |109j for a 
list of such terms in second order in a'. We note that while the list of second order 
beta function terms is not the complete list from [109J, we have taken all possible 
combinations of the parameters R pi/ , H^x, R^ U a P and and derivatives thereof, 
within this list. Alternative contractions of the combinations of the parameters 
simply involve the metric g^ v or its inverse and do not change the power of X° for 
that term. For example, a term of the form: 

Roa/3yRo al31 (6.39) 

is of the order (X ) -2 . However this simply a component of the tensor: 

R m ^Rf\ (6.40) 

and the square of this is: 

Rfia/j-yRp 01131 R P 5e\Rif eX = Rpa^yR^ 1 9^ RujSexRj^ , (6.41) 

which is of order (x°)~ 2+2 ~ 2 = (X°)~ 2 as the inverse metric g puJ is of order (X ) 2 . 
Thus in taking higher powers of the terms one must contract the indices carefully 
with the metric or its inverse in order to preserve covariance. We thus conclude 
that all powers of the terms containing indices are of order (X°)~ 2 . The terms 
from |1(J9] along with the powers of X° they imply based on our solution (6.38) 
are computed below. 
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. for/^ and/^ = l,2: 



p p _ " / v0\-2 



4/? 2 

rtaPpa tt tt _ 
-K tlOaP-tlOpa — ) 



• for f3^ 3 : All the contributions vanish. 



• for /3f 2 : 



it r H ia/ 3H ipa — (A J 



Ah A 

H p<Tl H^H%W m = -^°(X )- 2 



V H afh V H°*> = ^(X )- 2 



k 

K 



(6.42) 



(6.43) 



• for 
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{H aM H°^)' 



36^ 



6 



2 



T>afipa TT TT A _ U 

ti n a B\tL n 



20h% 



V x H a ^V x H^ = ^ (6.44) 



Thus to this order in a' the beta functions are homogeneous within each beta 
function, in powers of X°. We are thus are satisfied that we have shown 
homogeneity of the beta functions in X° to second order in a'. While we have 
not done a similar exercise for third order or higher beta functions, we note we 
have taken all possible combinations of the parameters R pu , H^x, Rp Vap which 
produce the correct tensor form of the beta function concerned. Higher powers 
of these combinations will need to be contracted with the metric or its inverse 
which will produce a similar dependence on X°. The leaves the scalar field and 
the scalar curvature R as potentially producing terms with differing powers of X°. 
|109] notes that the dilaton field only occurs in the beta functions as a derivative 
(i.e. V„0). R, which is constant for our configuration, does not appear in any of 
the first or second order beta terms and we assume it is absent from higher order 
terms. Thus, while we do not prove third order or higher terms are homogeneous 
in X° we work with this assumption. 



6.1.4.4 Field Re-definition and Conformal Invariance 

Thus the beta functions using our solution are homogeneous to first and second 
order in a'. We extrapolate this based on the fact that to all orders, by design, 
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beta functions are homogeneous in the spacetime indices. Given our consistent 



configuration (6.38), we expect this to be true for all orders as we have shown that 



all possible contractions of the spacetime indices in the beta functions produce 
terms homogeneous in X°. We can write down a general power dependence for 
the beta functions in terms of X° in terms of powers of a' and X°. 



A 9 



00 



fix 



01s 



fa 



P 4 



0\2 



n=0 



En 



(I»)^ (n(a) "'" (X° 

oo 
n=0 

j oo 

— y^Vnia'T = E 4 , 



(6.45) 



n=0 



where the coefficients £, n , (m Xm $n, Vn are independent of a . The terms -£^,1,2,3,4 
are constants and it is seen that ^ and /3| 3 are independent of X°. We need 
now to make the step from the beta functions being homogeneous to all orders, to 



conformal invariance. We do this using local field redefinitions. In (5.58) we noted 



that the vanishing of the beta functions is not affected by a transformation which 
consists of a re-definition of the beta functions. 



^ ^ fa = ^ + 



where 5f3 i is defined as (<f - g*) in eq. (|6.49|) as (3 i -f3 i 



(6.46) 



notes 



that this property is due to target space diff invariance, a symmetry essential in 
general relativity, and denotes the transformed f3 l as the Weyl anomaly coefficients. 
This point is developed explicitly for use with a graviton and dilaton background 
in |lU4j with reference to work conducted related to second order beta functions in 
[109] . In [106] this is extended to beta functions in the action under consideration 



in this thesis, i.e. (6.12). In [106J it is noted that the target space physics in 
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the one loop model (that is one order in a') is unchanged under the following 
transformations. 

9nv = 9p,v + ag^u {a\R + a 2 d P 4>d p (f) + a 3 V 2 + a 4 H p ^ u H pfny ) 
4- = B^ + a (b 1 V p H ppv + b 2 d p <pH p ^) 

= + a (c x R + c 2 d p (f>d p (f> + c 3 V 2 + c 4 H p ^H ppu ) , (6.47) 

where (ai, £>i, ci, ...) are dimensionless, unspecified parameters. These 
transformations do not change the X° dependence of the fields in our configuration 



given in eq. (6.38 ) as can be seen by computing the relationship between the initial 



and transformed fields. 

B pv = A 2 B pv 

= + ^3, (6.48) 

where Ai,A 2 ,A 3 are constants. However as noted in |lU9j . the transformations 



(6.47) do change the beta functions as follows: 



A = A + iff ' 9 3 ){~ - l J . );iJ (,</' - .</') • (6-49) 



In the above (6.49) the factors g l represent the fields g^, B^ u and <p in turn (with 
Pi the associated beta functions), and the j index notation represents integration 
over X°. One may then fine tune the parameters (oi, 6j, ci, ...) (which 
are arbitrary) to make the transformed beta functions Pi vanish. We note that 
the terms Pi remain homogeneous in X° and homogeneous to the original beta 
functions pi as no operations raising or lowering powers of X° have been conducted 



in (6.49). [106J notes that given the homogeneity of the beta functions to all 



orders and the ability to arbitrarily rescale the fields, the following conditions are 
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necessary to ensure vanishing of the $i functions to second order in a'. 

E 
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oo 



33 







Ex 



(X 



0\2 







E 2 = 
E 3 



E 4 = 0, 



(6.50) 



where the constants £0,1,2,3,4 are linear functions of the parameters (01, 61, Ci, .. 



The equations (6.50) consist of five parameters £0,1,2,3,4 which are in turn linear 



functions of the arbitrary constants a^, b, L and q in (6.47). There are ten such 



constants and six linear equations and thus solutions may always be found. By 
induction, |106] argues that this can be extended to all orders and thus we conclude 



that our non-perturbative evolution equation (6.20) with the conditions (6.22) can 



be shown to have non vanishing beta functions to all orders, therefore conformally 
invariant and of interest in the string cosmology context. 



6.2 Cosmology 



We have demonstrated that there is a viable (conformally invariant) solution to 
the non-perturbative effective action on the world-sheet. We can now utilise the 
results in target space and suggest cosmological properties. We first consider what 
kind of physical universe our solution represents and then couple it to a real-world 
field, the electromagnetic field. 
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6.2.1 The Dynamics of the Universe 



In the comments surrounding eq. (5.61) it was noted that the Einstein frame 



for a worldsheet effective action and ensuing results is a convenient one in terms 
of describing the real world picture of the universe's evolution. The generalised 



expression for the Einstein frame line element is quoted in eq. (5.60). We have 
in D = 4 using, our configuration (6.38) (with our scale factors = a(t) now 



generalised to allow variation with spatial direction): 



g E (t)dx^dx l 



dt 2 - a\(t)(dx 1 ) 2 - a\{t){dx 2 ) 2 
exp {— 20(x )} g fJjU (x )dx tl dx l/ 
' dx°\ 2 ( dx 1 



(£°)- 2< ^0 



X" 



a\[t)(dx 



dx 1 



3\2 



(dx 



3\2 



X 



(.6.51) 



and 



Here, the x^ are the zero modes of X 11 (here the classical field, i.e. X c 
di(t) are the scale factors in the different space directions. In this, t is the Einstein 
frame time dimension (or cosmic time). We have then 



dt 
dx^ 



K {X 



0\-l- 



such that 



100 



KQ 'x )-^. 



(6.52) 



(6.53) 



(We note here we only consider the solution such that ^ > as we require a 
cosmic time which does not diminish). The scale factors are: 



ai (t) = a 2 (t) = a t 1+1/(t>0 
«3(0 = a t 



(6.54) 



where is a constant. The linear sigma string cosmology model for which 
B„.„ = 0, (5.63) shows linear expansion in all spatial directions and a flat Minkowski 



' [IV 



metric. Our model shows linear expansion in the direction for which we fixed the 
antisymmetric tensor field strength in the X 3 direction in eq. (6.27) with dilaton 
(0o) dependent expansion in the x\ and x 2 directions for which H piiv is not fixed. 
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We also note than when (j) = —1 and when B„ u = we revert to this currently 



observed isotropic, Minkowski universe, (5.63). We will couple this configuration 



to an electromagnetic field in the next section. It was discussed in section (5.2.3.2) 
that our configuration as above gives rise to a power law dependence by the scale 
factors a,- with cosmic time, in the Einstein frame. We noted that in the a-model 



frame as described by eq. (5.59) a de Sitter, inflationary universe is possible, 
but the introduction of the exp(— 20(x )) factor in front of the metric in the 
Einstein frame removes, following Big Bang, the possibility for an exponential 
relationship representing inflation, i.e. a(t) ~ exp(Ht) as our configuration for 
<p involves a logarithm (0(X°) = 0ohi(X )). [96J notes the importance of a non 
trivial dilaton field in accounting for inflationary and other expanding universes in 



string cosmology, as seen in our configuration for the scale factors in eq. (6.54). 
[9"6] also notes Pre Big Bang (PBB) models in which the singularity is replaced by 
a dilaton-potential barrier such that before and after the dilaton becomes strongly 
coupled. The weakly coupled stage prior to big bang has the solution (not derived 
here), [96]: 

= -(1 - EiCti) ln(-t), = 1, t < 0. (6.55) 
Here a,- are the scale factors, t is cosmic time, and d> is the dilaton field. We 



note this is similar to our configuration of the dilaton field in eq. (6.38). [96] 
qualitatively argues that at cosmic time t < in this PBB scenario and after 
crossing the dilaton-potential barrier, one expects an inflationary phase followed 
by a graceful and ensuing normal expansion to modern times. While we will not 
analyse the PBB scenario in detail, we note that our configuration may encompass 
an inflationary phase within the context of PBB theory. 

It is noted in passing that the treatment of the bosonic string in this thesis including 
the tachyon is conducted in |82j and the configuration in this study allows for a de 
Sitter universe with the condition 20+ T = 0, where T is the tachyon field. We also 
acknowledge that there is an issue of how to emerge from anisotropic expanding 



universe described by the conditions (6.54) to a isotropic Minkowski model. While 



we do not cover this point in detail, it is possible that a time dependence on h in 
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(6.27) could address this issue, i.e. the string axion field, which relaxes to zero in 



cosmic time: h(t) — > 0. 



6.2.2 Optical Anisotropy 



We described recent research in section (5.2.4) which dealt with how the coupling 



of a pseudo-scalar field such as the string axion to the photon may result in optical 
anisotropy on a universal scale. We wish to test this with our configuration. 



First we reformulate our field configurations in (6.38) by considering the (fio = —1 



solution such that the universe exhibits linear expansion in the X 3 direction but 
static in the other spatial directions. We rescale xo — > so that t = x$ 

where t is the Einstein frame cosmic time. We take a' = A = 1 and the solution 
becomes: 

g^{t) = diag(l, -1, -1, -t 2 ) 

-#012 (* J — 

<f>(t) = -Int. (6.56) 

Motivated by our desire to couple the string axion h to a photon represented by a 
gauge field A p we introduce a modified field strength, (in a similar to manner to 
that done in |115j ). H ptlv : 

H Pfa/ = H pilv + — AfrFp,,], (6.57) 

where M is a parameter with mass dimensions dependent on the string model and 
the electromagnetic tensor. In this solution we assume the string lives in D = 4 
spacetime dimensions and we refer again to our non-perturbative demonstration 
of conformal invariance which does not limit the value of D. a' is the only 
dimensionful parameter in the bare string action and since it has units y/l we 
regard M cx a/o 7 and as the string mass scale. We could like to couple the photon 



field strength F pu to the string background fields defined by the solution in (6.56) 
and hence identify any anisotropies in F^ in the low energy effective action. We 
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note that while our configuration in (6.56) is based on a non-perturbative effective 
action equation, we are not able to formulate such a non-perturbative approach 
for an effective action involving the background fields represented in ( 6.1[ ). We 
point out that our results have been designed to apply to a late (post inflationary) 
era in the universe and we thus justify at this stage using a one loop model in 
the Einstein frame, (zeroth order in a'), [88] and used in [107j. We note that this 
action, quoted here without derivation, is a perturbative effective action. While our 



model and configuration (6.38) is based on non-perturbative methods, we assume 
this configuration to be valid for our purposes in weak field configurations, R, H 



and F, 



flU- 



S, 



eff 



- / d*xy/\g\ 



AA 



e^ + R- 2d L 



-44 _ _ 

n hpv _ p pv 

11 PV-v 11 ^ 1 \iv L 



(6.58) 

As discussed surrounding our eq. (5.61) from |HS], there is a relative change of 
sign between the curvature and dilaton kinetic term, when transforming from the 



a-model frame to the physical Einstein frame. Inserting our conditions (6.56) into 
this and considering only the electromagnetic sector we have: 



S 



EM 



dtdx | - jF^F*" + etA [0 F 12] j 



(6.59) 



where e 



This part of the action as it is is valid only for a specific 
gauge transformation — > + d^ip as this transformation will also impact 

^ F for 



the antisymmetric tensor Bn V which must transform as B^v — > B 



the action (6.58) to remain invariant. The term etA\^Fi2\ arises in order to make 
the action (6.59) invariant under the transformation of B^ v . The factor t 2 in front 
of F^ V F^ U arises from our configuration (j)(t) = — Int. We consider this effective 
action in terms of the equations of motion of the electromagnetic field, that is: 



dC 



d{d u A, 



dC 

dA,, 



0. 



(6.60) 



6.2 Cosmology 



134 



where £ 



F^F"" + etA l0 F 



[0-^12]- 



VB 


— 









V-E 




t 






V x K 

V /\ l J 




-d t B 






(V x B)j 




d t E 1 + 




2e „ 2e , 

T ^ + ^ 2 , 


(V x B) 2 




d t E 2 + 


I £2 + 

—E 3 
t 3 




(V x B) 3 




d t E 3 + 




E=(E! 


E 2 


E 3 ), 


B = 


(Bi, B 2 , B 3 ) 


E = {Ex 


E2,t 2 E 3 ), 


B = 


(t 2 B 1 ,t 2 B 2 ,B 3 



(6.61) 



where we define 



(6.62) 

Here Ei, Bj are the components of the electric and magnetic fields in the 
"conventional" electromagnetic Lagrangian £ = -\F^F» V . We note the presence 
of the gauge fields A\ and A 2 in (6.61 ) is a result of the term etA\ Fi 2 ] in the action 
(6.59). We notice from (6.61) that the the X 3 direction is anisotropic, as expected 
from our choice of pseudo-scalar axion h = h^X 3 . We hence look at the time (t) 
varying (represented by the overdot) propagation of an electromagnetic wave in 
the X 3 direction arising from (6.61). 



» 3. 2e- 
E+ -E + — Ej 

t t 



0. 



(6.63) 



where we have truncated to include only powers of 1/t, justified by the relatively 
large cosmological time scales we are concerned with. Here E = (Ei,E 2 ,0) and 
Ej_ = (— E 2 , Ei, 0). This has a solution, by inspection, of 

Kf{x 3 ^ 



E 



t 2 



exp (— 2ie\nt) 



(6.64) 



where K is a constant complex number and (x 1 , x 2 ) are in the complex plane. f(x 3 ) 
is an arbitrary function of the x 3 spatial coordinate. From the V x B expressions 



in (6.61) we also have: 



V x B 



Kf(x 3 
t 3 



exp (— 2zelni) 



(6.65) 
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We take B and E as representing an electromagnetic wave lying in the (xi,^) 
plane and propagating in the 23 optic axis. If we take f(x 3 ) = cos(x 3 ) we have 
solutions for the B and E fields: 

,3\ 



E = Kexp(-2ielnt) 
B = zKexp(-2ielnt) 



cos^x 

t 2 

sin(x 3 ) 
t 3 ' 



(6.66) 



Using this we obtain for the term 



-^-F fJjl/ F f11 ' in the effective action: 



-20 



t 1 



E • El + IB • B| 



K 

2t 2 ' 



(6.67) 



4 p 2 

In arguments that will not be detailed here, |107j . it is conjectured that the 



configuration (6.56) coupled to the gauge field A„ is conformally invariant. This is 



done by applying again the argument of homogeneity combined with the potential 
to redefine the fields to achieve the vanishing of the beta functions. As seen from 
eqs. (6.83), the effect of the dilaton (i.e. the factor e~ 2< ^ = t~ 2 in front of F^F^) 



is to damp the amplitude of the electric field: 



|E| 



|K| 



(6.68) 



Another aspect of eqs. (6.66) of interest is the change in polarisation which is 



generated by the term yEj_ in the equation of motion of the electric field where 
the direction of the electric field (in the xi,X2 plane) rotates with the angle: 

AK h 



A(t) 



arg 



arg 



4e hit 



M 



Int. 



(6.69) 



This contrasts with the relationship postulated by |118j in (5.80) which shows a 
simple proportionality between the change in polarisation and cosmic time (in this 
source denoted by rj and r respectively). 



At this stage we make a comment on the term ^e 2 ^ present in our effective 



spacetime action in the Einstein frame, eq. (6.58) which derives from the 
D = 4 case in the general expression in (5.61). We note this implies a negative 



cosmological constant relative to the Einstein curvature term. We do not explore 
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a solution for this and note that the phenomenological value of our model may 
be limited if this negative term cannot be compensated for by other terms in the 
effective action. 



6.3 Generalisation of the Antisymmetric Tensor 



6.3.1 Introduction 



In the last section we chose an anisotropic configuration of the antisymmetric tensor 



field strength, H Pf±u , using constraints (6.27). We now return to the generalised 
form of ( |6\24| : 

H pilv = u ■ (X°) m e ppua (X°)d°h, (6.70) 
Previously we found that R 00 is homogeneous to Vo Vo0 only if 

<f>(X°) 



Mn(X°) 

Kq 



(6.71) 



and thus a requirement for homogeneous beta functions /3q and ffi as denoted in 



(6.23). As in (6.29) we assume a form for the components of the metric Ti(X° 

-Kq 



(6.72) 



We now set up a more general ansatz for H pflu , using similar notation as in (6.28) 



H oij = -K h k (X Y m ~ 1+( - n ^ +n ^e i:jk . (6.73) 



where we have taken u>q ■ (X°) m 
invariance of this configuration. 



(X°) m . We again consider the conformal 



6.3.2 Conformal Invariance 

As in the previous section, we subject our configuration and model to the two- 
step test of conformal invariance. First we ensure that the configuration can be 
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homogeneous to all orders in a 1 in the beta functions in the variable X°. Then 
we redefine the fields as before to make the beta functions vanish. We consider 



the first order beta functions. We use the configuration of Hq^ in (6.73) in the 



beta functions in (6.23), (without summation on the index i, where i,j,k all run 
between 1 to 3). 

3 . 3 



o9(l) 
^00 



4(X 



0\2m+n 4 -2 



i=l 



i=l 



9(1) 



n, 



4(X 



0)rn 



11; 



i=l 



+ 



0\2m+n,— ri; 



■</(!) 



s(1) _ u; 2 (^)(*°) 2m 



z^j = l,2,3 







B(i) 



i) 



3a' 4k 



rii rij rik m 
T + T ~~ T ~~ 2" ~" 



ft; 



i=l 



E 

i=l 



1.2 2 



(6.74) 



For these to be homogeneous, the constraint must apply: 2m + rij = for those 
i which /ij ^ 0. We are interested in looking at the most general, isotropic 
configuration so we take n\ = ri2 = by symmetry, but for the time being 
leave the notation separate as n^^- We now look at representative terms of 
the next order beta functions (again quoted from |109j ) for homogeneity (here 
k = 1,2,3). 



for f3 9 m : 



-Ro«/37-Rq 



a/37 



H pa0 H aaP H p ^ l H 1 a0 



16k,q 



1 3 

R a ^H 0af3 H 0pa = - — V h\n t n 3 (X 

4Kn 



0\2m-2+n fc 



i,j,k=l 



«0 



i,j,k=l 



1=1 



V H aPl V H a ^ = J2 h 2 k (n l + n J -n k -2m) 2 (X°) 



2 / -^0\2m+n fe -2 



i,j,k=l 
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• for /3fj (no summation on i): 



117 



k=i 



• for j3fj, with i ^ j: 



r>af3pa tt tt 
-fl- ia.fi 11 j pa 



K ° 1=1 



III 



for /3f j: with i ^ j ^ k: 



njhk 
8k 



(rii 



{m + rij -n k - 2m)(X°)-^ +n i- nk - 2m)/2 
+ rij -n k - 2m) \hl(X°y- ni - n i +3nk+6m V 2 



2k 

+fl k fl^x0yn l -n J +n k +6rn)/2 



+ -{ni - rij + n k - 2m)h k h 2 (X°Y~ n > +n J +n * +6m V 2 
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for ^: 



R a ^H apx H pa A 



36 



Am+2rii 



i=l 



1 3 



.i,J=l 



2kq 



' • , i 



£ n l (n i + n J )^(I ) 2m+nfc 
i,j,k=l 

3 



0\4m 



+ £ /i 2 [3/i 2 (X°) 2n < + 4/ i 2 (X°) n * +r ^ 



^ /,2 ( n . + n . _ nfc _ 2m ) 2 (X°) 



i,j,k=l 



In terms of homogeneity in X° the above do not provide any new constraints 



compared with those in the first order in (6.74). We thus assume, as with 



the anti-symmetric tensor-anisotropic case presented in section (6.1.4) that our 



configuration is homogeneous to all orders and as in the anisotropic case, we 
may rescale the fields to make the beta functions vanish and hence demonstrate 
conformal invariance to all orders in a'. We utilise similar methods and arguments 



to section (6.1.4.3) in assuming conformal invariance to all orders following 



examination of the second order terms. 
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6.3.3 Cosmology Discussion 



As in the eq. (6.51) the spacetime metric in the Einstein frame is related to 
our metric g^ u \ 



ds 2 = g^ v {t)dx^dx v = dt 2 - a\(t){dx 1 ) 2 - a 2 2 (t)(dx 2 ) 2 - a 2 3 (t)(dx 3 f 
= esxp{-2(f>(x )}g IJU ,(x )dafdx v . 



(6.75) 



We consider the representative case where h% = h 2 = and h 3 = h ^ 0. For 
brevity we denote n\ = n<i = n and n 3 = —2m. Using our configuration in eq. 



(6.71 ) we have: 

3 

dt 2 -J2 a2 i(t)(dx l ] 



2 = ( x *y^ Ko 



i=i 



dx°\ 2 A (dx*) 2 

X° J (xO)nt 



and as in (6.53) we have: 



I0o| 



> u )- 



(6.76) 



(6.77) 



This configuration (analogous to (6.38)) is in the string frame: 



g, v (X°) = Ko diag((X )- 2 ,-(X )-",-(X )-™,-(X°) 2m ) 
Hoi2(X°) = -K,h{X Q Y l ' n 
<j>(X°) = o ln(X°). 



(6.78) 



From eqs.(6.76) and (6.77), the scale factors in the Einstein frame read then 



ai(t) = a 2 (t) 
a 3 (t) 



a t 1 +«/(20o) 
S t l - m ^° 



(6.79) 



where a , d are constants, and the corresponding target space is a power- law 
expanding Universe. For f/012 we have set hi = h 2 = 0, /13 = h ^ and take n\ = 
712 = n, but note that similar configurations (which by symmetry will be physically 



identical) exist for H 013 and H 02 3- We now note that the (6.91) scale factors no 
longer necessarily represent a linearly expanding universe, when O — — 1> as did 



the scale factors in (6.54). In that case, the specific m = condition was imposed. 
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Here we consider the more general isotropic configuration of H ptlv with dependence 
on n, m and 4>q. We again couple to the electromagnetic field, and focus on, for 
(6.88) a configuration leading to a Minkowski geometry, i.e. m = —1, </>o = — 1 
and n = 0. We note that since n 3 does not appear in the physical configurations 
below, the constraint 2m + n.j = (for those i which hi ^ 0) does not impact the 
result and m becomes somewhat arbitrary. 

KqJi 



9, 



i in/ 
H 012 (t) 
<t>{t) 



= H 013 
-Int. 



Hi 



023 



t 3 



(6.80) 



We now couple the configuration (6.80) to the photon field via the effective action 



in (]6.58|) and the gauge field action in (|6.59|), and the following modified Maxwell's 

(6.81) 



equations are arrived at. 



S7-B = 









V-E = 






t A 


V x E = 


-d t B 


Vx5 = 


d t E + 



where E± = (—E 2 ,Ei,0) and A± = (—A 2 ,Ai,0) and as in eqs. (6.61) we 
have ignored terms of order 1/t 2 or higher. As previously, motivated by cosmic 
birefringence, we consider a plane wave in the (x\,x 2 ) plane moving in the 23 
direction. Using the complex notation E = E\ + iE 2 , we obtain the following 
equation of motion 

E - E" + = 0, (6.82) 



t 

where overdot is t derivative and prime is spatial derivative. By inspection, E 
E r cos(t 



X" 



(where E and a are complex constants), satisfies equation (6.82) 



for a 



■1 — ie/2, neglecting terms of 1/t 2 or higher. The solution is then: 
Eq 



E 



t 



exp 



-i— lnt ) cos(t — x" 



(6.83) 



As in our previous solution, the dilaton provides for a damping effect (the 1/t term 
in (6.82) while the antisymmetric tensor, via the field strength, provides a varying 
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polarisation of the electric field in the (xijXz) plane, in a similar manner to the 
result (6.69) we have for the isotropic H Pflu case: 

K h 



A(t) = |arg(E) - arg(E* 



M 



(6.84) 



where t can be seen as the time interval between emission and observation of the 
electromagnetic waves, which we note provides the same optical anisotropy effect 



as in the case where the field strength tensor was not generalised as in (6.69). 



6.4 Summary and Discussion 

The primary result of this section of the thesis is to derive a non-perturbative 
formulation for the effective action of the bosonic string, on the worldsheet. We 
then showed that a viable solution of this was conformally invariant, again, non- 
perturbatively, and therefore suitable for consideration in string cosmology. To 
achieve this we commenced with a bare string action on the worldsheet with 
background fields 77^, a^ v and <I>(X ) representing the graviton, antisymmetric 
tensor and dilaton fields respectively. Here A is a varying parameter which equals 



l/a' in the full quantum theory of (6.85), where a' is the Regge slope. 



Sx = ^LJ d ^ i X ^ b%v + daX " dbXU + Ri2) ®( X °)}- ( 6 - 85 ) 

We then applied techniques with roots in the exact renormalisation group (ERG) 



methods, outlined in section (2.3.3 ) of this thesis, to derive a form for the Legendre 



effective action associated with (2.3.3): 

+ > { ^ + ^ w4 (mSxm) } ■ (6 - 86) 

Here the overdot on T indicates a derivative with respect to A, which we use 
in an analogous manner to which the high energy cutoff is used in the exact 



renormalisation group evolution equations such as that derived in Appendix (8.1) 
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We then derive fixed point or f = solutions. We then plug into (6.86) an assumed 
form for the effective action T in which the three background fields vary only with 
the X° component of spacetime and for which the spacetime metric is given by: 

g^ u = diag (k(X°), ti(X°), t d _i(X )), where D is the spacetime dimensionality. 

At this point we note that the method employed thus far is designed to show 
that a conformally invariant configuration is possible for our non-perturbative 



expression (6.86). We acknowledge that the configuration in (6.87) is intended as 
a demonstrative solution only and for the reasons below, its applications in string 



cosmology may be limited. Given this, we find an acceptable solution to (6.86) is 
given by the following constraints where prime denotes a derivative with respect 
to X°. Here <fio and are constants and n,rii are parameters to be determined. 

<f>'(X°) 
Ti(X°) 



(I )"' 

( X 0)n+1 

(X )™* ' 



(6.87) 



We then set out to show through our assumptions that the configuration (6.87) 
is conformally invariant by showing the beta functions representing our effective 



action vanish to all orders in a'. The beta functions are given in (6.23). We do 



this by first showing that the beta functions with (6.87) are homogeneous in X° 
to first order in a'. We then show the same for selected terms from a complete list 
of the second order beta functions which make up all possible combinations of the 



field parameters. Using arguments outlined in section (6.1.4.3) we deduce, but do 
not comprehensively prove, homogeneity in X° to be true to all orders in a'. We 



find that such conformally invariant configuration of (|6.87|) is given by: 

K 



diag 



-Kq 



-Kq 



B^X ) = (5^-5^) 



(X ) 2 ' (X ) 2 ' (X ) 2 ' 
K h 



-Kq 



2(x°y 



0(X C 



Mn(X°). 



(6.88) 



In this we have assumed a form for the antisymmetric field strength such that 
Hcii2 = — / K °n°n and thereby limiting this parameter to one particular spacetime 



(X0)3 
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configuration (hence anisotropic for the time being). We then make use of a 
property of the beta functions under which the following transformation leaves 
the underlying physics invariant, including the fact that they may be found to 
vanish. 

9nv = 9nu + a'g^u (aiR + a 2 d p (j)d p (f) + a 3 V 2 + a 4 H p ^ u H ppu ) 

= + a ( Cl R + c 2 d p (f>d p (f) + c 3 V 2 + c A H ptlv H^ v ) , (6.89) 
where (aj, &i, cj, ...) are dimensionless, arbitrary parameters. Thus we can 



choose these parameters in such as way that the beta functions in (6.89) vanish 
and hence we have shown conformal invariance to all orders in a'. This proof 
relies on the fact that we have also shown them to be self homogeneous (and 
therefore the sums of terms within each beta function may cancel) in X° for our 



non-perturbative solution to (6.87). 



In terms of cosmology in the Einstein frame, (6.88) represents a metric: 

g^ v {t)dx p dx u = dt 2 - a\(t)(dx 1 ) 2 - a 2 2 {t){dx 2 ) 2 - a 2 3 {t){dx 3 ) 2 , (6.90) 

where 

a x {t) = a 2 {t) = a t 1+1/</, ° 

a 3 (t) = a t (6.91) 

The solution shows linear expansion in the direction for which we fixed the 
antisymmetric tensor field strength in the X 3 direction in eq. (6.27) with dilaton 



(</>o) dependent expansion in the X\ and x 2 directions for which H pflll is not fixed. 
When <J)q = —1 and when B^ u = we revert to the currently observed isotropic, 
Minkowski universe. Firstly, we acknowledge that our model does not take into 
account (nor allow for) an inflationary scenario (i.e. a(t) ~ exp(Ht)) and that our 
results are concerned with post inflationary phenomena. We do however note that 
our configuration can allow for inflation in Pre Big Bang theories. Secondly, we 
acknowledge that there is an issue of how to emerge from anisotropic expanding 
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universe described by the conditions (6.54) to a isotropic Minkowski model. While 



we do not cover this point in detail, it is possible that a time dependence on h in 



(6.27) could address this issue, i.e. the string axion field, which relaxes to zero in 



cosmic time: h(t) — > 0. 



We then reformulate (6.88) to read in the Einstein frame: 



H 012 (t) 



diag(l,-l,-l,-t 2 ) 

_Koho 
t 3 

-Int. 



(6.92) 



Motivated by optical anisotropy, we insert this configuration (6.92 (into a one loop 



effective action in target space where the antisymmetric tensor has been coupled 
to the electromagnetism via a modified field strength such that H plw = H pfiu + 
j^A^F^ where M is the string mass scale and F^ v is the electromagnetic tensor. 
This results in modified Maxwell's equations which admit as a solution: 



E = Kexp(-2ielnt) 
B = zKexp(-2ielnt) 



COS I" 



t 2 



sin ix 



t 3 



(6.93) 



Here E and B represent an electromagnetic wave lying in the (xi,X2) plane and 
propagating in the X3 optic axis whose angle of polarisation changes with cosmic 
time t. 



A(t) = |arg(E) -arg(E)*| =4elnt= ^^ln/. 



M 



We then generalise the configuration for H piiv 

n mn/ — [y\. 



<-piiV 



- pfiucr 



from H 012 
(X°)d a h, 



to: 



(6.94) 



(6.95) 



and repeat the process of demonstrating conformal invariance, devising the 
Einstein metric and coupling to the electromagnetic field in a spacetime effective 
action. The result is a similarly logarithmic dependence by the polarisation on 



cosmic time as in (6.94). This phenomena is also known as cosmic birefringence as 



described in section (5.2.4) and efforts are underway to examine evidence for this in 
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the cosmic microwave background and from studies of radiation from distant extra- 
galactic objects which emit polarised radiation, [1 13j . |114j . |117j . Theoretical 
treatments similar to that conducted in this thesis result in a cosmic time 
dependence by the polarisation proportional to h'r, |115j and tan -1 - and tan -1 4, 
with the variation due to the configuration of the Einstein metric considered. 
The logarithmic dependence of our solution arises as a result of the power law 



constraints we established in (6.87), which in turn were motivated by the need to 
prove conformal invariance in a non-perturbative manner. 



Chapter 7 



Conclusion 



7.1 Full Quantization of the Axion 



In section (4.4) we described how the spinodal instability in the QCD axion's 



potential flattened it when interactions are not considered. We derived a non- 
perturbative expression for the evolution of the effective potential with the scale 



factor / as in (4.15) and found a flat-potential solution consistent with the well 



known property of convexity outlined in eq. (2.59). We then assigned a boundary 
condition to this expression such that the potential when / = A is some value U\ 
leading to: 

Physically, at energy level A, which represents the upper limit cutoff of our theory, 
we consider that the evolution of the axion field is in its earliest describable form. 
In this sense, for very small 9 (the dynamical axion parameter), the commonly used 



cosine form of the axion potential (3.62) may be equated to a double well scalar 



potential as in (4.18). We consider that such a double well potential represents 



the origin of the spontaneously broken U(1)pq symmetry responsible for the axion 



field's development. We consider that our result (7.1 ) describes the system on the 



verge of spontaneous symmetry breaking. At this point the axion has yet to evolve 
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fully and the effects of the spinodal instability serve to flatten the potential, prior 



to any interactions. We consider that the term U\ in (7.1 ) is a one loop correction 
to this. We identify it with the one loop correction to the quantized form of (4.19) 
at fi = 0, which we compute using an established result quoted in j2]. For A = /, 
the result which can be expressed in terms of the axion mass m a and the ratio of 

4 

j4 where m is an arbitrary mass scale less than A, and / is the scale factor which 
we here set equal to A. 

4 

U effJ=A ~ 4 x 10- 2 -pmi, (7.2) 

4 

where j± is a ratio less than one. We interpret this as an energy density associated 
with the axion at the earliest phase in its development, prior to significant 
interactions, and influenced by spinodal instability effects. (At later stages the 
axion acquires mass and resolves the CP problem as well as providing a candidate 
for cold dark matter). With the caveat that the ratio is not determined by 
our theory and using 10 _4 eV < m a < 10 _1 eV, [26] it could be that the result in 



eq. (7.2) it is not orders of magnitude away from required values for dark energy 



{U DE ~ (io-W) 4 , M)- 

We acknowledge the limitations in this result and treatment as follows. Firstly, 
our model and assumption of the flattening of the potential does not take into 
account interactions of the axion scalar field with other particles, in particular 



with gravity. Section (4.4.0.3) shows how interactions will impact this flattening. 
Secondly, our result contains indeterminate energy scales (e.g. j^) which lessen 
the usefulness in phenomenological comparisons. Further work in this area would 



involve an extension of section (4.4.0.3) to devise and quantize a more realistic 



model of the axion scalar field in the early universe including coupling to gravity, 
including perhaps, spacetime curvature which may be relevant at the energies and 
epochs we are considering. This could then lead to a more detailed investigation of 
the flattening effect of the axion's effective potential at high energy scales and lead 
to more justifiable allocation of values to the at-present arbitrary energy scales in 
the derivation and results. 
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7.2 The Bosonic String Axion 

The primary result of this section of the thesis is to derive a non-perturbative 
formulation for the effective action of the bosonic string including the string axion, 



on the worldsheet, (6.86). We showed that a viable solution of this (6.87)) is 
conformally invariant and therefore suitable for consideration in string cosmology. 
We achieved this demonstration of conformal invariance in a non-perturbative 
manner, by demonstrating homogeneity of the beta functions with time; combined 
with field rescaling to allow for cancellation of all terms. This contrasts with the 
conventional method in string cosmology of canceling a perturbative expression for 



the beta functions up to some order in a'. We find that a configuration for (6.87) 
is given by: 

VH A ) - aiag l (xoj2' (xop' 
'Q\ - fx x x x \ K oh 



<f>(X°) = o ln(X°). (7.3) 

In this we have assumed a form for the antisymmetric field strength such that 
^012 = — p an d thereby limiting this parameter to one particular spacetime 
configuration (hence anisotropic). In this ho is the constant field strength of the 
bosonic string axion. Motivated by optical anisotropy, we insert this configuration 



(7.3)into a one loop effective action in target space where the antisymmetric 
tensor has been coupled to the electromagnetism via a modified field strength 
such that Hp^ u = H pfiu + j^A^F^ where M is the string mass scale and F^ v is the 
electromagnetic tensor. The result is a logarithmic dependence by the polarisation 



on cosmic time as in (6.94). This phenomena is also known as cosmic birefringence 



as described in section (5.2.4). The logarithmic dependence of our solution arises 



as a result of the power law constraints we established in (6.87), which in turn 
were motivated by the need to prove conformal invariance in a non-perturbative 
manner. We make several comments on this section of the thesis with a view to 
future work. 
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Section (5.2.3.2) noted that the potential for tachyon modes arising in the 
theory can be removed by applying supersymmetry. This was not done in 
our work and is an obvious extension of the model. 



It was discussed in section (5.2.3.2) that our configuration as above gives rise 
to a power law dependence by the scale factors with cosmic time, in the 
Einstein frame. The predictive powers of the model may be incomplete as 
the resulting metric does not support an inflationary scenario other than in a 
Pre Big Bang scenario. Further work could involve relaxing our constrain of 
considering D = 4 spacetime dimensions and/or considering tachyon modes 
in our non-sypersymmetric string model as was done in [82J. 

We note the ^e 2 ^ present in our effective spacetime action in the Einstein 



frame, eq. (6.58) which derives from the D = 4 case in the general expression 
in (5.61). We note this implies a negative cosmological constant relative to 
the Einstein curvature term. We do not explore a solution for this and note 
that the phenomenological value of our model may be limited if this negative 
term cannot be compensated for by other terms in the effective action, the 
form of which we have assumed. 

We acknowledge that there is an issue of how to emerge from anisotropic ex- 



panding universe described by the conditions (6.54) to a isotropic Minkowski 



model. While we do not cover this point in detail, it is possible that a time 



dependence on h in (6.27) could address this issue, i.e. the string axion field, 



which relaxes to zero in cosmic time: h(t) — > 0. 

The next steps in this part of the thesis would be as follows. We would 
develop, as far as possible, a more rigorous demonstration of conformal 
invariance of our configuration via closer examination of the higher order 
beta functions. We could further generalise the bosonic string action to 
incorporate supersymmetry. 



Chapter 8 
Appendix 



8.1 Derivation of Exact renormalisation Group 
Equations 



The Wilsonian effective action in (2.46) can be written as: 



exp(-5 fc _ 5fc ($)) = J Difex P (-S k ^ + ip)). (8.1) 

The field (f)(p) = $(p) + <p(p) is split into two sectors, with p the momentum 
variable: 

<f)(p) = (p(p) for k — 5k < \p\ < k (the UV modes) 

(p{p) = <&(p) for \p\<k-5k (the IR modes), (8.2) 



where <p(p) is zero outside the regions specified in (8.2). Expanding Sk(& + V 9 ) 



around $ gives, (where the phase space volume within p — » A is V = J j^u)'- 



(8.3) 



In (8.3), the term is the derivative of the functional 6Sk(&) with respect to 



the function $(p). The integrals J", are defined in D dimensional Euclidean space, 
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and where dVLp is the solid angle in D dimensions. 

d D p f k dVL D p D - l dp 



k-Sk 



(27T 



k-8k 



(2n) 



D 



U) 



If there is an expression f(p 2 ) under the integral in (8.4) such that the integral 
only depends on the magnitude of the momentum (p 2 ), we can perform the dVto 
integral to give: 



k dtt D p D ~ l dp M/ 2 



-HP 2 



VL D k D ~ l 8k- 



1 



\D 



5) 



, k _ 5k (2tt) d "-(2* 
Thus with these conditions each J k contributes a factor proportional to 8k. In 



(8.3) the terms higher than ^72 J k J k fip) si^s^^ viv) are a ^ l eas ^ quadratic in ip 



and therefore the J k integrals can be evaluated to give a factor of at least (8k) 2 
in the terms. Hence we use the notation 0(8k) 2 in (8.3). By substituting the 



expansion in (8.3) into (8.1), and expanding exp(0(8k) 2 as (1 + 0(8k) 2 ) we obtain 
the following. 



exp(5 fc ($) - S k ^ Sk {$)) 



Dip exp — 



Vj k 5*(p) 



<p(p) 



1 f f 8 2 S k ($) 



[l + 0{8kf) (8.6) 



The right-hand-side of (8.6) can be evaluated by considering it as a Gaussian 
integral with and A4 as matrices, and using the following, where K and K' are 
constants. 



(11* J cZ£ fc )exp- Q&-B i j& + A&J 



x exp (-iTr(lni^) + ' . 1,) . (8.7) 
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We use the following equivalences: 

n fc 



<p(p) 
1 



we compute (8.6): 



exp(S k ($) - S k . Sk m = exp(--Tr In 



k Jk 



s 2 s k 



(g) \ 



2 " V.""" 5<S>(p)5$(q) J 
r \ 5$(p) \5$(p)5$(q) J 5<S>(q) 



+ 




In the right-hand-side of (|8.9|), we assume that the factor ^^^^ 



DyO{8kf (8.9) 
in the logarithm 



has been multiplied by a unitary constant of units k to avoid taking the logarithm 



of a dimensionful parameter. In (8.9), the trace taken over the products of the 



matrices and A4 is defined as Tr[...] = f k dp f, dq(2Tr) D 5 D (p + <?)[....]. 



We now evaluate the factor in (Js - 9 [) . The effective action can be expressed 



based on the gradient expansion approximation: 

1 



eff 



d u x -Z(^)d^d^ + U eff (^) 



(8.10) 



where U e ff(<&) is the effective potential containing the mass and interaction terms 
and Z($) contains all fluctuations involving derivatives of the field $. Considering 



5S k ($) 

my) 



d D x[U' eff (<f>(x)) + ^Z\^(x))d^d^ - Zmx^d^^ix - y). 



(8.11) 

where U' e fj(<&) is the partial derivative with respect to $(?/). Evaluating S D (x — y) 
under the integral, using integration by parts on the second term on the right- 
hand-side and taking J dPx (S(x — y)d fJ- §) — > in the integration limits, gives: 



5S k ($) 



U'efAHy)) + -Z\<$>{y))d^d^ - Zmy))d^, 



(8.12) 
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In a constant field configuration $ , the factor d^d^Q vanishes and the expression 
becomes, in momentum space: 

(27T) 

5 D (p)U' eff ($o). (8.13) 



5$(p) 



1*0 



X= J d D yexp(-ip.y)U' eff ($ ) 



The integral in (8.6) is between the values of p = k — 5k and k. The S D (p) factor 



in (8.13) indicates evaluation only at p = and hence the entire expression for 



Uo van i snes when integrated in the shell for $ = constant as does the entire 



term 



ss k ($) ( <s 2 s fc (4) \ 1 sstm] ■ 



54(p) ^54(p)54(g) J 8<S>(q) 



in (8.9), resulting in: 



«p(fc(*) - ^(*)) = exp ( -i T r(ln jgg^ 



DipO{8k) 2 . (8.14) 



Evaluating the trace in (8.14) as: 



Tr(/(p,g)) = / dp / dg(27r)^(p + g)(/(p,(/)) ) 



(8.15) 



and the delta function in (8.15): 



Tr(/)= I 5 D (0)f(p,-p) 



(8.16) 



we use the identity S D (p) = J d D xexp(ipx) and that S D (0) = J dPx = Vd where 
Vd is the volume integration over space-time in the D dimensions (assumed to be 
finite): 

Tr(f) = V D I f(p,-p). (8.17) 
Jp 



Using f{p, q) = In 



6 2 S k (t>) 
8<S>(p)S<S>(q) ■ 



(8.14) becomes: 



exp(S fc ($) - 5 fc _ 5fc ($)) = exp 



V D 5 2 S k ($) 



In 



2 J p 5$(p)5$(-p) 



+ / D<fO(Skf 
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The f v in (8.18) can be evaluated as in (8.5) as the field is quadratic in p and 

5 2 S k (<S>) 



substituting into (8.14) gives: 



' " l 5k. V ",\n 



(2tt) d 5<$>(p)5$(-p) 

D<pO{5k) 2 (8.19) 



exp(S , fc ($) - S k sk(®)) = exp ( --Q D k 



We now take the logarithm of both sides, using on the right-hand-side \n(x + e) = 
In x + e given x ^> e (or in this case 5k ^> (5k) 2 ). We then divide both sides by a 
factor of 5k and finally take the limit as 5k — > 0. This limit produces the derivative 
with respect to k on the LHS of ( 8.14[ ), that is: 

(S fc ($o)-S fc _ 5fe ($ )) 



d k S k ($ ) = lim 

(5fc-s>0 



5k 



(8.20) 



On the right-hand-side of (8.19), the 0(5k) 2 factor vanishes in this limit to give an 
exact form for the variation of the action with cutoff momentum or energy scale: 



d k S k (^ 



5 2 M$ ) 



(8.21) 



(2tt) d l5$(p)5<S>(-p) 

where here k 2 = p 2 . This is known as the Wegner-Houghton equation, derived in 
the path integral formalism, [8]. 



Considering the effective action in (8.10), we differentiate (8.12) with respect to 
$(2;) with the objective of evaluating ^(^a^-p) i n terms of the U e ff and Z($). 

5^y)5%) = [Um + \ Z "^v)) d ^ 9 ^ ~ Z'(Hy))d,d^ 

- Z^{y))d^]5 D {y -z) + Z' {<$>(y))d,{5 D (y - z))9>*. (8.22) 

Evaluating s ^^t\z) a ^ a cons ta n t field value, results in the terms containing 
c?^$ vanishing and: 

5 2 S k (<5>) 



This can be expressed in momentum space as, using —d^d^ = k 2 : 
5 2 S k ($) 



5$(p)5$(q) 



\* = K{* )+Z{* )k*\6 D (p + q). 



(8.23) 



(8.24) 
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Thus in a constant field configuration <3> = <3> , (8.21) reduces to, with a static field 
potential C/ e //($o) and a configuration £4(0). 



d k U eff (<S> ) - d k U eff (0) 



Vtnk D ~ l 



2(2vr 



In 



^($0) + Z($ )A; 2 



U' e ' ff (0) + Z(0)k< 



(8.25) 
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